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1.  INTRODUCTION:  The  blood-brain  barrier  renders  the  central  nervous  system  a  sanctuary  site  for  disease. 
In  the  era  prior  to  the  development  of  effective  systemic  therapies  this  fact  was  not  of  much  clinical  relevance  as 
most  patients  succumbed  to  visceral  metastatic  disease  before  their  brain  metastases  became  symptomatic. 
However,  as  systemic  breast  cancer  therapies  improve  and  are  able  to  successfully  control  non-Central  Nervous 
System  disease,  the  brain  is  increasingly  the  first  site  of  relapse.  Both  conventional  chemotherapeutic  agents  and 
targeted  monoclonal  antibodies  do  not  cross  the  blood-brain  barrier  (BBB)  at  concentrations  sufficient  to 
successfully  treat  metastatic  disease.  The  objective  of  the  proposed  work  is  to  use  an  active  form  of  transport  of 
therapeutics  across  the  BBB;  a  transport  mechanism  that  does  not  rely  on  passive  diffusion  or  receptor-mediated 
transcytosis.  We  hypothesize  that  monocytes/macrophages  are  actively  recruited  to  metastases  by  cytokines 
elaborated  by  the  tumors  and  that  these  cells  can  be  utilized  to  transport  therapeutics  directly  to  brain  metastases 
of  breast  cancer.  This  hypothesis  is  supported  by  a  pilot  study  we  published  that  demonstrated  delivery  of 
fluorescent  microspheres  to  brain  metastases  by  macrophages  in  a  mouse  model  [  1  ] .  Our  obj  ective  in  this  proposal 
is  to  deliver  a  therapeutic,  doxorubicin,  to  brain  metastases  using  a  nanoshell-double-stranded  DNA  (dsDNA)- 
drug  complex  loaded  within  monocytes/macrophages.  Once  present  within  the  metastasis,  we  aim  to  release  the 
drug  by  transcranial  irradiation  at  near-infrared  (NIR)  wavelengths. 

2.  KEYWORDS:  Brain  metastases,  nano  shells,  macrophage,  monocyte,  blood-brain  barrier,  photothermal 
therapy,  light  triggered  release,  Adriamycin,  doxorubicin,  metastatic  breast  cancer,  Trojan  horse 

3.  OVERALL  PROJECT  SUMMARY:  Summarize  the  progress  during  appropriate  reporting  period  (single 
annual  or  comprehensive  final).  This  section  of  the  report  shall  be  in  direct  alignment  with  respect  to  each  task 
outlined  in  the  approved  SOW. 

Task  1 

a.  Current  objectives 

The  original  proposal  suggested  the  immobilization  and  near-IR  triggered  delivery  of  lapatinib  from  gold 
nanoshells.  The  goal  was  to  load  lapatinib  onto  the  nanoshell  surface  by  intercalating  the  drug  within  of  dsDNA 
sequences  bound  to  the  nanoshell  surface.  The  dsDNA  sequences  would  have  one  strand  that  was  thiolated  and 
covalently  linked  to  the  particle  surface,  while  the  complementary  strand  was  not.  Upon  irradiation,  the 
complementary  sequence  would  be  released  from  the  particle  surface  thereby  releasing  lapatinib  locally.  We 
proposed  to  load  these  nanocomplexes  into  macrophage/monocyte  cells  that  could  then  be  used  to  deliver  them 
into  the  brain  metastases. 

In  the  first  year’s  report,  we  detailed  the  difficulties  we  encountered  with  lapatinib  solubility  and  light-triggered 
release.  As  a  potential  solution,  we  decided  to  focus  on  the  delivery  of  doxorubicin  (DOX).  DOX  is  a 
chemotherapeutic  agent,  used  in  the  conventional  treatment  of  cancers  including:  Hodgkin's  lymphoma,  myeloma, 
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lung,  ovarian,  gastric,  thyroid,  breast,  sarcoma,  and 
pediatric  cancers.  Its  mechanism  of  action  is  to 
intercalate  between  the  double  strands  of  DNA.  Our 
modified  approach  is  to  intercalate  DOX  into  a  GC  rich 
double- stranded  DNA  (dsDNA)  sequence  and  under  laser 
illumination  at  the  plasmon  resonance  of  the 
NS@DNA@DOX  complex,  dehybridize  the  dsDNA  and 
release  DOX  (Scheme  1). 

b.  Results 

In  order  to  deliver  DOX,  we  synthesized  a  nanocarrier  from  which  DOX  could  be  released  by  DNA 
dehybridization.  The  dsDNA  sequences  used  in  this  study  were  engineered  specifically  to  bind  to  the 
chemotherapeutic  DOX:  5’-HS-(CH2)6-CAA  TCA  ATA  GCT  ATC  GTT  CG-3’  and  3’-GTT  AGT  TAT  CGA 
TAG  CAA  GC-5’.  Prior  to  attachment  to  the  nanoshell  (NS)  surface  the  DNA  was  hybridized  by  mixing  equal 
moles  of  the  two  DNA  strands  with  33  mM  NaCl,  heating  the  mixture  to  100  °C  for  4  min  and  slowly  cooling  to 
room  temperature  overnight.  Nanoshells  were  functionalized  with  double-stranded  DNA  (dsDNA)  via  a  Au-S 
bond.  The  dsDNA-NS  mixture  was  incubated  overnight  with  continuous  shaking.  The  solution  was  washed 
several  times  to  remove  the  dsDNA  excess.  The  NS@dsDNA  complex  was  incubated  with  DOX  for  24  h  and 
excess  DOX  was  removed  via  centrifugation. 

Vis-NIR  spectroscopy,  transmission  electron 
microscopy  (TEM),  zeta  potential,  and  dynamic 
light  scattering  (DLS)  measurements  were 
utilized  to  provide  evidence  of  the  formation  of 
the  DOX  nanocarrier.  The  extinction  maximum 
of  bare  NS  slightly  red-shifted  after  dsDNA 
functionalization,  due  to  the  change  in  the 
dielectric  environment  around  the  NS  (Fig.  1A). 

No  changes  in  the  peak  shape  or  width  of  the 
extinction  spectra  was  observed,  indicating  that 
the  NS  @  dsDNA  nanocomplexes  do  not  aggregate 
during  the  functionalization  process.  TEM  images 
illustrate  the  thin  layer  of  DNA  formed  on  the  nanoshell  surface  (Fig.  IB).  Zeta  potential  measurements  further 
confirm  modification  of  the  NS  surface  (Figure  1C).  The  zeta  potential  of  the  particles  decreased  ~18  mV  after 
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Figure  1.  Formation  and  characterization  of  doxorubicin 
nanocarrier.  (A)  Extinction  spectra  of  NS  (black), 
NS@DNA  (green),  and  NS@DNA@DOX  (red);  (B)  TEM 
images  of  NS  (black)  and  NS  @ DNA  (green);  (C)  zeta 
potential  of  NS,  NS@DNA,  and  NS@DNA@DOX. 


NS@dsDNA@DOX  NS@ssDNA  ssDNA  +  DOX 

Scheme  1.  Near-IR  light-triggered  release  of 
doxorubicin  (DOX)  from  nanoshells  by 
dehybridization  of  dsDNA. 
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DNA  functionalization  as  a  result  of  the  negative  charges  on  the  DNA  phosphate  backbone.  The  zeta  potential 
increased  slightly  after  incubation  with  DOX  due  to  screening  of  the  DNA  charges. 


Cellular  uptake  of  the  NS@DNA  complex  was  investigated  using 
confocal  and  reflectance  (emission  and  excitation  at  640  nm) 
microscopy.  For  the  initial  experiments,  RAW  264.7  a  mouse 
macrophage  cell  line  from  an  Abelson  murine  leukemia  virus-induced 
tumor  was  utilized.  This  cell  line  has  the  advantage  of  being 
immortalized  and  is  easy  to  propagate  in  cell  culture  in  comparison  to 
the  monocytes  isolated  from  human  blood  and  utilized  previously  [1,2]. 

Briefly,  the  RAW  264.7  macrophage  cells  were  cultured  for  24  h  in  an 
8-  well  micro-chamber  slide  and  then  incubated  with  the  NS@dsDNA. 

After  20  h,  the  cells  were  rinsed  three  times  with  PBS  and  fixed  in  a  4% 
paraformaldehyde  solution.  The  cells  were  washed  twice  with  PBS  and 
incubated  with  a  Hoechst  33342  solution  to  stain  the  cell  nucleus.  After 
10  min  of  incubation  and  protected  from  light,  the  cells  were  rinsed 
twice  with  PBS.  All  images  were  collected  on  a  Nikon  Al-Rsi  Confocal 
microscope  and  the  particles  were  imaged  in  the  reflection  mode.  The 
3D  maximum  intensity  projection  image,  corresponding  to  a  slice  from 
the  middle  of  the  cell,  confirms  internalization  of  the  NS@DNA 
nanoparticles  (red  dots,  Fig.  2).  The  details  in  the  orthogonal  YZ  and 
XZ  planes  illustrate  the  location  of  the  nanoparticles  within  the  cell  boundary,  around  but  not  inside  the  nucleus. 


Figure  2.  Confocal  image  of  cellular 
uptake  of  the  nanoshell-DNA 
nanocomplex  in  a  RAW  264.7 
macrophage  cell.  Orthogonal  view 
of  the  3D  maximum  intensity 
projection  image  of  cells  incubated 
overnight  with  nanocarrier.  Nucleus 
is  stained  blue  with  Hoechst  33342. 
Nanoparticles  were  located  by 
reflection  (emission  and  extinction 
at  640  nm).  Scale  bar  =  10  pm;  60x 
oil  imersion  objective;  NA  =  1.4. 


In  response  to  the  previous  challenges  we  encountered  with  regard  to  demonstrating  light-triggered  release  of 

lapatinib,  we  decided  to  take  a  step  back  and  investigate  in  greater  detail  the  mechanism  of  light-triggered  release 

in  order  to  better  control  and  tailor  the  release  for  various  molecular  cargos.  In  this  study,  we  used  NS 

functionalized  with  dsDNA;  one  strand  was  thiolated  and  the  complementary  strand  tagged  with  fluorescein 

[5’HS-(CH2)6-TAT  GAT  CTG  TCA  CAG  CTT  GA-3’  and  5’-  ATA  CTA  GAC  AGT  GTC  GAA  CT-Fluor-3’]. 

The  dsDNA  dehybridization  process  was  investigated  by  fluorescence  spectroscopy.  It  had  previously  been 

demonstrated  that  when  a  high  concentration  of  nanoparticles  was  present,  release  of  cargo  using  a  continuous 

wave  (CW)  laser  could  be  achieved  at  temperatures  below  the  DNA  dehybridization  temperature  [3].  Based  on 

this  observation,  a  non-thermal  release  mechanism  by  hot-electrons  was  proposed.  In  this  study,  we  reduced  the 

particle  concentration  to  better  simulate  the  concentration  of  particles  that  is  realistically  achievable  within  a 

tumor.  We  found  by  heating  the  nanoparticle  solution  with  a  temperature-controlled  stage  that  thermally  induced 

dehybridization  of  our  specific  dsDNA  sequence  occurred  at  61.1  ±  2.3  °C  (Fig.  3A).  Dehybridization  of  the  same 
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sequence  using  CW  light  occurred  at  59.4  ±  5.1  °C  (Figure  3B). 

In  these  experiments,  we  observed  that  release  does  not  occur 
until  the  laser  power  is  high  enough  to  increase  the  bulk 
temperature  above  the  dehybridization  temperature.  This  is 
evidence  that  the  dsDNA  dehybridization  occurs  as  a  result  of 
the  heating  from  the  nanoparticles  as  opposed  to  hot-electrons, 
which  would  be  linearly  dependent  on  the  laser  power.  The  high 
nanoparticle  concentrations  used  in  earlier  experiments  resulted 
in  optical  light  trapping  due  to  multiple  scattering  of  the 
nanoparticles,  allowing  the  local  temperature  to  increase  and 
dehybridize  the  DNA  before  a  significant  rise  in  the  global/bulk 
temperature  was  measured.  Unfortunately,  that  nanoparticle 
concentration  is  not  achievable  in  a  tumor.  Alternatively,  a 
femtosecond  pulsed  wave  laser  enables  significant  temperature 
increases  on  the  nanoparticle  surface  without  heating  the 
solution  (Figure  3C).  By  irradiating  our  nanoparticle  solution 
with  an  ultrafast  laser,  we  were  able  to  demonstrate  release  at 
43.9  +  1.1  °C,  ~16  °C  below  the  dehybridization  temperature. 

We  followed  up  these  observations  by  embarking  on  an  in  depth 
investigation  of  the  specific  mechanisms  of  DNA  release  from 
gold- silica  nanoshell  complexes  using  the  two  lasers:  continuous 
wave  (CW)  and  femtosecond  pulsed.  As  detailed  in  our  recently 
published  paper  entitled:  “Understanding  Resonant  Light- 
Triggered  DNA  Release  from  Plasmonic  Nanoparticles” 

(Appendix)  [4],  CW  laser  illumination  enables  photothermal  by 
dehybidizing  the  doubled-stranded  DNA  and,  thereby,  releasing 
single-stranded  DNA,  while  pulsed-laser  excitation  results  in 
double  stranded  DNA  release  by  cleavage  of  the  Au-S  bond  at  the  nanoshell  surface,  with  negligible  local 
heating.  [4] 

Subsequently,  we  examined  the  intracellular  release  of  DOX  by  both  CW  and  pulsed  lasers.  Also,  armed  with 
our  better  understanding  of  DNA  release  as  a  function  of  the  laser,  we  revisited  lapatinib  release  and  cytotoxicity. 
DOX  and  HER2-targeted  lapatinib,  were  delivered  to  MDA-MB-231  and  SKBR3  (over  expresses  HER2)  breast 
cancer  cells  and  compared  to  release  in  non-cancerous  RAW  264.7  macrophage  cells.  Continuous  wave  laser- 
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Figure  3.  Comparison  of  the 
dehybridization  temperature  achieved  by 
various  DNA  release  methods;  (A)  thermal 
release,  (B)  continuous  wave  light-triggered 
release,  and  (C)  ultrafast  pulsed  wave  light- 
triggered  release. 


Temperature  (°C) 


Temperature  (°C) 


induced  release  of  DOX  from  a  DNA  scaffold 
showed  increased  cell  death,  but  also  resulted  in  non¬ 
specific  cell  death  from  the  laser  itself  (Figure  4A). 
For  lapatinib  release,  the  pulsed  laser  treatment  using 
a  human  serum  albumin  protein  scaffold  to  envelope 
the  drug  showed  increased  cell  death  in  cancerous 
cells,  while  non-cancerous  macrophages  were 
unaffected.  These  results  show  that  a  pulsed  laser  can 
induce  drug  release  in  vitro  using  low  levels  of  near- 
IR  light  without  causing  non-specific  cell  death, 
unlike  continuous  wave  light- triggered  release.  Using 
the  hydrophobic  pocket  of  a  protein  scaffold  to  load 
and  release  hydrophobic  drug  molecules  is  a 
favorable  strategy  for  controlled  drug  delivery  that 
can  be  extended  to  myriad  of  cancers  and 
chemotherapy  drugs.  We  have  now  come  full  circle 
and  have  devised  a  method  in  which  to  deliver 
lapatinib  using  gold-silica  nanoshells. 


Figure  4.  Nanocomplex  formation.  (A)  Formation  of 
NS@DNA@Drug  nanocomplex  by  functionalization  of 
NS  with  thiolated  dsDNA  through  Au-thiol  bonds, 
followed  by  attachment  of  drug  molecules.  Release 
triggered  by  CW  irradiation  results  in  dsDNA 
dehybridization  and  release  of  the  drug,  while  pulsed 
irradiation  results  in  Au-S  bond  breakage  and  release  of 
the  dsDNA-drug  complex.  (B)  Formation  of 
NS@HSA@Drug  by  coating  NS  with  HSA,  and 
attachment  of  drug  molecules.  Drug  release  from  the 
protein  is  triggered  by  CW  or  pulsed  irradiation. 


Task  2 

a.  Current  objective 

To  optically  release  the  drug  attached  to  the  surface  of  the  gold  nanoparticle,  it  is  necessary  to  deliver  and  validate 
the  instantaneous  fluence  (the  radiant  energy  per  unit  surface  area)  at  the  location  of  the  tumor.  Given  that,  there 
is  little  data  on  the  amount  and  distribution  of  near  infrared  light  through  the  skull  and  brain  tissue.  To  provide 
an  estimate  of  the  energy  transport  of  near  infrared  light  through  the  skull  and  within  brain  tissue,  a  number  of 
diagnostic  tools  were  developed.  First,  a  nearly  isotropic  optical  dosimetry  probe  was  fabricated  to  obtain  local 
photon  fluence  measurements  with  1 .5  mm  spatial  resolution  and  high  sensitivity.  Calibration  of  the  probe  was 
performed  in  air  and  water-based  phantoms  using  a  calorimeter  (considered  a  Gold  Standard).  Second,  a  voxel- 
based  Monte  Carlo  simulation  code  written  in  CUDA  for  GTX  Nvidia  (Dr.  Fang,  Harvard)  was  further  developed 
for  different  optical  beam  properties,  such  as  super  Gaussian  beams  with  parallel,  converging,  and  diverging  beam 
profiles,  and  with  the  ability  to  run  on  multiple  Nvidia  cards  (n=3).  Third,  an  empirical  method  to  simulate  photon 
transport  was  developed  to  provide  sufficiently  accurate  solution  at  accelerated  rates.  Monte  Carlo  simulations 
for  photon  transport  are  computational  intensive  and  therefore  are  typically  not  practical  for  clinical  use.  These 
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codes  provide  solutions  in  real-time  during  the  developmental  stage  of  the  treatment  plan.  These  tools  were  used 
to: 


2a.  Estimate  of  the  attenuation  of  near  infrared  irradiation  by  the  skull  using  the  fluence  probe;  and 

2b.  Estimate  of  the  attenuation  of  near  infrared  irradiation  by  the  skull  and  brain  tissue  phantom  using  the 
fluence  probe. 


b.  Results 


i.  Calibration  of  the  Optical  Dosimetry  Probe  (ODP)  [5,  6]:  The  objective  was  to  fabricate  and  calibrate  the 
optical  dosimetry  probe  (ODP)  to  a  calorimeter  (Nova  II;  Ophir).  The  ODP  consists  of  a  1.5  mm  diameter  sphere 
(epoxy  and  titanium  dioxide)  molded  onto  the  end  of  a  200pm  optical  fiber.  The  output  of  the  fiber  was  coupled 
to  a  photodiode  (ET-2030,  Electro  Optics  Technology  Inc),  the  output  of  which  was  recorded  by  a  digital 
oscilloscope  (Tektronix  TDS3052B).  The  digitized  signal  was  uploaded  to  a  PC,  where  the  signal  was  passed 
through  a  high-pass  butterworth  filter  to  reduce  background  noise  and  the  area  underneath  the  pulse  calculated  to 
determine  the  photon  fluence.  To  calibrate  the  ODP,  the  photon  fluence  in  air  was  measured  by  both  the  ODP  and 
calorimeter  using  a  uniform  laser  beam  formed  by  passing  the  output  of  the  laser/OPO  system  (Nd:YAG  class  4 
laser/  MagicPRISM  OPO,  Opotek  Inc.)  (tuned  to  800nm)  through  an  integrating  sphere  (Melles  Griot).  The  results 
are  displayed  in  Figure  1.  To  determine  the  fluence  in  water  (or  tissue),  the  change  in  the  index  of  refraction 
causes  the  reflective  properties  between  the  probe  media  and  water  to  change,  thus  requiring  an  additional 
correction  factor.  The  probe  was  placed  within  a  macrocuvette,  a  plexiglass  container  which  was  approximately 
5.0  x  5.0  x  5.0  cm3  (5mm  thick  walls).  The  ODP  signal  and  calorimeter  measurements  were  acquired  in  air  as 
before.  The  macrocuvette  was  then  filled  with  water,  and  the  ODP  was  positioned  along  the  central  axis  of  the 
uniform  beam.  From  this  data,  a  conversion  factor  was  applied  to  the  ODP  signal  to  convert  the  calibration  curves 
in  air  (Figure  1)  to  tissue.  This  conversion  factor  also  took  into  account  the  variation  in  the  beam  width  when 
traversing  from  plexiglass-to-air  versus  plexiglass-to-water  and  photon  attenuation  for  optical  path  lengths  within 
water,  both  of  which  were  relative  small  (~3%,  the  average  error  in  the  laser  output  power). 


Equatorial  3nfj  Azimuthal  Uniformity  (Pro be  IjQtnqprcity) 


Probe  Linearity 


~  0.4 


5 

T 


J 


I  Equatorial 

Variation 

I  Azimuthal 
Variation 


Figure  1 :  (Left)  The  output  of  the  ODP  in  mVolts  plotted  as  a  function  of  laser  fluence  measured  by  a 
calorimeter.  (Right)  The  angular  dependency  (azimuthal  and  rotation)  of  the  ODP  to  the  direction  of  laser 
light  was  calculated  as  a  function  of  the  absolute  value  of  the  residuals  and  how  it  compares  to  nylon 
probes. 
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ii.  Phantoms  Reproducing  the  Optical  Properties  of  the  Head:  The  objective  was  to  develop  a  formula  that 
could  be  used  to  create  water-based  (and  agar-based)  phantoms  with  the  desired  optical  properties  (pa  and  ps’) 
for  the  skull,  white  matter  and  gray  matter  (at  800nm,  the  wavelength  of  stimulated  drug  release).  India  ink  was 
used  as  the  absorber  (with  some  but  little  scattering),  while  intralipid  was  used  as  a  scatterer  (with  some  but  little 
absorption).  To  determine  the  absorption  coefficient,  various  concentrations  of  india  ink  in  water  were  placed 
within  the  photospectrometer,  from  which  a  calibration  curve  plotted:  15.07  x  %  (india  ink)  +  0.0049  (r2=0.99) 
(data  not  shown).  The  reduced  scattering  coefficient  was  determined  by  placing  different  concentrations  (from  1 
to  7%)  of  intralipid  in  water,  exposing  the  solution  to  a  uniform  (1cm  diameter)  laser  beam  using  an  integrating 
sphere  as  described  in  the  previous  section.  The  photon  fluence  was  measured  using  the  ODP  at  various  path 
lengths  (or  depths)  along  the  central 
axis  of  the  beam.  To  determine  the 
reduced  scattering  coefficient  in  each 
sample,  Monte  Carlo  simulations 
were  performed  for  various  reduced 
scattering  coefficients  until  a 
minimum  in  the  chisquare  was 
reached.  The  data  was  plotted  and  fit  to  a  line:  1 1.30  x  (%intralipid)  -  2.17  (r2=0.99)  (data  not  shown).  From  these 
measurements,  the  white  and  gray  matter  phantoms  were  fabricated  based  the  optical  properties  measured  by  Van 
der  Zee  et  al  [7],  and  the  skull  phantom  optical  properties  were  based  on  measurements  by  Firbank  et  al  [8],  see 
Table  I. 


Table  I:  Optical  properties  of  the  skull  bone,  white  matter,  and  gray 
matter. 


(US 

g 

ps' 

intralipid 

[cm"1] 

[cm'1] 

white 

0.05 

0.85 

82.5 

7.5 

0 

gray 

0.35 

0.965 

24.5 

2.36 

0.023 

skull 

0.24 

1.82 

0.0156 

iii.  Simulating  Photon  Transport  in  Phantoms/Tissue  [5,  6,  9]:  The  objective  was  to  develop  and  validate  a 
voxel-based  Monte  Carlo  code  of  photon  transport  in  the  head  and  brain.  A  voxel-based  3-D  Monte  Carlo  (MC) 
code  (Fang  et  al  [10])  developed  in  CUDA  for  NVidia  cards  was  modified  to  allow  the  user  to  create  laser  sources 
with  various  aperture  profiles,  such  as  uniform  and  super-Gaussian,  and  directionality,  such  as  diffusive, 
converging,  and  diverging  beams.  Supporting  code  (Matlab)  was  developed  to  allow  the  user  to  design  their  own 
objects,  which  in  this  study  was  based  on  phantom  designs  and  medical  images  (CT). 

A  critical  advantage  of  a  voxel-based  code  is  not  only  to  include  tissue  heterogeneities  but  to  integrate  medically 

acquired  images  into  the  simulations.  To  date,  the  MC  code  has  not  been  validated  with  measured  data;  therefore, 

a  number  of  experiments  were  performed.  First,  the  photon  fluence  distribution  was  measured  using  the  ODP  in 

phantoms  representing  white  matter,  gray  matter,  and  skull  bone.  The  ODP  probe  was  secured  to  a  3-D  stage 

micrometer  to  obtain  measurements  along  the  central  axis  of  the  beam  at  various  depths  in  the  phantom  (axial 
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profile)  and  perpendicular  to  the  central  axis  (lateral  profile).  These  measurements  were  taken  using  a  uniform 


and  a  super  Gaussian  laser  beam.  Next,  Monte  Carlo  simulations  were  performed  and  compared  to  measurement. 
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Figure  2.  The  fraction  of  the  input  fluence  as  measured  by  the  ODP  at  different  depths  in  white  and  gray  matter 
phantoms  for  a  uniform  laser  beam  (1.0cm  in  diameter)  (left)  and  a  super  Gaussian  beam  profile  (right).  Included 
in  these  plots  are  the  Monte  Carlo  simulated  results,  which  agree  to  within  3.0%  of  the  measured  values. 


In 


Figure  2,  the  normalized  fluence  along  the  central  axis  of  the  flat  beam  was  plotted  for  both  white  and  gray  and 
for  ODP  and  MC  data.  For  depths  up  to  2  cm,  the  agreement  is  within  3%.  Similar  experiments  were  performed 
for  the  super  Gaussian  beam,  where  for  depths  up  to  3  cm,  the  MC  simulated  results  matched  the  ODP 
measurements.  The  comparison  of  the  lateral  profiles  for  the  super  Gaussian  beam  is  shown  in  Figure  3,  and 


demonstrates  the  excellent  agreement  between  ODP  and  MC  results. 


Figure  3.  Plotted  are  the  relative  photon  fluence  as  measured  by  the  ODP  and  as  simulated  by  the  MCcode  for  a  super  Gaussian 
laser  beam  in  air  (top  left),  white  matter  ( top-right ),  gray  matter  (bottom-left),  and  skull  (bottom-right). 
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iv.  Sheep  Skull  and  Brain  Phantom:  In  combination  with  the  phantom  study  from  the  previous  section,  an 
estimate  of  the  attenuation  of  near  infrared  irradiation  by  the  skull  and  brain  tissue  phantom  was  made.  The  laser 
beam  emanating  from  the  optical  fiber  bundle  was  measured  in  air  using  the  ODP  and  used  to  determine  the  beam 
profile,  e.g.,  super  Gaussian  function,  and  laser  fluence.  This  laser  beam  was  positioned  as  close  to  normal  to  the 
sheep  skull  as  possible.  The  thickness  of  the  bone  where  the  laser  beam  was  exposed  was  measured  to  be  1.1cm 
using  calipers  and  validated  by  CT  scans.  To  measure  the  laser  beam  profile  at  the  inner  surface  of  the  skull,  the 
dosimetry  probe  was  axially  positioned  within  4  mm  of  this  inner  surface,  and  using  a  3-D  micrometer,  the  ODP 
was  translated  perpendicular  to  the  laser  beam  to  measure  the  fluence  profile  and  along  the  axis  of  the  beam 
profile  to  determine  the  peak  beam  intensity  as  a  function  of  the  distance  from  the  inner  surface  (see  Figure  4). 
After  applying  the  calibration  factors  (Section:  Calibration  of  ODP),  the  percent  of  the  laser  fluence  transmitted 
through  the  skull  was  measured  to  4.41%  (SA2.A).  MC  simulated  data  were  consistent  with  these  ODP 
measurements  (see  Figure  4). 
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Figure  4.  ( Top  row)  Displayed  are  the  photos  of  the  sheep  skull  and  the  ODP  used  in  the  phantom  experiments.  The  picture  on 
the  right  has  the  laser  beam  exposing  the  frontal  bone  (left  to  right)  and  the  ODP  embedded  in  the  brain  tissue  phantom.  ( Bottom 
row)  Plotted  are  the  lateral  profile  near  the  inner  surface  of  the  skull  and  the  axial  profile  along  the  central  axis  of  the  laser  beam 
as  measured  by  the  ODP  probe  and  simulated  using  the  MC  code. 


With  this  same  setup,  the  skull  cavity  was  filled  with  a  4%  intralipid  solution.  The  ODP  was  used  to 
measure  the  attenuation  of  the  beam  (Figure  4.  Top  right,  bottom  right):  the  transmitted  fluence  was  measured  to 
be  7.25%  per  cm,  which  is  consistent  with  Figure  2.  If  a  7.5%  of  intralipid  (white  matter)  was  used,  the  measured 
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transmitted  fluence  was  3.1%  per  cm,  and  gray  matter  was  measured  to  be  0.31%  per  cm  (see  Figure  2).  To  relate 
these  measurements  to  a  human  patient,  average  values  for  the  thickness  of  the  cortical  gray  matter  and  skull  bone 
is  determined.  Cortical  gray  matter  thicknesses  can  range  from  2.0  to  5.0  mm,  averaging  around  4.7-4.85  mm  [11, 
12].  Skull  bone  thickness  can  range  from  5.0  to  11.0  mm  in  an  adult  human  head,  averaging  around  7.8  mm. 
Based  on  the  above  measurements,  the  average  fluence  transmitted  through  8.0  mm  of  the  skull  bone  would 
be  10.3%  and  through  4.8mm  of  gray  matter  would  be  9.1%,  or  0.93%  on  average  through  both  tissues. 
Therefore,  a  tumor  4.5cm  deep  within  the  brain  would  receive  approximately  8.5xl0'9  of  the  laser  fluence 
(SA2.B),  Depending  on  the  location  of  the  tumor,  the  positioning  of  the  optical  sources  could  take  into  account 
the  heterogeneity  of  the  skull  and  gray  matter.  Optimally,  a  region  of  the  skull  and  gray  matter  that  that  is  5mm 
(24.3%  transmitted)  and  2.5mm  (20.0%  transmitted)  thick  would  increase  the  fluence  to  1.9xl0"8,  an  improvement 
of  a  factor  of  4.4.  This  suggests  that  a  multiple  number  of  sources  optimally  positioned  for  treatment  based  on 
the  patient’s  anatomy  may  be  beneficial. 

v.  Empirical  Model  of  Photon  Transport  [6,  9]:  An  empirical  method  was  developed  to  simulate  photon 

transport  in  tissue  based  on  the  acquired  medical  images  of  an  individual  patient  (CT  and  MRI)  with  sufficient 

accuracy  and  at  accelerated  rates  for  clinical  use.  This  approach  is  based  on  the  assumption  that  the  photon  fluence 

in  a  voxel  within  a  particular  layer  (or  slice  of 

an  image)  is  the  weighted  sum  of  the  fluence 

from  its  neighboring  voxels  in  the  previous 

layer.  These  weights  determine  the  fraction  of 

photons  within  a  photon  packet  that  is 

transported  from  one  voxel  to  its  neighbors 

relative  to  its  average  (or  weighted)  direction 

of  propagation.  It  depends  not  only  on  a 

voxel-based  geometry  of  the  image  but  the 

voxel  size  and  optical  properties  of  the 

medium.  It  has  been  shown  to  faithfully 

reproduce  the  local  shadowing  effects  from 

high  absorbers  (e.g.,  blood  vessels)  and  high 

contrast  media  (gold  nanoparticles),  and 

general  attenuation  in  homogeneous  tissue 

types.  To  determine  the  weights,  MC  simulations  were  performed  in  known  media,  such  as  white  matter,  gray 

matter,  skull  bone,  and  an  Astrocytoma  (tumor).  Using  these  fluence  profiles,  the  weights  were  modified  until  the 

photon  fluence  distribution  from  the  empirical  model  matched  (minimum  in  Chisquare)  MC  simulated  results  for 

different  voxel  sizes.  A  comparison  between  the  photon  fluence  within  each  individual  tissue  type  and  for  a 
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Figure  5.  Displayed  is  a  diagram  depicting  the  local  photon  transport 
of  the  empirical  model  for  a  voxel  within  tissue.  Note,  this  takes  into 
account  the  general  direction  of  the  photon  packet  and  the  scattering 
properties  of  the  tissue. 


software  phantom  of  the  skull  and  brain  are  shown  in  Figure  6.  There  remains  a  relative  small  systematic  effect 
(over  estimate)  of  the  photon  fluence  which 
requires  further  investigation,  but  the 

agreement  is  very  good  to  develop  and 
optimize  delivery  systems  based  on  an 
individual  patient’s  anatomy.  When 
implemented  on  CUD  A,  the  time  it  takes  to 
simulate  a  human  head  with  a  broad  beam 
source  (e.g.,  super  Gaussian  beam)  was  on  the 
order  of  seconds  to  low  minutes  where  the  MC 
would  run  in  excess  of  24  hours. 


Figure  6.  Comparison  of  the  Empirical  algorithm  with  MC 
simulations  in  homogenous  (a, b)  and  heterogeneous  (c,d) 
phantoms  that  resemble  white  matter,  gray  matter  and 
astrocytoma.  The  plots  in  panels  (a),  (b),  and(d)  display  the 
exponentialdecreaseinfluencealongthecentralaxisofa 
Gaussian  beam. 

vi.  NIR  transmission  in  the  human  skull  as  a  function  of  the  thickness  and  density  of  bone 

(1)  Introduction 

The  skull  forms  an  optical  attenuator  (absorber  and  scatterer)  in  the  transmission  of  near-infra-red  light  to  the 
brain.  The  success  of  any  NIR  based  imaging  or  therapeutic  technique  depends  on  an  accurate  and  precise 
determination  of  the  optical  transmission  properties  of  the  skull.  The  skull  bone  consists  of  osteocytes  and 
calcium  deposits  in  complex  geometries  that  result  in  large  density  variations  (or  heterogeneities)  at  different 
anatomical  locations  within  the  skull.  Previous  studies  have  looked  at  the  optical  transmission  properties  of 
different  tissues  [13-19]  using  collimated  laser  beams  with  powers  ranging  from  10mW/cm2  to  25mW/cm2[13, 
14].  These  studies  included  NIR  transmission  through  human  cadaver  skulls. [14]  In  our  study,  we  used  an  NIR 
laser  source  (17.64mW/cm2)  with  a  complex  super-Gaussian  broad  beam  profile.  Our  objective  is  to  investigate 
if  the  relationship  between  voxel-based  HU  and  the  optical  attenuation  in  the  skull  bone  exists,  allowing  NIR- 
based  treatment  planning  of  drug  release  based  on  CT  scans.  Simulated  Monte  Carlo  data  was  compared  to 
measured  beam  intensity  and  distribution  implementing  the  Ti02-based  optical  dosimetry  probe.  This  data  was 
used  to  investigate  the  relationship  between  the  HU  from  CT  images  and  the  optical  attenuation  by  measuring  the 
transmitted  NIR  light  through  human  cadaver  skulls.  Critical  to  these  measurements  is  the  ability  to  quantify  low 
power  fluence  transmitted  through  the  skull  afforded  with  the  Ti02  optical  dosimetry  probe  compared  to  the 
calorimeter  used  in  previous  bone  transmission  studies  [15,  16]. 
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To  date,  studies  [8,  20-22]  have  measured  the  optical  properties  of  the  brain  tissues,  mapped  their  heterogeneity 
and  determined  their  inter-sample  variations.  In  these  studies,  a  comparison  of  the  optical  properties  of  skull  bone 
[8,  20-22]  is  shown  in  Table  II.  It  shows  a  wide  variation  in  the  optical  properties  reported  for  the  human  skull 
bone.  A  comprehensive  work  by  Bashkatov  et  al  [20]  mapped  these  variations  in  optical  properties  of  the  skull 
bone  at  different  wavelengths.  Note,  these  studies  assume  an  average  overall  measure  of  the  skull’s  attenuation. 
Our  aim  is  to  further  understand  this  complexity  caused  by  variations,  or  heterogeneity,  in  the  bone  density  and 
bone  thickness  to  map  out  an  accurate  fluence  map  of  the  Monte  Carlo  in  the  brain.  The  use  of  Computed 
Tomography  images  provides  a  3D  map  of  the  heterogeneous  bone  structure  and  allows  us  to  accurately  generate 
the  boundary  conditions  and  density  maps  using  the  Monte  Carlo. 


Table  II  Optical  properties  of  skull  bone  measured  at  800nm  wavelength  in  different  studies.  [8, 20-22]  The 
anisotropy  factor  (g)  is  0.9. 


Study  By: 

Absorption  Coefficient  (pa)  (cm4) 

Scattering  Coefficient  (ps)  (cm-1) 

Firbank  et  al 

0.24 

18.4 

Bashkatov  et  al 

0.11 

19.48 

Pifferi  et  al 

0.07 

12.5 

Ugryumova  et  al 

0.25 

21 

(2)  Materials  and  Methods 

(a)  Optical  Bone  Dosimetry  Setup 

10  Different  skull  bone  samples  were  used  for  the  optical  transmission  study.  The  light  source  consisted  of  the 
optical  fiberguide  connected  to  a  laser  source.  The  NIR  output  of  the  fiberguide  was  measured  using  a  calorimeter 
to  be  19.94mW  or  (17.64  mW/cm2)  at  800nm  wavelength.  For  transmission  studies,  the  optical  fiberguide  was 
placed  against  the  outer  surface  of  the  skull  bone,  such  that  its  aperture  was  flat  against  the  bone  surface,  thus 
providing  a  normal  exposure.  The  fluence  transmitted  through  the  bone  was  measured  using  the  TiC>2  dosimetry 
probe  by  placing  it  against  the  inner  surface  of  the  skull  bone,  opposite  the  fiberguide  on  the  other  side  of  the 
skull,  which  was  ensured  by  locating  the  point  of  maximum  fluence.  The  energy  fluence  for  5  points  or  locations 
per  skull  was  acquired,  and  at  each  location  the  measurement  was  repeated  three  times.  A  total  of  10  cadaver 
skull  bones  were  measured,  which  gave  a  total  of  50  measurement  data  points.  The  setup  for  these  experiments  is 
shown  in  figures  7.1  and  7.2. 
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Figure  7.1  Setup  to  measure  the  optical  transmission  of  NIR  light  through  the  skull  bone.  The  source  used  is  an  optical 
fiberguide  connected  to  a  laser.  The  NIR  fluence  collected  by  the  probe  tip  is  carried  by  an  optical  fiber  to  a  PIN  photodiode 
circuit,  whose  electrical  output  is  read  using  an  oscilloscope  (right). 


For  the  lateral  beam  profile  measurements,  each  skull  was  placed  in  a  large  Plexiglas  cuvette,  such  that  the  point 
of  illumination  (marked  previously  on  the  skull  surface)  was  directly  in  contact  with  the  cuvette  bottom.  The 
transparent  optical  bench  helped  in  accurately  aligning  the  skull  with  the  source.  The  experimental  setup  is  shown 
in  figures  7.2  and  7.3.  The  lateral  profile  was  first  measured  in  air  and  in  water  0.5cm  from  the  inner  skull  surface. 
While  measurement  in  air  was  for  further  validation  of  optical  transmission,  the  measurement  in  water  was  for  a 
more  realistic  scenario,  where  the  refractive  index  of  water  is  similar  to  that  of  soft  tissue. 


Figure  7.2  Skull  bone  sample  (top  portion)  used  in  the  experiment.  The  dosimetry  probe  is  also  seen. 
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Figure  7.3  Optical  dosimetry  setup  for  lateral  beam  profile  measurement  in  skull  sample. 


Figure  7.4  Optical  dosimetry  setup  for  lateral  beam  profile  measurement  in  a  skull  sample.  This  shows  the  illumination  setup 
using  the  optical  fiberguide  source. 


(b)  Acquisition  and  Processing  of  CT  Images  of  the  Skull  Bone 
Computed  Tomography  (CT)  was  used  to  map  optical  boundaries  and  density  of  the  bone  on  a  Philips  Brilliance 
clinical  CT  scanner.  The  CT  images  were  acquired  at  300  mAs  and  120  kVp,  and  reconstructed  with  a  slice 
thickness  of  0.67mm,  a  pixel  size  (xy  plane)  of  0.2344mm,  and  using  the  bone  kernel.  The  points  on  each  skull 
phantom  at  which  the  laser  source  was  applied  were  labeled  using  CT  contrast  beads  and  marked  with  a  pencil 
(see  Figure  7.5). 


Figure  7.5  shows  the  use  of  CT  contrast  beads  to  mark  the  point  of  illumination  on  the  skull  samples.  The  CT 
images  were  displayed  and  processed  using  Matlab  based  software,  which  was  custom  developed  to  identify  the 
source  and  detector  positions  and  simulate  the  bone  boundaries  on  CT  images,  as  shown  in  figure  7.5.  The  point 
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of  measurement  for  optical  transmission  lies  on  the  inner  skull  boundary  on  the  line  while  lateral  profile  is 
measured  at  0.5  cm  from  the  inner  skull  surface. 


In  order  to  identify  the  source  voxel  for  illumination,  another  line  was  drawn  tangential  to  the  skull  surface, 
passing  through  the  intersection  point  of  the  bead  and  the  outer  skull  surface  as  shown  in  figure  7.5  (B).  The  cross 
product  of  the  illumination  vector  and  the  source  tangential  vector  is  the  other  vector  which  also  lies  tangential  to 
the  skull  surface  along  the  z  axis  [0  0  1].  The  source  voxels  are  then  identified  as  those  lying  on  the  3D  circle 
defined  by  the  two  tangential  vectors.  The  parametric  equation  used  of  a  3D  circular  disk  [23]  can  be  then  used 
to  represent  the  fiberguide  source  surface  as  described  by  the  two  tangential  vectors.  Thus  one  can  accurately 
represent  any  broad  source  in  the  Monte  Carlo  at  any  arbitrary  angle  as  used  in  the  experimental  setup. 


Figure  7.5  The  CT  image  of  a  skull  bone  is  processed  in  Matlab  to  extract  the  source  directional  vector  (A)  and  in-plane  source 
vector  (B)  by  drawing  lines  (blue)  on  the  image  and  extracting  the  end  points  of  the  lines  to  calculate  the  vectors.  The  third 
vector  is  the  cross  product  of  these  two  vectors  and  is  in  the  plane  of  the  source. 
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(3)  Results 


(a)  Optical  Transmission  Measurements  in  the  Skull 

The  result  for  the  optical  transmission  for  10  different  phantoms  is  shown  in  figures  7.6,  7.7  and  7.8.  Note  that 
the  phantoms  are  numbered  from  1-4  and  6-11  as  phantom  5  was  not  used  (incomplete  CT  scan).  The  optical 
transmission  measurements  show  a  wide  variation  in  the  optical  power  transmitted,  due  to  variation  in  skull 
thickness  and  (heterogeneity  in)  bone  densities.  In  order  to  map  the  Monte  Carlo  fluence  to  the  %  transmission 
measurements  of  the  probe,  we  used  CT  images  to  determine  the  optical  path-length  through  the  bone  (i.e.  skull 
thickness  in  the  direction  of  illumination).  The  skull  thicknesses  measured  using  CT  images  showed  different 
skull  thicknesses  both  within  and  across  the  skull  samples.  The  measured  skull  thicknesses  from  the  CT  images 
are  shown  in  figure  7.9.  We  also  measured  the  lateral  beam  profiles  in  the  skull  phantoms  in  air  and  water  as 
shown  in  figure  7.10. 


Figure  7.6  Percentage  fluence  transmitted  for  phantoms  1-4  measured  using  the  probe  calorimeter  setup. 
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30 1 


□  Phantom  7 


Figure  7.7  Percentage  fluence  transmitted  for  phantoms  6-9  measured  using  the  probe  calorimeter  setup. 


Figure  7.8  Percentage  fluence  transmitted  for  phantoms  10-11  measured  using  the  probe  calorimeter  setup. 
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Figure  7.9  Skull  thicknesses  measured  across  skull  samples  using  CT  image  data.  Measurement  points  are  categorized  as 
follows:  Phantoml  (1-5),  Phantom2  (6-10),  Phantom3  (11-15),  Phantom4  (16-20),  Phantom5  (21-25),  Phantom6  (26-30),  Phantom7 
(31-35),  Phantom8  (36-40),  Phantom9  (41-45),  Phantoml  0  (45-50). 


Figure  7.10  Lateral  beam  profile  (fraction  transmitted)  measured  using  the  dosimetry  probe  in  air  and  water  0.5mm  from  the 
inner  skull  surface. 


(b)  Comparison  with  the  Monte  Carlo 

Based  on  the  three  different  sets  of  optical  properties  of  the  skull  bone  shown  in  Table  II,  three  Monte  Carlo 
profiles  were  obtained  with  a  super-Gaussian  beam  distribution.  The  fluence  profiles  along  and  perpendicular 
(lateral  profile)  are  shown  in  figures  7.1 1  and  7.12  respectively.  The  Monte  Carlo  simulations  show  that  different 
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optical  properties  of  the  skull  bone  will  give  rise  to  different  transmission  characteristics,  which  could  result  in 
different  transmission  factors  and  lateral  beam  profiles.  The  lateral  beam  profiles  due  to  different  optical 
properties  also  show  the  variations  between  the  four  models. 
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Figure  7.11  MC  simulated  fluence 
profiles  along  beam  axis  for  a 
homogeneous  skull  medium  generated 
using  Monte  Carlo  simulations  for  three 
sets  optical  properties  of  the  skull  bone 
(Ugrynova  etal ,  Pifferi  etal,  Firbank  etal , 
and  Bashkatov  et  al).  A  broad  beam 
super-Gaussian  NIR  source  (800nm 
wavelength)  was  simulated. 


Figure  7.12  MC 
simulated  lateral 
fluence  profiles  at 
different  depths 
(z=052,4,6mm)  for  a 
homogeneous  skull 
medium  generated 
using  Monte  Carlo 
simulations  for  three 
sets  optical  properties 
of  the  skull  bone 
(Ugrynova  etal,  Pifferi 
etal,  Firbank  etal,  and 
Bashkatov  etal).  A 
broad  beam  super- 
Gaussian  NIR  source 
(800nm  wavelength) 
was  simulated. 
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The  %  power  transmitted  using  the  probe-calorimeter  setup  was  then  compared  to  the  %  power  transmitted  using 
the  Monte  Carlo  simulations.  The  simulations  in  the  homogeneous  skull  phantom  were  used  in  the  form  of  a 
lookup  table  to  derive  the  MC  generated  fluence  using  interpolation.  The  skull  thicknesses  measured  using  CT 
images  were  used  to  look  up  and  interpolate  the  corresponding  normalized  fluences  in  the  MC  table.  The 
percentage  power  transmitted  is  the  normalized  fluence  on  the  beam  axis  at  that  specific  depth  multiplied  by  100. 
Figure  7.13  shows  the  comparison  between  the  transmitted  values  measured  by  the  probe  and  those  simulated  by 
the  Monte  Carlo  code. 
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Figure  7.13  Comparison  of  % 
transmitted  fluence  between 
the  calorimeter-probe 
measurements  and  the  Monte 
Carlo.  The  error  bars  in  the 
probe  measurements 
represent  the  maximum  and 
minimum  intensities  in  the  3 
averages,  while  the  error  bars 
in  the  Monte  Carlo  are  the 
minimum  and  maximum 
fluences  between  the  3 
optical  models  (see  Table  II). 


Figure  7.13  shows  that  the  Monte  Carlo  simulations  compare  closely  to  some  of  the  probe  measurement  (example 
points  1-5  in  phantom  1)  while  underestimating  the  transmission  %  for  other  points.  Thus  Monte  Carlo 
overestimates  the  attenuation  of  fluence  in  the  skull,  while  probe  readings  show  a  higher  %  power  transmitted 
through  the  skull.  The  reason  behind  this  difference  is  due  to  the  differences  in  the  bone  density  in  the  skull 
samples  which  were  not  included  in  the  MC  simulations.  The  MC  simulations  were  done  in  a  homogeneous  brain 
phantom  and  assumed  uniform  density. 


In  order  to  demonstrate  the  feasibility  of  using  CT  image-guided  skull  boundaries  for  Monte  Carlo  simulation, 
we  used  the  CT  images  to  map  the  heterogeneity  in  distribution  of  bone  voxels  along  with  the  geometry  of 
boundaries.  Figure  7.14  shows  the  comparison  between  the  lateral  fluence  profiles  (transmitted  through  the  skull 
immersed  in  water)  measured  using  the  probe  versus  the  Monte  Carlo.  The  boundary  conditions  derived  from 
the  CT  images  is  shown  in  figure  7.15.  The  fluence  profiles  were  normalized  to  the  maximum  value  measured  at 
a  depth  of  5mm  from  the  inner  surface  of  the  skull  bone.  The  results  show  a  close  match  between  the  beam  spread 
simulated  by  the  Monte  Carlo  and  measured  by  the  probe.  This  demonstrates  the  accuracy  with  which  the 
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scattering  properties  of  the  bone  were  modelled,  along  with  the  CT  image  enhanced  boundary  conditions.  The  CT 
images  of  the  skull  sample  also  show  the  changing  Hounsfield  CT  numbers  across  the  bone  which  is  the  result  of 
varying  density  of  the  bone  structure.  Thus  compared  to  previous  studies,  we  have  used  the  CT  images  to  account 
for  the  different  optical  path  lengths  along  with  boundary  conditions.  We  have  however,  not  accounted  for  the 
varying  densities  of  the  bone,  which  results  in  differences  between  the  probe  and  the  Monte  Carlo  fluences. 
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Figure  7.14  Comparison  of  lateral  profiles  measured  using  dosimetry  probe  with  the  Monte  Carlo  in  water. 
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NIR 

source 

Figure  7.15  Boundary  conditions  simulated  using  the  CT  image-guided  Monte  Carlo. 


(4)  Discussion 

The  optical  transmission  study  demonstrates  the  difficulties  encountered  in  optical  bone  transmission  studies.  The 
Monte  Carlo  can  accurately  simulate  photon  attenuation  in  the  skull  only  when  all  the  optical  properties  and 
boundary  conditions  are  simulated  correctly.  While  the  skull  bone  has  been  measured  and  modeled  as  a  relatively 
homogeneous  medium  in  optical  simulation  studies,  the  use  of  CT  images  shows  that  it  is  far  from  homogeneous 
with  wide  variations  in  bone  densities  at  different  anatomical  regions  and  across  subjects.  This  presents  a 
complicated  scenario  to  model  the  optical  propagation  of  photons  through  the  skull. 

We  demonstrated  the  accuracy  of  boundary  conditions  determined  using  CT  images  by  comparing  the  lateral 
profiles  between  the  Monte  Carlo  and  the  dosimetry  probe.  The  close  match  between  the  profiles  shows  the 
accuracy  of  our  simulations.  Our  study  thus  also  reduced  the  error  in  measuring  skull  thicknesses  using  calipers, 
which  often  goes  unnoticed  in  complex  tissue  structure.  For  example:  an  error  in  correctly  measuring  the  skull 
thickness  by  a  millimeter  can  result  in  errors  of  up  to  40%  in  magnitude  (seen  from  our  MC  data  in  figures  7.1 1, 
7.12),  depending  on  the  boundary  and  illumination  conditions.  Thus  any  imaging  or  therapeutic  technique  needs 
to  reduce  the  errors  due  to  mismatch  in  the  boundary  conditions  by  using  an  anatomical  imaging  technique  (such 
as  CT)  which  can  map  different  tissue  boundaries. 
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To  improve  the  match  between  the  Monte  Carlo  and  the  measured  fluence  in  skull  models,  we  hypothesize  that 
using  bone  density  measurement  technique  such  as  quantitative  CT  can  significantly  improve  the  localization  of 
high  scattering  bone  voxels  in  the  skull  bone.  The  mapping  of  Hounsfield  values  of  CT  number  can  be  directly 
linked  to  the  bone  density  values  by  imaging  bone  density  phantoms  in  the  CT  scanner.  This  would  provide 
volumetric  images  of  bone  densities  which  can  be  converted  (scaled)  to  the  corresponding  scattering  absorption 
coefficients.  The  optical  coefficient  scaling  will  also  need  to  be  calibrated  by  measuring  the  optical  properties  of 
bone  samples  of  known  density. 

Currently,  we  are  developing  iterative  techniques  to  identify  a  relationship  between  HU  and  reduced  scattering 
coefficient  in  bone  on  a  voxel-wise  basis.  This  is  anticipated  to  improve  the  accuracy  and  precision  in  the 
transmission  and  lateral  beam  profile  through  the  skull  bone.  A  second  line  of  investigation  is  the  change  in 
wavelength.  This  has  the  potential  to  reduce  the  effects  of  scattering  but  at  the  cost  of  water  absorption.  This 
tradeoff  may  help  focus  the  beam  for  smaller  depths  beyond  the  skull  within  the  brain  tissue.  However,  the 
nanoparticle  geometry  will  significantly  change  in  order  to  provide  stimulated  drug  release,  and  whether  this  is 
possible  within  macrophages  is  yet  unexplored. 

Tasks  1  d  &e,  2c,  and  Task  3:  These  goals  were  unable  to  be  met  because  of  the  issues  with  drug  solubility 
(laptinib),  light-triggered  release,  and  probable  leakage  of  doxorubicin  from  the  NS@DNA@DOX  nanocarriers. 
They  all  required  the  delivery  of  NS@DNA  with  stably  loaded  drug  (DOX),  which  could  be  released  by  laser 
irradiation.  The  work  around  with  the  protein  scaffold  and  femtosecond  laser  was  not  realized  until  late  in  the 
grant  cycle. 

4.  KEY  RESEARCH  ACCOMPLISHMENTS:  Bulleted  list  of  key  research  accomplishments  emanating  from 
this  research.  Project  milestones,  such  as  simply  completing  proposed  experiments,  are  not  acceptable  as  key 
research  accomplishments.  Key  research  accomplishments  are  those  that  have  contributed  to  the  major  goals 
and  objectives  and  that  have  potential  impact  on  the  research  field.  If  there  is  nothing  to  report,  simply  state 
"Nothing  to  report. " 

•  Nanocarriers 

Earlier  experiments  by  the  Rice  Nanophotonics  Group  suggest  that  DNA  dehybidization  and  cargo  release 

occurred  by  a  non-thermal  release  mechanism  involving  hot-electrons. [3]  The  collective  oscillation  of  surface 

conduction  electrons  in  metallic  nanoparticles  induced  by  light  at  a  resonant  frequency  is  known  as  a  localized 

surface  plasmon.  Upon  excitation,  surface  plasmons  can  decay  radiatively  by  re-emitting  photons  of  light  or 

non-radiatively  by  the  generation  of  local  heat  and/or  energetic  electrons.  These  energetic  electrons  are  known 

as  “hot-electrons”  and  the  Rice  group  had  hypothesized  that  the  two  DNA  strands  on  the  NS  surface,  which 

are  negatively  charged,  were  dehybridized  by  electrostatic  interaction  with  the  hot-electrons.  It  is  now  clear 
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from  the  experiments  conducted  as  a  part  of  this  research  that  at  particle  concentrations  achievable  within  a 
tumor,  dehybridization  of  the  DNA  double  helix  is  a  thermal  phenomenon.  This  resolves  a  fundamental 
question  as  to  how  the  NS  @  DNA  nanocarriers  function  to  release  cargo.  It  has  also  led  to  the  testing  of  the 
femtosecond  pulsed  laser  as  a  work-around,  which  releases  cargo  by  breaking  the  the  Au-S  bond  at  the 
nanoshell  surface  but  results  insufficient  bulk  heating  to  cause  protein  denaturation  and  subsequent  tissue 
damage. 

Many  targeted  therapies  are  small  hydrophobic  molecules;  of  note-this  includes  almost  the  entire  family  of 
tyrosine  kinase  inhibitors.  We  have  developed  a  new  strategy  which  uses  the  hydrophobic  pocket  of  a  protein 
scaffold  on  the  nanoshell  to  load  and  from  which  to  release  hydrophobic  drug  molecules  for  controlled  drug 
delivery  that  can  be  extended  to  myriad  of  cancers  and  chemotherapy  drugs. 

•  Delivery  of  infrared  irradiation  to  brain 

As  pointed  out  earlier,  prior  to  our  work  there  was  little  data  on  the  amount  and  distribution  of  near  infrared 
light  able  to  be  delivered  through  the  skull  and  brain  tissue.  We  have  been  able  to  model  this  using  both  the 
sheep  skull  and  phantom  as  well  as  the  Monte  Carlo  simulations.  This  has  led  us  to  proposing  that  multiple 
simultaneous  sources  of  IR  irradiation  with  or  without  intracranial  intervention  methods  is  likely  to  be 
required  to  deliver  a  sufficient  fluence  of  irradiation  to  effect  drug  release. 

We  have  demonstrated  that  using  CT  to  map  skull  geometry  and  determine  bone  density  enables  a  better 
understanding  of  the  attenuation/scattering  of  NIR  transmission  through  the  skull.  Based  on  our  preliminary 
data  and  NIR  transmission  results,  MC  simulated  results  utilizing  the  average  optical  properties  of  the  skull 
(accounting  for  thickness)  obtained  from  literature  underestimated  the  measured  light  distribution,  primarily 
due  to  the  heterogeneity  in  bone.  Further  work  will  be  performed  to  provide  a  voxel-based  method  to  account 
for  this  heterogeneity  and  improve  our  ability  to  predict  the  photon  fluence  rate  distribution  within  the  brain. 

5.  CONCLUSION:  Summarize  the  importance  and/or  implications  with  respect  to  medical  and  /or  military 
significance  of  the  completed  research  including  distinctive  contributions,  innovations,  or  changes  in  practice  or 
behavior  that  has  come  about  as  a  result  of  the  project.  A  brief  description  of  future  plans  to  accomplish  the  goals 
and  objectives  shall  also  be  included. 

The  delivery  of  effective  therapeutics  across  the  blood-brain  barrier  remains  as  critical  a  clinical  challenge  as  it 
was  when  we  first  proposed  our  study.  We  remain  convinced  that  the  active  delivery  by  macrophages  of 
nanobound  therapeutics  to  intracerebral  malignancies  is  a  promising  strategy.  Our  studies  have  provided 
fundamental  information  on  a  number  of  aspects  of  the  delivery  of  therapeutics  by  NS@DNA.  1.  It  is  of  key 
importance  that  the  therapeutic  is  soluble  and  stable  in  aqueous  solution  if  the  DNA  duplex  is  to  be  used  as  a 
carrier;  2.  Dehybridization  of  DNA  duplexes  on  the  nanoshell  surface  is  a  thermal  phenomenon  when  the  NIR 
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source  is  a  continuous  wave  laser.  A  femtosecond  pulsed  laser  cleaves  the  Au-S  at  the  nanoshell  surface  releasing 
the  DNA.  3.  Although  doxorubicin  is  a  DNA  intercalating  agent,  it  may  distort  the  double  helix  such  that  it  is 
not  a  secure  cargo.  4.  Other  strategies  are  required  to  complex  to  the  nanoshells  small  molecules  that  are  poorly 
soluble  in  water.  We  have  demonstrated  that  the  combination  of  a  protein  scaffold  in  which  to  envelope  the  drug 
on  the  nanoshell  and  a  femtosecond  pulsed  laser  is  a  promising  strategy  for  controlled  drug  delivery.  5.  We  have 
estimated,  to  the  best  of  our  knowledge  for  the  first  time,  the  amount  and  distribution  of  near  infrared  light  able 
to  be  delivered  through  the  skull  and  brain  tissue.  We  have  also  begun  to  understand  how  inhomogeneity  of  the 
skull  bone  with  regard  to  thickness  and  density  scattering  NIR  and  how  this  may  be  modeled  to  optimize 
transmission.  While  we  have  demonstrated  substantial  attenuation  of  the  IR  irradiation  on  its  pathway  through 
the  skull  and  brain  phantom,  our  results  have  prompted  us  to  consider  utilizing  multiple  NIR  sources. 

Although  this  marks  the  end  of  our  DoD  funding,  we  hope  to  continue  to  pursue  this  novel  delivery  method. 
Should  we  secure  future  funding,  we  will  determine  if  and  to  what  extent  DOX  is  leaking  from  the  DNA  duplex; 
presuming  that  DNA  duplex  distortion  is  responsible  for  the  leakage,  we  will  search  for  an  alternate  drug  that  can 
be  maintained  within  the  dsDNA@NS  complex  until  release  is  triggered;  and  we  will  follow  up  on  our  observation 
that  hydrophobic  drugs  can  be  delivered  in  and  released  from  protein  pockets.  With  regard  to  delivering  sufficient 
IR  irradiation  to  the  brain,  the  following  are  strategies  and  fields  of  investigation  that  we  hope  to  pursue: 

(1)  Place  multiple  sources  to  avoid  thicker  regions  of  gray  matter  and  overall  skull  density. 

(2)  Increase  the  number  of  sources  around  the  skull  in  order  to  get  the  fluence  to  plateau  (instead  of 
exponentially  decay)  at  larger  depths  in  the  brain. 

(3)  Examine  the  possibility  of  directional  scattering  in  neuronal  fibers  as  a  means  to  increase  NIR  depth 
penetration. 

(4)  Given  the  diffusive  nature  of  light,  the  beam  surface  area  can  be  increased  to  reduce  surface  heating  near 
the  surface  of  the  brain  without  too  much  loss  in  the  fluence  at  depth. 

(5)  Examine  different  wavelengths.  For  example,  look  at  some  IR  windows.  The  advantage  is  less  scattering 
at  the  cost  of  water  absorption.  The  goal  is  to  determine  how  much  of  an  improvement  on  fluence  localization 
can  be  achieved. 

(6)  Continue  to  study  the  dispersive  effects  for  pulse  laser  pulses. 
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ABSTRACT 

Purpose:  The  purpose  of  this  study  is  to  develop  an  alternate  empirical  approach  to  estimate  near-infra-red  (NIR)  photon 
propagation  and  quantify  optically  induced  drug  release  in  brain  metastasis,  without  relying  on  computationally 
expensive  Monte  Carlo  techniques  (gold  standard).  Targeted  drug  delivery  with  optically  induced  drug  release  is  a  non- 
invasive  means  to  treat  cancers  and  metastasis.  This  study  is  part  of  a  larger  project  to  treat  brain  metastasis  by 
delivering  lapatinib-drug-nanocomplexes  and  activating  NIR-induced  drug  release.  The  empirical  model  was  developed 
using  a  weighted  approach  to  estimate  photon  scattering  in  tissues  and  calibrated  using  a  GPU  based  3D  Monte  Carlo. 
The  empirical  model  was  developed  and  tested  against  Monte  Carlo  in  optical  brain  phantoms  for  pencil  beams  (width 
1mm)  and  broad  beams  (width  10mm).  The  empirical  algorithm  was  tested  against  the  Monte  Carlo  for  different  albedos 
along  with  diffusion  equation  and  in  simulated  brain  phantoms  resembling  white -matter  (p*’ =8. 25  mm"1,  pa=0. 005mm'1) 
and  gray-matter  (ps’=2.45mm1,  pa=0. 035mm'1)  at  wavelength  800nm.  The  goodness  of  fit  between  the  two  models  was 
determined  using  coefficient  of  determination  (R- squared  analysis).  Preliminary  results  show  the  Empirical  algorithm 
matches  Monte  Carlo  simulated  fluence  over  a  wide  range  of  albedo  (0.7  to  0.99),  while  the  diffusion  equation  fails  for 
lower  albedo.  The  photon  fluence  generated  by  empirical  code  matched  the  Monte  Carlo  in  homogeneous  phantoms 
(R2=0.99).  While  GPU  based  Monte  Carlo  achieved  300X  acceleration  compared  to  earlier  CPU  based  models,  the 
empirical  code  is  700X  faster  than  the  Monte  Carlo  for  a  typical  super-Gaussian  laser  beam. 

Keywords:  Empirical  algorithm  for  light  propagation,  NIR  photon  propagation,  Monte  Carlo  simulation,  Optical 
therapy,  Optically  induced  drug  release,  Tissue  Optics,  Light  Scattering 


1.  INTRODUCTION 

Optical  therapy  using  NIR  light  delivery  requires  the  estimation  of  photon  fluence  in  tissues  and  heterogeneous  media. 
This  study  is  a  part  of  a  larger  project  which  seeks  to  selectively  deliver  macrophages  laden  with  lapatinib -gold- 
nanocomplex  drug  molecules  to  brain  metastasis  with  over-expressed  HER2+  receptors  and  activate  targeted  drug 
release  using  NIR  photons.1,2  The  effectiveness  of  achieving  a  targeted  drug  delivery  and  optically  activated  drug  release 
with  a  Monte  Carlo  based  optical  therapeutic  protocol  has  been  demonstrated  in  previous  studies.1'5  The  design  of  an 
optical  therapeutic  protocol  requires  the  evaluation  of  various  illumination  conditions  such  as  the  angles  of  incidence  and 
source  positions  for  optimum  drug  release  in  the  metastatic  lesions  of  the  brain.5  While  the  Monte  Carlo  is  the  gold 
standard  for  photon  estimation  in  3D  media,  it  is  computationally  expensive  to  be  used  over  a  large  number  of  iterations. 
Hence  there  is  a  need  for  a  computationally  efficient  evaluation  algorithm  which  can  reduce  the  possible  illumination 
conditions  to  a  few  iterations,  while  minimizing  error  in  fluence  estimation.  This  algorithm  should  also  be  able  to 
interface  seamlessly  with  any  imaging  modality  system  to  effectively  localize  fluence  distributions  with  respect  to 
anatomical  and  physiological  information  provided  by  imaging  modalities  such  as  CT,  Ultrasound,  PET,  MRI,  etc.  In  the 
current  work  we  have  developed  a  fast  voxel  based  empirical  approach  which  can  estimate  photon  fluence  in  the  brain 
and  evaluate  different  illuminations  before  using  Monte  Carlo  as  a  final  tool  for  validation.  Being  a  3D  voxel  based 
method;  this  algorithm  can  be  easily  integrated  with  other  imaging  modalities  to  evaluate  fluence  distributions  in  the 
brain  and  other  tissues,  by  co -registering  3D  fluence  images  with  images  from  CT,  MR,  Ultrasound,  etc. 

Light  propagation  studies  in  3D  media  rely  on  the  use  of  Monte  Carlo  simulations  as  the  gold  standard,  to  estimate  the 
photon  fluence  distribution  in  three  dimensional  units  of  volume  known  as  voxels.  The  Monte  Carlo  model  uses 
probability  distribution  to  generate  a  random  walk  of  photons  in  three  dimensional  media  and  relies  on  generating  a  large 
number  of  events  to  produce  statistically  significant  fluence  distributions  in  optical  media  with  specific  absorption  and 
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scattering  properties.6'10  The  photon  fiuence  deposited  (and  traversing)  within  each  voxel  depends  on  the  specified 
optical  properties  assigned  by  the  user,  such  as  the  absorption  coefficient  (pa),  scattering  coefficient  (ps),  anisotropy 
factor  (g)  and  refractive  index  (n).  The  traditional  Monte  Carlo  software  is  computationally  expensive  as  it  relies  on 
generating  millions  of  photons  and  tracking  each  photon  as  it  propagates  through  the  medium.  Recent  research  in  GPU 
based  Monte  Carlo  by  the  Fang  group  has  resulted  in  speeds  up  to  300X  faster  than  the  traditional  CPU  based  versions.6 
However  for  large  optical  datasets,  such  as  those  obtained  and  translated  from  CT  or  MR  images,  and  under  wide  beam 
illumination  conditions,  the  Monte  Carlo  suffers  from  slow  computation  speeds  due  to  large  photon  propagation  times  in 
highly  scattering  tissues.  This  results  in  a  drawback  and  slowdown  in  the  process  of  planning  targeted  light  delivery  to 
specific  regions  within  the  human  body.  In  this  study  we  have  designed  an  alternate  empirical  approach  to  light 
propagation  which  can  estimate  the  photon  fiuence  in  tissues  in  complex  three  dimensional  media  to  provide  a  fast 
preliminary  approach  to  light  estimation  while  planning  targeted  NIR  light  delivery.  The  purpose  of  this  study  is  to 
devise  a  fast  means  to  estimate  NIR  photon  fiuence  under  different  illumination  conditions  in  an  optical  treatment  plan 
before  using  Monte  Carlo  for  final  verification  and  corrections. 

The  empirical  approach  estimates  the  fiuence  distribution  in  a  voxelated  medium  by  using  a  weighted  approach  to 
fiuence  propagation  from  one  layer  to  the  next.  The  assumption  used  is  that  the  photon  fiuence  in  a  voxel  is  a  weighted 
sum  of  the  fluences  in  the  neighboring  voxels  of  the  preceding  layer.  The  weights  are  arbitrarily  assigned  to  each  voxel 
based  on  its  optical  properties,  mainly  the  scattering  coefficient.  The  accuracy  of  the  empirical  approach  depends  on  the 
optimization  method  used  to  estimate  the  weights  for  each  set  of  optical  properties.  The  empirical  approach  simulates 
photon  scatter  by  dividing  the  fiuence  into  scattering  vectors  in  the  forward  and  diagonal  directions,  while  photon 
absorption  is  simulated  using  the  Beer  Lambert’s  rule.11  An  optimization  routine  is  used  to  estimate  the  weights  assigned 
to  each  voxel  by  comparing  with  a  pencil  beam  simulated  by  the  Monte  Carlo.  In  this  study  we  estimated  the  weights  for 
typical  brain  tissues  such  as  Gray  Matter,  White  Matter,  Skull  bone  and  Astrocytoma  with  known  optical  properties.12'14 
These  weights  were  then  tested  in  homogeneous  and  heterogeneous  brain  phantoms  under  pencil  beam  and  broad  beam 
illumination.  The  Monte  Carlo  has  been  used  as  the  gold  standard  for  evaluating  the  performance  of  this  algorithm.  The 
3D  GPU  based  Monte  Carlo,  developed  by  the  Fang  group,  has  an  acceleration  of  300X  over  the  traditional  CPU  based 
version6,  while  the  Empirical  approach  was  seen  to  be  700X  faster  than  the  GPU  based  Monte  Carlo  for  a  super- 
Gaussian  beam. 


2.  MATERIALS  AND  METHODS 


2.1  Algorithm  Design 

The  empirical  approach  approximates  the  photon  scatter  from  one  voxel  to  the  next  by  assigning  weights  to  the 
scattering  directions.  The  main  goal  of  this  algorithm  is  to  approximate  the  total  scattering  events  within  a  voxel  into  14 
component  vectors.  Since  tissues  are  mostly  forward  propagating  with  anisotropies  approximately  equal  to  or  greater 
than  0.9,  it  can  be  assumed  that  most  photon  scatter  directions  can  be  approximated  into  forward  and  diagonal  vectors. 
Hence,  each  scattering  event  is  divided  into  4  scattering  components:  one  forward  component  with  the  wtl  and  4 
diagonal  components  with  weight  wt2  as  shown  in  the  figures.  The  fiuence  enters  a  voxel  from  14  neighboring  voxels. 
Each  incoming  fiuence  component  is  then  scattered  into  one  forward  component  and  4  diagonal  components.  The  total 
fiuence  in  each  individual  direction  is  accumulated  and  absorbed  within  a  voxel  based  on  the  absorption  coefficient.  The 
remaining  fiuence  then  exits  the  voxel  into  the  14  neighboring  voxels.  It  is  important  to  note  that  the  fiuence  entering  a 
voxel  is  equal  to  the  sum  of  fiuence  absorbed  within  the  voxel  and  the  fiuence  exiting  the  voxel.  This  is  achieved  by 
ensuring  that  the  sum  of  the  weights  used  is  equal  to  one  [i.e.  (wtl  +  4  x  wt2)  =1]. 

The  steps  of  the  algorithm  can  be  summarized  as  follows: 

1)  Assign  the  incident  photon  fiuence  (e.g.  Gaussian  beam)  to  the  first  layer  of  voxels.  Here  we  have  assumed  that  the 
photons  enter  the  volume  (of  dimensions  dimlxdim2xdim3)  at  z=l  (first  layer). 

2)  Begin  the  propagation  from  layer  z=l  to  z=dim3. 
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3)  For  every  voxel  in  a  layer  z=k: 

a)  Calculate  the  14  directional  components  from  the  incident  photon  fluence.  For  the  first  layer  the  fluence  is 
assumed  to  be  in  the  forward  direction.  The  fluence  from  neighboring  voxels  is  accumulated  in  each  direction 
as  follows: 

Forward  Component  =  If(x,y,z)  =  £  (forward  components  x  wtl)  +  X(forward  diagonal  components  x  wt2) 

Forward  Diagonal  Component  =  Id(x,y,z)  =  X  (forward  components  x  wt2)  +  X(  forward  diagonal  components 
x  wtl)  +  +  X(side  components  x  wt2) 

Side  Component  =  Is(x,y,z)  =  £(  forward  diagonal  components  x  wt2)  +  £(side  components  x  wtl) 

+  X(  backward  diagonal  components  x  wt2) 

Backward  Diagonal  Component  =  Idb(x,y,z)  =  JXside  components  x  wt2)  +  ^(backward  diagonal  components 
xwtl)  +  H  backward  component  x  wt2) 

Backward  Component  =  Ib(x,y,z)  =  X  (backward  component  x  wtl)  +  ^(backward  diagonal  components  x  wt2) 
Total  Fluence  =  I(x,y,z)  =  If(x,y,z)  +  XId(x,y,z)  +  JTs(x,y,z)  +  XIdb(x,y,z)  +  Ib(x,y,z) 

The  forward  component  (If)  refer  to  the  direction  of  propagation  of  the  main  fluence  beam  (or  ‘forward 
scatter’),  while  the  backward  component  (lb)  refers  to  the  180  degree  scatter  (or  backscatter)  which  is  opposite 
to  the  incident  photon  fluence  vector.  The  diagonal  components  are  divided  into  4  forward  diagonals  (Id)  and  4 
backward  diagonals  (Idb).  Thus  each  direction  is  associated  with  a  neighboring  voxel.  The  total  fluence  within 
a  voxel  is  obtained  by  summing  up  the  fluences  in  each  component  directions.  Note  that  the  above  directions 
can  be  arbitrarily  chosen  according  to  the  experimental  setup  or  at  the  convenience  of  the  user. 

b)  Simulate  photon  absorption  by  depositing  the  photon  fluence  into  the  voxel  based  on  its  absorption  coefficient 

(Ma): 

Absorbed  Fluence  =  Iabs(x,y,z)  =  I(x,y,z)  x  { l-exp(-pa  x  dl) } 

(where  dl  =  path  length) 

4)  Once  fluences  have  been  accumulated  over  the  current  layer,  increment  z  to  proceed  to  the  next  layer.  Thus  by 
repeating  this  process  for  all  layers  and  updating  the  directional  components  of  fluences  for  every  voxel,  the 
empirical  approach  propagates  the  fluence  from  one  layer  to  the  next  across  the  volume. 

2.2  Optimization  Phase  (Pencil  beam  simulations) 

Optical  properties  of  the  medium  (scattering  coefficient  ps,  anisotropy  factor  g,  absorption  coefficient  pa,  albedo  a’) 
define  the  photon  fluence  generated  by  the  empirical  approach.  The  weights  assigned  to  each  voxel  depend  on  its 
scattering  property,  while  the  absorption  coefficient  determines  the  amount  of  fluence  deposited  and  energy  absorbed. 
The  weights  were  derived  during  the  optimization  phase  by  comparing  the  exponential  attenuation  profiles  of  the 
Empirical  approach  to  the  Monte  Carlo  using  root  mean  square  error  and  coefficient  of  correlation  criteria.  Previous 
studies  have  shown  the  applications  as  well  as  the  limitations  of  the  diffusion  equation  in  estimating  photon  fluence  for 
media  with  low  albedos.10,15'22  In  this  study  we  first  tested  the  performance  of  the  Empirical  approach  with  respect  to  the 
Monte  Carlo  and  the  diffusion  equation  for  media  with  different  albedos  from  0.7  to  0.99.  The  albedo  (a’)10  is  given  by 
the  equation: 


(l-g^s+Ma 

Next,  the  empirical  algorithm  was  trained  to  derive  the  forward  (wtl)  and  diagonal  (wt2)  weights  for  brain  tissues  such 
as  white  matter  (p^O.OOSmm1,  ps=55rnrrf g=0.85,  a’=0.9994),  gray  matter  (p^O.OSSmm1,  ps=70mm"',  g=0.965, 
a’=0.9859),  skull  bone  (pa=0. 024mm'1,  p^lS^mm'1,  g=0.9,  a’=0.9871)  and  tumor  (astrocytoma)  (p^O.Smm1, 
ps=500mm_1,  g=0.88,  a’=0.995).12'14  A  pencil  beam  was  simulated  using  the  Monte  Carlo  using  the  known  optical 
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properties  of  each  of  the  above  brain  tissues.12'14  The  Empirical  algorithm  was  then  trained  to  emulate  the  same 
exponential  decay  profile  by  incrementing  the  forward  scattering  weight  from  0  to  1,  while  the  diagonal  was  adjusted  to 
ensure  that  wtl  +  (4  x  wt2)  =1.  The  weights  derived  from  the  optimization  phase  were  used  for  the  broad  beam 
simulations  in  homogeneous  and  heterogeneous  phantoms. 

2.3  Broad  beam  simulations  in  Homogeneous  and  Heterogeneous  Brain  phantoms 

The  empirical  weights  obtained  from  the  pencil  beam  training  phase  were  used  for  further  simulations  using  broad 
beams  in  homogeneous  and  heterogeneous  brain  phantoms.  The  Gaussian  beam  generated  by  the  Monte  Carlo  was 
matched  to  that  of  the  Empirical  algorithm.  An  optical  property  map  was  generated  in  which  every  voxel  was  assigned  a 
specific  weight  and  absorption  coefficient  based  on  the  tissue  property.  The  Gaussian  beam  distribution  was  first  tested 
in  homogenous  phantoms  of  each  of  the  brain  tissues:  skull  bone,  gray  matter,  white  matter  and  tumor  (astrocytoma). 
The  weights  derived  from  the  pencil  beam  optimization  phase  were  used  in  the  broad  beam  simulation  to  test  their 
validity.  Next,  a  software  phantom  of  a  human  head  with  15  cm  diameter  was  created  using  concentric  spheres  to  denote 
the  different  tissue  types.  The  incident  beam  had  a  Gaussian  profile.  The  fluence  profiles  and  absorption  energy  maps 
generated  Monte  Carlo  and  Empirical  approach  were  compared.  The  execution  times  for  generation  of  a  broad  beam 
profile  was  compared  between  the  two  models. 


3.  RESULTS 


3.1  Pencil  beam  simulations  (Optimization  Phase) 


In  this  study,  the  empirical  software  was  trained  in  an  iterative  process  by  comparing  its  exponential  decay  profiles  to 
that  of  the  Monte  Carlo  using  pencil  beams  in  optically  homogeneous  media.  The  Empirical  approach  had  a  better  fit  to 
the  Monte  Carlo  as  compared  to  the  diffusion  equation,  especially  for  low  albedos.  The  weights  wtl  and  wt2  were 
derived  and  optimized  for  each  of  the  brain  tissues  (white  matter,  gray  matter,  skull  bone  and  tumor)  by  using  the  profile 
which  best  matched  the  Monte  Carlo.  The  absorption  coefficient  assigned  was  the  same  as  has  been  obtained  from 
previous  studies.  Figure  3  shows  the  pencil  beam  simulation  profiles  obtained  by  the  Monte  Carlo  and  the  Empirical 
algorithm. 


|  | 
Pencil  Beam  Simulations  (Training  Phase) 


Gaussian  broad  beam  Simulation  (Exponential  profile 
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Figure  3.  (a)  Training  Phase.  The  figure  shows  pencil  beam  profiles  along  z  axis  between  the  Monte  Carlo  and  the 
Empirical  algorithm.  The  empirical  algorithm  is  trained  to  emulate  the  same  exponential  fluence  attenuation  as  that  of 
the  Monte  Carlo  by  adjusting  and  optimizing  the  weights  wtl  and  wt2.  (b)  Broad  beam  simulations:  The  empirical 
algorithm  was  compared  to  the  Monte  Carlo  for  incident  Gaussian  beam  profiles.  Overall  the  fluence  profile  of  the 
Empirical  closely  followed  that  of  the  Monte  Carlo. 
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Figure  3(b)  compares  the  Gaussian  beam  attenuation  along  beam  axis  between  the  Empirical  and  the  Monte  Carlo.  The 
Gaussian  beam  analysis  shows  that  the  empirical  approach  can  emulate  the  Monte  Carlo  fluence  for  broad  beam 
simulations.  The  acceleration  provided  by  the  Empirical  algorithm  was  approximately  700X  as  compared  to  the  Monte 
Carlo  for  a  broad  beam.  This  however  depends  on  the  tissue  optical  properties  as  the  Monte  Carlo  simulations  take 
longer  execution  times  for  high  scattering  tissues  and  are  relatively  fast  for  low  scattering  media.  The  differences 
between  the  Monte  Carlo  and  the  Empirical  approaches  can  be  significant  at  larger  depths  and  hence  need  to  be 
accounted  for. 
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3.3  Heterogeneous  Brain  phantom  simulations  (Gaussian  beam) 

Figure  5  shows  the  fluence  profiles  and  absorption  energy  maps  generated  Monte  Carlo  and  Empirical  approach  for  the 
head  phantom.  The  errors/differences  induced  by  the  tissue  heterogeneity  are  seen  between  the  two  models  and  needs  to 
be  studied  further.  However  it  is  seen  that  the  exponential  profile  of  the  Empirical  more  or  less  follows  that  of  the  Monte 
Carlo.  In  the  absorbed  energy  images  generated  in  Figure  5  (b)  and  (c),  the  profile  of  the  tumor  absorption  is  clearly 
visible  in  the  Empirical  and  can  be  demarcated  within  the  Monte  Carlo.  The  Monte  Carlo  also  generates  a  large  amount 
of  statistical  noise  at  larger  depths  which  is  not  seen  in  the  Empirical  model.  The  time  for  execution  of  the  Monte  Carlo 
for  a  broad  beam  in  the  heterogeneous  head  phantom  is  3.01  seconds  as  compared  to  9  hours  and  47  minutes  required  by 
the  Monte  Carlo.  Thus  the  Empirical  model  gives  a  smoother  approximation  of  the  fluence  profile  as  compared  to  the 
Monte  Carlo  at  a  fraction  of  the  time  (11701X  acceleration)  and  thus  could  provide  an  alternative  to  the  Monte  Carlo  in 
cases  where  an  approximate  idea  of  fluence  distribution  is  needed  while  planning  optical  therapy. 


4.  DISCUSSION 


4.1  Implications  and  Future  Research 

The  current  study  has  looked  at  providing  an  alternate  voxel  based  approach  of  estimating  photon  fluence  in  3D 
heterogeneous  tissues  besides  the  diffusion  approximation  and  Monte  Carlo  methods.  The  Monte  Carlo  is  by  far  the 
most  accurate  method  for  simulating  photon  deposition  in  tissues.  However  due  to  the  large  execution  times,  even  after 
being  implemented  on  parallel  processing  platforms,  there  is  a  need  for  a  method  for  fast  preliminary  estimation  of 
fluence  within  a  certain  error  limit.  The  empirical  approach  is  one  such  method  which  is  fast  and  can  approximate  the 
fluence  in  heterogeneous  tissues.  By  generating  fluence  in  a  3D  voxelated  format,  it  can  be  easily  integrated  with 
multiple  imaging  modalities.  More  study  is  however  necessary  to  understand  the  errors  induced  by  the  Empirical 
algorithm  at  large  depths  as  seen  in  the  pencil  beam  and  broad  beam  fluence  profiles.  These  errors  seen  can  be  large  at 
depths  greater  than  2cm  in  heterogeneous  phantoms.  Hence  the  empirical  approach,  at  this  time,  cannot  be  used  for 
absolute  quantification  of  photon  fluence.  It  can  however  give  a  good  approximation  of  the  relative  fluence  distribution 
in  tissues  with  respect  to  the  input  photon  fluence.  The  errors  in  quantification  can  be  reduced  by  further  research  into 
optimization  methods  to  match  the  Monte  Carlo.  The  training  methodology  used  to  derive  the  weights  could  be 
improved  to  account  for  the  differences  in  the  fluence  profiles.  Iterative  routines  and  algorithms  with  merit  functions  can 
significantly  enhance  the  accuracy  of  the  Empirical  approach. 

4.2  Conclusion 

The  ability  of  the  Empirical  approach  to  approximately  emulate  the  fluence  distribution  of  the  Monte  Carlo  has  been 
studied  in  homogeneous  and  heterogeneous  media.  Even  though  the  empirical  algorithm  cannot  be  used  for  absolute 
quantification,  it  can  be  used  in  optical  therapy  to  evaluate  the  merits  of  different  angles  of  illuminations  and  source 
positions  which  is  a  prerequisite  for  determining  the  placement  of  an  optical  source.  The  3D  GPU  based  Monte  Carlo 
can  then  be  used  for  absolute  quantification  and  validation  of  photon  energy  distribution.  The  acceleration  provided  by 
the  Empirical  algorithm  can  be  used  to  reduce  the  iterations  of  the  Monte  Carlo  by  identifying  only  a  few  possible 
illumination  possibilities.  Thus  the  Empirical  algorithm  has  the  potential  to  be  used  as  a  fast  preliminary  tool  to  devise 
optical  therapies  by  estimating  relative  fluence  distribution  in  tissues  prior  to  validation  and  correction  by  using  Monte 
Carlo. 
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An  Optical  Therapeutic  Protocol  to  treat  brain  metastasis 
by  mapping  NIR  activated  drug  release:  A  Pilot  Study 

Akshay  Prabhu  Verleker,  Qianqian  Fang,  Mi-Ran  Choi,  Susan  Clare  and  Keith  M.  Stantz 


Abstract-Treatment  of  Central  Nervous  System  (CNS) 
metastasis  poses  a  critical  clinical  challenge  due  to  limitations  in 
drug  uptake  across  the  blood  brain  barrier  and  blood- 
cerebrospinal  fluid  barrier.  Recent  research  has  shown  the 
efficacy  of  using  macrophages  as  drug  carriers  to  target 
metastatic  sites  in  the  brain,  which  can  then  be  activated  by 
illuminating  with  near  infrared  radiation.  The  goal  of  this 
research  is  to  develop  an  optically  targeted  therapeutic  treatment 
of  metastasis,  specifically  for  metastatic  breast  cancer  of  the 
brain.  As  a  first  step  towards  accomplishing  this,  we  developed  a 
3D  Monte  Carlo  photon  transport  code  capable  of  simulating 
phantoms  with  optical  properties  of  brain  and  tumor;  and 
validated  the  Monte  Carlo  simulated  photon  fluence  within  brain 
phantoms  using  an  optical  dosimetry  probe.  The  phantom  studies 
showed  good  correlation  (correlation  coefficient  R=0.977) 
between  the  probe  measurements  and  the  Monte  Carlo 
simulation  in  a  white  matter  phantom  (reduced  scattering 
coefficient  ps^S.lSmm"1,  absorption  coefficient  pa=0.005  mm'1). 
Our  future  steps  will  be  to  implement  the  Monte  Carlo  to  map 
out  photon  energy  distribution  in  the  brain,  and  subsequent  drug 
release,  by  segmenting  &  translating  head  CT  image  volumes  to 
corresponding  optical  properties  of  brain  tissues.  To  access  the 
therapeutic  response,  changes  in  the  vascular  physiology  of  the 
brain  due  to  Her2  inhibition  will  be  measured  using  dynamic 
contrast-enhanced  imaging  (e.g.  DCE-CT),  and  with  Monte  Carlo 
based  optical  fluence  maps.  An  optical  treatment  plan,  using  fast 
Monte  Carlo  software,  optimized  with  CT  segmented  image 
volumes,  would  significantly  reduce  the  treatment  time  and  allow 
targeted  drug  activation  while  sparing  healthy  tissues. 

I.  Introduction 

reatment  of  Central  Nervous  System  (CNS)  metastasis 
poses  a  critical  clinical  challenge  due  to  limitations  in  drug 
uptake  in  the  brain  (across  the  blood  brain  barrier  and  blood- 
cerebrospinal  fluid  barrier)  and  adverse  neurotoxic  effects  of 
mainstay  therapies  such  as  whole  brain  radiation  therapy 
(WBRT)  and  stereotactic  radiosurgery.1  With  non-CNS 
metastasis  being  treated  successfully  with  receptor  targeted 
drug  delivery,  aided  by  improved  imaging  and  localized 
radiotherapy  techniques,  the  brain  remains  a  sanctuary  for 
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metastasis  (e.g.  triple  negative  and  HER2+  breast  cancer 
metastasis  with  median  survival  rates  of  2  to  16  months).1  This 
study  is  part  of  a  project  which  demonstrated  a  macrophage 
based  “Trojan  Horse”  delivery  of  fluorescent  molecules  into 
brain  metastasis2  and  subsequent  release  of  the  fluorescent 
molecules  on  activation  in-vivo  by  NIR  radiation.3  The  larger 
goal  of  this  collaborative  project  is  to  deliver  lapatinib-gold 
nanocomplex  laden  macrophages  to  brain  metastasis  and 
design  a  therapeutic  protocol  to  optically  stimulate  drug 
release  in  target  tissues  in  the  brain.  In  the  current  project,  the 
purpose  is  to  design  an  optical  simulation  protocol  (using 
Monte  Carlo  simulations)  to  effectively  predict  the  photon 
distribution  and  subsequent  drug  activation  in  the  brain  tissue, 
and  validate  it  by  using  an  optical  dosimetry  probe.  With  the 
rate/quantity  of  drug  release  being  directly  dependent  on  the 
optical  power  delivered,  the  Monte  Carlo  simulation  protocol, 
with  the  ability  to  predict  photon  distribution  in  heterogeneous 
brain  tissue,  can  be  used  to  generate  an  equivalent  drug  release 
map.  The  optical  properties  of  the  brain  can  be  mapped  by 
converting  the  CT  image  densities  of  the  brain  into  the 
corresponding  optical  properties  of  the  tissue  (e.g.  white/gray 
matter,  CSF,  etc.),  which  would  serve  as  an  input  to  the  Monte 
Carlo  routine.  The  optical  drug  release  map  generated  by 
Monte  Carlo  can  be  co-registered  with  CT  images  to  improve 
the  therapeutic  efficacy  of  the  protocol. 

In  order  to  predict  the  photon  energy  distribution  in 
heterogeneous  brain  tissue,  we  used  a  GPU  based  3D  Monte 
Carlo  simulation  code4,  which  was  modified  to  generate  user 
defined  beam  profiles  to  emulate  the  NIR  (800  nm)  laser 
source  used  in  our  lab.  Monte  Carlo  simulation  has  been  used 
as  the  gold  standard  for  photon  propagation  studies  (in  3D 
heterogeneous  media),  with  the  GPU  version  having 
acceleration  speeds  100  to  300  times  faster  than  corresponding 
CPU  versions.4  In  order  to  validate  and  quantify  the  Monte 
Carlo  generated  fluence  in  brain  tissue  phantoms,  we  used  an 
optical  dosimetry  probe  to  measure  the  photon  fluence  in  low 
optical  fluence  regions.  An  optimized  treatment  plan  using  a 
fast  Monte  Carlo  software  to  predict  in  vivo  drug  release, 
would  significantly  reduce  the  treatment  time  and  allow  the 
targeted  drug  activation  while  sparing  healthy  tissues. 

In  addition  to  designing  the  Optical  therapeutic  protocol,  we 
also  looked  into  the  feasibility  of  using  Dynamic  Contrast 
Enhanced  CT  (DCE-CT)  imaging  to  monitor  and  quantify  the 
efficacy  of  treatment.  The  vascular  structure  of  a  tumor  is 
complex  and  heterogeneous,  and  is  seen  to  vary  between 
individual  tumors  and  during  different  phases  of  tumor 
growth.5  This  heterogeneity  in  tumor  vasculature, 
characterized  by  hypoxia,  anoxia  and  hypoglycemia,  leads  to 
differing  response  to  treatment  due  to  formation  of  barriers  to 
transport  and  distribution  of  drug  molecules.5’9  Following 
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treatment  with  anti-angiogenic  or  cytotoxic  drugs,  the  tumor 
vasculature  undergoes  complex  physiological  and  biological 
changes  some  of  which  are  adaptive  in  nature.5’9  DCE-CT  has 
been  successfully  used  to  simultaneously  quantify 
physiological  parameters  such  as  perfusion  and  vascular 
plasma  volume  in  tumors  after  treatment  with  therapeutic 
drugs.10'12  Past  research  has  proved  the  feasibility  of  using 
DCE-CT  in  measuring  the  physiological  response  of  breast 
tumors  to  anti-angiogenic  therapies  using  the  two- 
compartmental  model.11’12  In  the  present  study  we  used  the 
above  model  to  obtain  parametric  maps  of  perfusion, 
permeability,  fractional  plasma  volume  and  fractional 
interstitial  volume  in  glioma  tumors  grown  in  mice;  thus 
demonstrating  the  feasibility  of  DCE-CT  imaging  to  monitor 
brain  tumor  response  to  treatment. 

II.  Methods 

Step  1:  The  optical  dosimetry  probe  (1.5mm  diameter)  was 
calibrated  using  a  calorimeter  (Gold  standard)  in  low  fluence 
regions.  The  probe  measurements  were  used  to  validate  Monte 
Carlo  generated  photon  fluence  in  brain  phantoms  resembling 
white  and  gray  matter,  made  with  specific  concentrations  of 
India  ink  (absorber)  and  Intralipid  (scatter),13'14  as  measured 
by  photo-spectrometer  (Gold  standard),  consistent  with  brain 
tissue  and  tumors  (values)  as  shown  in  Fig  1.  The  source  of 
light  used  was  an  integrating  sphere  connected  to  a  pulsed 
NIR  laser  with  wavelength  800  nm,  frequency  200MHz,  and 
energy  lOmJ/pulse.  The  broad  beam  NIR  output  power 
emanating  from  the  integrating  sphere  was  0.2083mW/mm2. 

Step  2:  In  order  to  apply  the  Monte  Carlo  (MC)  simulation 
protocol  to  the  brain,  we  designed  a  segmentation  software  to 
determine  tissue  boundaries  in  the  brain  using  Computed 
tomography  images.  The  CT  image  densities  in  Hounsfield 
units  (HU)  of  the  brain  were  translated  into  corresponding 
optical  properties  (absorption  coefficient,  scattering 
coefficient,  and  anisotropy  factor)  using  lookup  tables  (Fig  2). 
The  lookup  tables  were  prepared  from  earlier  studies  which 
derived  the  Hounsfield  units  and  optical  properties  of  brain 
tissues  such  as  white/gray  matter.13'15  The  segmentation 
software  generated  an  optical  map  of  the  brain,  which  was 
given  as  input  to  the  Monte  Carlo.  The  photon  energy  map 
generated  by  the  Monte  Carlo  can  be  co-registered  with  the 
CT  image  to  visualize  the  efficacy  of  drug  release. 

Step  3:  “Optical  treatment  protocol”,  voxel-based  optical 
properties  of  a  human  head  (obtained  from  CT  scans)  were 
input  into  MC  code  to  simulate  3D  fluence  maps  from 
multiple  (optically  coupled)  light  sources.  This  study  also 
looked  into  the  feasibility  of  monitoring  effectiveness  of 
targeted  drug  release  by  monitoring  the  physiology  at  target 
sites  and  surrounding  tissues.  An  Orthotopic  model  of 
glioblastoma  tumors  were  grown  in  mice  and  parametric 
images  of  perfusion,  permeability,  fractional  plasma  volume 
and  fractional  interstitial  volume  were  derived  using  the  two 
compartmental  model.12 


Probe  measurements  v/s  Monte  Carlo 
(white  matter  phantom) 

1  1 

4> 

2  I 

S  v 

QJ 

4 

♦  Peak  area  normalized 

© 

4 

► 

c 

►+ 

(Integrating  Sphere) 

© 

(WM  lip  phantom) 

.= 

Qh 

o  ni 

♦, 

* 

4 

►+ 

N  U.l 

+  MC  broad  beam 

s 

© 

Normalized 

£ 

< 

H- 

u 

© 

(Plexiglas  s2) 

S3 

OX 

© 

n  ni 

Distance  (cm) 

Fig.  1.  Validation  of  Monte  Carlo  simulated  photon  fluence  (l/mm2)  with 
probe  measurements  (peak  area  in  mV)  in  a  white  matter  phantom.  Log 
normalized  fluence  plotted  v/s  depth  in  a  white  matter  phantom.  Correlation 
R=l. 


III.  Observations  &  Results 


The  Monte  Carlo  generated  photon  fluence  was  validated 
by  probe  measurements  in  optical  brain  phantoms  (e.g.  white 
matter  in  Figure  1).  The  optical  fluence  measured  by  the  probe 
was  in  turn,  calibrated  using  a  calorimeter  (Gold  standard)  and 
showed  a  linear  and  isotropic  response  to  input  NIR  fluence  at 
different  incidence  angles  (<5%  variation).  Our  MC 
simulation  studies  showed  that  99%  of  photons  are  attenuated 
by  the  bone  (1cm),  which  was  confirmed  by  probe 
measurements  in  an  empty  goat  skull.  The  segmentation  code 
(with  iterative  Monte  Carlo)  can  be  used  to  identify  regions  of 
least  bone  thickness  and  determine  the  path  of  least 
attenuation  to  couple  maximum  power  to  brain  tissues.  This 
study  proved  that  an  optimized  Monte  Carlo  could  be 
successfully  used  to  quantify  photon  propagation  in  the  brain 
and  thus  predict  the  rate  of  drug  release  in  target  tissues  to 
depths  up  to  5cm  in  brain  tissues. 
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Fig.  2.  Monte  Carlo  based  treatment  protocol.  Head  CT  image  densities  (in 
Hounsfield  Units/HU)  are  segmented  and  translated  into  optical  properties  of 
the  brain  tissue  (white/gray  matter,  CSF,  skull  bone,  etc.).  Monte  Carlo 
routine  is  then  used  to  simulate  light  propagation  through  the  brain.  [CT 
image  reference:  http://www.mpoullis.com/thorcd/rad/23.htm]16 
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Our  preliminary  longitudinal  study  has  indicated  that  DCE-CT 
is  an  effective  tool  to  study  tumor  progression  and  physiology 
(Fig  3)  of  glioma  tumors.  We  propose  to  monitor  these  effects 
by  studying  the  change  in  fractional  plasma  volume  and 
fractional  interstitial  volume,  both  of  which  are  excellent 
indicators  of  vascular  physiology  (Fig  3). 
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Fig.  3.  Monitoring  cancer  therapy  through  Dynamic  Contrast  Enhanced  CT 
(DCE-CT).  (A)  Schematic  of  two  compartmental  model.12  (B)  Parametric 
images  of  glioblastoma  (orthotopic  model)  grown  in  mice.  In  this  longitudinal 
study,  parametric  images  were  obtained  of  perfusion,  permeability,  fractional 
plasma  volume  and  fractional  interstitial  volume  of  glioma  tumors  grown  in 
mice  using  the  two  compartmental  model.  These  parameters  can  be  used  to 
quantify  the  efficacy  of  treatment,  thus  proving  effectiveness  of  using  DCE- 
CT  to  monitor  targeted  dmg  therapy. 


IV.  Discussion 

In  this  study,  we  have  assumed  that  the  quantity  of  dmg 
release  is  proportional  to  the  optical  power  coupled  to  the 
tissues  (as  shown  in  previous  in-vitro  studies).  Future  studies 
will  focus  on  quantifying  the  optical  dmg  release  threshold  by 
determining  the  relationship  between  optical  power  coupled 
and  amount  of  drug  released.  Further  studies  will  involve 
adjusting  beam  profiles  to  optimize  the  coupling  and 
localization  of  optical  energy  in  the  brain  for  maximal  dmg 
release.  In  order  to  monitor  tumor  response  to  therapy,  DCE- 
CT  will  be  used  to  further  study  vascular  normalization17  in 
tumors  caused  by  anti-angiogenic  effects  of  lapatinib  dmg, 
which  can  help  in  optimizing  therapeutic  efficacy. 


V.  Conclusion 

This  study  has  successfully  demonstrated  the  potential  of 
using  a  fast  voxel  based  Monte  Carlo  based  routine  to  predict 
the  rate  of  dmg  release  in  heterogeneous  brain  tissue  by  using 
Computed  Tomography  images  to  segment  and  translate  tissue 
boundaries  of  the  brain.  An  optimized  therapeutic  protocol 
with  a  fast  Monte  Carlo  routine  can  be  used  to  determine  the 
best  irradiation  conditions  for  targeted  dmg  release  in  brain 
tissue  within  clinically  relevant  time  frames. 
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A  three-dimensional  photon  dosimetry  in  tissues  is  critical  in  designing  optical  therapeutic  protocols  to  trigger 
light-activated  drug  release.  The  objective  of  this  study  is  to  investigate  the  feasibility  of  a  Monte  Carlo-based 
optical  therapy  planning  software  by  developing  dosimetry  tools  to  characterize  and  cross-validate  the  local  pho¬ 
ton  fluence  in  brain  tissue,  as  part  of  a  long-term  strategy  to  quantify  the  effects  of  photoactivated  drug  release  in 
brain  tumors.  An  existing  GPU-based  3D  Monte  Carlo  (MC)  code  was  modified  to  simulate  near-infrared  photon 
transport  with  differing  laser  beam  profiles  within  phantoms  of  skull  bone  (B),  white  matter  (WM),  and  gray 
matter  (GM).  A  novel  titanium-based  optical  dosimetry  probe  with  isotropic  acceptance  was  used  to  validate  the 
local  photon  fluence,  and  an  empirical  model  of  photon  transport  was  developed  to  significantly  decrease  ex¬ 
ecution  time  for  clinical  application.  Comparisons  between  the  MC  and  the  dosimetry  probe  measurements  were 
on  an  average  11.27%,  13.25%,  and  11.81%  along  the  illumination  beam  axis,  and  9.4%,  12.06%,  8.91% 
perpendicular  to  the  beam  axis  for  WM,  GM,  and  B  phantoms,  respectively.  For  a  heterogeneous  head  phantom, 
the  measured  %  errors  were  17.71%  and  18.04%  along  and  perpendicular  to  beam  axis.  The  empirical  algorithm 
was  validated  by  probe  measurements  and  matched  the  MC  results  (R2  >  0.99),  with  average  %  error  of  10.1%, 
45.2%,  and  22.1%  relative  to  probe  measurements,  and  22.6%,  35.8%,  and  21.9%  relative  to  the  MC,  for  WM, 
GM,  and  B  phantoms,  respectively.  The  simulation  time  for  the  empirical  model  was  6  s  versus  8  h  for  the  GPU- 
based  Monte  Carlo  for  a  head  phantom  simulation.  These  tools  provide  the  capability  to  develop  and  optimize 
treatment  plans  for  optimal  release  of  pharmaceuticals  in  the  treatment  of  cancer.  Future  work  will  test  and 
validate  these  novel  delivery  and  release  mechanisms  in  vivo .  ©  2016  Optical  Society  of  America 

OC/S  codes:  (170.0170)  Medical  optics  and  biotechnology;  (170.3660)  Light  propagation  in  tissues;  (170.3890)  Medical  optics 
instrumentation;  (260.3060)  Infrared;  (040.0040)  Detectors. 

https://doi.org/10.1364/AO.55.009875 


1.  INTRODUCTION 

Optical  therapy  is  a  non-invasive  technique  that  can  be  uti¬ 
lized  to  treat  cancer  through  photodynamic  therapy  (PDT) 
techniques  [1]  and  the  triggered  release  of  drugs  from  nano¬ 
complexes  in  metastatic  tissues  [2,3]-  Near-infrared  (NIR)  light 
is  preferred  for  deep  tissue  activation  due  to  its  low  optical  at¬ 
tenuation  and  thus  high  penetration.  This  study  is  part  of  a 
project  which  seeks  to  develop  optically  activated  drug  release 
in  order  to  implement  macrophage-based  “Trojan  horse”  drug 
delivery  to  brain  metastasis  [2,3]-  Given  that  the  rate  of  release 
and  quantity  of  the  released  drug  are  likely  to  be  affected  by  the 
optical  power  delivered  [2,3] ,  the  development  of  optical 


simulation  and  dosimetry  tools  that  can  accurately  determine 
the  photon  fluence  in  the  brain  and  ultimately  be  integrated 
into  a  clinical  environment  is  critical.  In  our  study  a  GPU- 
based  3D  Monte  Carlo  (MC)  code  is  used  to  simulate  the  pho¬ 
ton  energy  deposition  in  the  brain  [4-10],  and  an  optical 
dosimetry  probe  was  fabricated  to  validate  the  photon  fluence 
in  heterogeneous  media  consistent  with  the  skull  and  brain. 
While  the  primary  application  is  in  the  field  of  controlled  drug 
release,  these  tools  can  also  be  applied  to  a  wide  range  of 
applications,  including  NIR-based  imaging  techniques,  such 
as  photoacoustic  tomography [1 1]  and  near-infrared  spectros¬ 
copy  [12],  and  in  infrared  neural  stimulation  [13]  of  the  coch¬ 
lea  and  neuroprosthesis. 
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Monte  Carlo  simulation  methods  have  been  used  as  the  gold 
standard  for  photon  propagation  studies  in  heterogeneous  me¬ 
dia  [4-8].  Traditional  CPU-based  Monte  Carlo  models  suffer 
from  slow  computation  speeds,  which  makes  them  non- 
viable  for  clinical  applications,  which  require  the  simulation 
of  a  large  number  of  photons  (excess  of  1 07 )  especially  in  highly 
attenuating  (scattering)  brain  tissues,  such  as  the  skull  and 
white  matter.  The  GPU-based  3D  Monte  Carlo  used  in  this 
study,  developed  by  the  Fang  group  [9],  can  achieve  complete 
brain  simulation  of  a  high-fluence  broad-beam  source  within 
hours  compared  to  days  for  CPU-based  models  (300  x  faster) 
[9].  This  new  model  features  a  voxel-based  3D  fluence  and 
energy  distribution  with  significant  improvements,  such  as 
efficient  parallel  random  number  generators,  reflection  at 
3D  boundaries,  time-resolved  simulations,  fluence  normaliza¬ 
tions,  and  ease  of  integration  with  anatomical  imaging  modal¬ 
ities  (e.g.,  CT,  MR,  ultrasound,  etc.),  thus  improving  its 
clinical  viability  and  efficacy  for  diagnosis  and  therapy.  To  date, 
this  code  has  not  been  validated  through  3D  measurements 
under  various  illumination  conditions. 

Application  of  a  3D  Monte  Carlo  as  an  optical  therapy 
planning  tool  in  clinical  practice  requires  rigorous  validation 
of  its  physics-based  model,  so  as  to  reduce  the  uncertainty 
in  light  dose  delivery  and  improve  diagnostic  accuracy.  Past 
studies  have  validated  the  Monte  Carlo  using  diffusion 
approximation  simulations  in  slab  geometry  and  under  limited 
heterogeneity  conditions  [8,9].  Only  a  few  studies  have  at¬ 
tempted  to  rigorously  validate  the  3D  MC  fluence  using 
dosimetry  measurements  in  phantoms  under  varying  illumina¬ 
tions  with  finite  beam  width  and  heterogeneity  conditions 
[14,13].  The  varying  degree  of  accuracy  can  be  attributed 
to  differences  in  estimation  of  optical  properties  of  the 
medium,  illumination  and  object  boundary  conditions,  and 
the  ability  of  the  Monte  Carlo  to  accurately  estimate  photon 
distribution  in  complex  heterogeneous  boundary  conditions. 
All  of  these  differences  need  to  be  analyzed  and  optimized 
in  order  to  get  maximum  efficiency  from  the  Monte  Carlo 
simulations.  The  use  of  optical  dosimetry  for  pre-clinical 
validation,  along  with  the  ability  to  verify  optical  properties 
of  phantoms  used  in  the  measurements,  is  essential  in  order 
to  ensure  therapeutic  accuracy  of  the  Monte  Carlo. 
Additionally,  any  mismatch  in  simulation  and  actual  measure¬ 
ment  conditions  can  be  addressed  by  this  exercise  to  improve 
therapeutic  confidence  in  the  tool. 

To  validate  the  MC,  the  localized  photon  distribution  in 
tissues  was  measured  using  a  novel  titanium  dioxide-based  op¬ 
tical  dosimetry  probe  (ODP).  ODPs  have  been  extensively  used 
in  PDT  to  measure  photon  fluence  as  well  as  to  deliver  light 
locally  to  tissues  [1,16-21].  The  small  size  of  the  dosimetry 
probe  (1.5  mm  diameter)  and  near-isotropic  response  at 
NIR  wavelength  (800  nm)  improves  the  accuracy  of  localized 
light  measurement  in  highly  scattering  tissues/phantoms,  com¬ 
pared  to  other  dosimetry  probes.  It  is  also  advantageous  over 
standard  photon  measurement  devices,  such  as  calorimeters, 
which  only  measure  fluence  in  air.  In  this  study,  the  novel 
titanium-based  probe  was  designed  and  evaluated  against  a 
standard  nylon-based  probe  [17,20]  and  their  linear  and 
isotropic  response  analyzed.  Its  linear  response  in  air  was 


calibrated  using  a  calorimeter  (gold  standard)  for  absolute  flu¬ 
ence  quantification. 

Translating  an  optical  therapy  planning  tool  into  the  clinic 
requires  that  simulations  be  performed  in  feasible  time  frames. 
While  the  GPU-based  Monte  Carlo  takes  a  few  hours  for  a 
complete  head  simulation  under  a  single  illumination  condi¬ 
tion,  optimizing  a  treatment  plan  would  require  iterations  with 
multiple  laser  sources.  While  a  cluster  of  GPUs  could  be  used  to 
solve  this  problem,  thereby  significantly  increasing  costs,  it 
prompted  us  to  design  an  empirical  model  of  photon  propaga¬ 
tion  [22]  to  approximate  3D  photon  fluence  in  tissue-like  me¬ 
dia  at  a  fraction  of  time  taken  by  the  GPU-based  Monte  Carlo. 
This  model  assumes  the  photon  fluence  in  a  voxel  is  the 
weighted  sum  of  the  fluences  in  the  neighboring  voxels  in 
the  preceding  layer  [22] ,  where  the  weights  are  determined  us¬ 
ing  the  Monte  Carlo  for  typical  brain  tissues,  such  as  gray  mat¬ 
ter,  white  matter,  skull  bone,  and  astrocytoma.  Previous  studies 
have  used  approximate  methods,  such  as  the  diffusion  approxi¬ 
mation  (DE)  [23,24],  hybrid  (MC  +  DE)  [25],  radiative  trans¬ 
port  equation  [26,27],  and  the  adding-doubling  method  [28], 
to  approximate  3D  photon  distribution.  The  empirical  method 
(EM)  is  an  addition  to  the  existing  list  of  methods  and  mainly 
differs  in  the  means  by  which  photon  scatter  is  calculated  and 
has  a  similar  speed  as  the  diffusion  approximation  [29]  (in  sec¬ 
onds).  While  other  methods  rely  on  physics-based  models 
based  on  approximations  to  the  radiative  transport  equation, 
the  empirical  approach  is  designed  to  be  calibrated  against 
the  Monte  Carlo  in  a  homogeneous  medium  for  a  fixed  set 
of  optical  properties,  before  being  applied  in  a  complex 
medium  under  different  illumination  conditions.  Hence  a 
calibration  step  with  a  fixed  set  of  optical  properties  is  necessary 
before  applying  this  model.  It  is  important  to  note  that  once 
calibrated,  this  model  does  not  require  additional  calibrations, 
which  allows  for  its  repetitive  use  for  standard  clinical  applica¬ 
tions  with  a  fixed  set  of  optical  properties  (e.g.,  head  simula¬ 
tion)  under  differing  anatomical  geometries  and  illumination 
conditions.  The  significant  improvement  in  the  speed  and 
reasonable  accuracy  inspired  us  to  integrate  it  into  the  optical 
therapy  planning  tool  along  with  the  Monte  Carlo  and  the 
titanium-based  optical  dosimetry  probe. 

In  this  paper,  a  titanium-based  optical  dosimetry  probe  was 
fabricated  and  calibrated  to  provide  localized  energy  fluence 
measurements,  and  used  to  validate  the  3D  GPU-based 
Monte  Carlo  code  under  differing  beam  profiles  and  hetero¬ 
geneous  tissue  phantoms.  To  provide  an  accelerated,  real-time 
solution  to  photon  transport,  an  empirical  model  was  proposed 
and  evaluated  against  MC  simulated  in  tissue  phantoms. 
Finally,  the  role  of  these  new  and/or  improved  tools  in  optically 
stimulated  drug  release  will  be  discussed. 

2.  MATERIALS  AND  METHODS 
A.  Optical  Brain  Phantoms 

Brain  phantoms  were  designed  to  mimic  the  optical 
absorption  and  reduced  scattering  coefficients  of  white  matter, 
gray  matter,  and  skull  bone  (at  800  nm)  [29-31]  using  India 
ink  (absorber)  and  intralipid  (scatterer)  [32-35]  in  water 
(fia  =  0.020  cm-1)  [36].  The  absorption  coefficient  of  India 
ink  (Higgins  non-water-proof  black  ink)  and  the  scattering 
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Intralipid  Calibration:  Scattering  &  Reduced  Scattering  Coefficient  v/s 
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Fig.  1.  Calibration  curve  of  spectrophotometer  measured  absorb¬ 
ance  versus  %  concentration  of  India  ink.  The  error  bars  represent 
minimum  and  maximum  absorbance. 


coefficient  of  intralipid  (Intralipid  20%  from  Fresenius 
Kabi),  at  800  nm  wavelength,  were  determined  and  validated 
using  a  spectrophotometer  (Shimadzu  UVmini  1240 
Spectrophotometer),  which  served  as  the  gold  standard  for 
absorbance.  The  scattering  properties  of  intralipid  were  deter¬ 
mined  using  the  spectrophotometer  for  low  concentrations 
(0.003%-0.012%),  where  \is  =  9.32  cm“1%concentration~1, 
compared  to  9.95  cm_1%concentration_1  as  shown  in  previous 
studies  by  Van  Staveren  [33],  which  is  within  the  4.32%  rel¬ 
ative  error.  Figures  1  and  2  show  the  calibration  curves  of  ab¬ 
sorbance  (or  attenuation)  versus  concentration  of  India  ink  and 
intralipid  obtained  using  the  spectrophotometer.  The  linear  re¬ 
lationship  between  concentration  and  absorbance  was  used  to 
design  optical  phantoms  with  specific  absorption  and  scattering 
properties  of  white  matter,  gray  matter,  and  skull  bone  [29,30]. 
These  optical  properties  (absorption  and  scattering  coefficients 
and  anisotropy  factor)  are  shown  in  Table  1  and  were  used  to 
design  the  optical  brain  phantoms  (Table  2)  using  the  equa¬ 
tions  derived  from  Figs.  1  and  2. 

The  linear  relationship  of  concentration  versus  total 
attenuation  is  particularly  useful  in  determining  intralipid 
concentrations  for  high  scattering  tissues  in  the  brain.  Both 
homogeneous  as  well  as  heterogeneous  phantoms  were  made 
using  various  concentrations  of  India  ink  and  intralipid.  The 
phantom  composition  for  the  homogeneous  phantom  was  a 
solution  of  India  ink  and  intralipid  in  water,  while  the  hetero¬ 
geneous  phantoms  were  comprised  of  layers  of  1%  agar 
(Fisher  Scientific  Laboratory  Grade  Agar  A3600-500)  doped 
with  India  ink  and  intralipid  to  represent  skull  bone  and  gray 
matter,  with  a  layer  of  white  matter  in  aqueous  solution.  The 
scattering  properties  of  agar  and  the  absorption  properties  of 
India  ink  were  also  measured  using  the  spectrophotometer 
for  each  individual  sample  used  in  the  study  while  the  scattering 
properties  of  intralipid  were  derived  from  the  linear  relationship 
in  Fig.  2. 


Intralipid  Calibration:  Extrapolated 
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Fig.  2.  Intralipid  scattering:  upper  figure  represents  measured  values 
of  scattering  coefficients  (jis  and  fi) )  versus  %  concentration  of  intra¬ 
lipid  using  a  spectrophotometer.  The  error  bars  represent  minimum 
and  maximum  values.  Lower  figure  shows  extrapolated  values  from 
the  linear  relationship  derived  from  the  measured  values  (with  relative 
and  statistical  uncertainties  of  4.32%  and  1.8%,  respectively). 


B.  Optical  Dosimetry  Probe  Design  and  Calibration 

The  optical  dosimetry  probe  consists  of  a  spherical  bulb 
of  highly  scattering  material  attached  to  an  optical  fiber 
[1,16,17,19].  The  construction  and  fluence  measurement  by 

Table  1 .  Optical  Properties  of  Brain  Tissues  [29,30] 

Reduced 

Absorption  Scattering  Scattering 

Tissue  Coefficient  Coefficient  Anisotropy  Coefficient 


Type 

Ha  (cm"1) 

Hs  (cm"1) 

Factor  g 

H's  (cm-1) 

White 

0.03 

550 

0.85 

82.5 

matter 

Gray 

0.35 

700 

0.965 

24.5 

matter 

Skull 

0.24 

184 

0.9 

18.4 

bone 
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Table  2.  Phantom  Composition  for  MC  Versus  Probe  Comparison 

Phantom  Type 

Tissue  Type 

Intralipid  % 

Ink  % 

Absorption 
Coefficient  fia  (cm1) 

Reduced  Scattering 
Coefficient  fi's  (cm-1) 

Liquid  phantom^ 

White  matter 

7.3 

0 

0.002 

71.84 

Liquid  phantom22 

Gray  matter 

2.36 

0.023 

0.379 

22.64 

Liquid  phantom'2 

Skull  bone 

1.82 

0.0156 

0.26 

17.44 

Solid  phantonf 

Gray  matter 

2.36 

0.0195 

0.35 

24.5 

Solid  phantonT 

Skull  bone 

1.92 

0.0134 

0.25 

18.399 

Thantom  composition  used  in  homogeneous  phantom  studies. 

^Phantom  composition  used  in  heterogeneous  (layered)  phantom  studies. 

Notes:  The  anisotropy  factor  of  intralipid  (g  =  0.636)  was  used  for  the  phantoms. 


the  probe  are  shown  in  Fig.  3.  The  optical  fiber  is  connected 
to  a  PIN  photodiode  circuit  (ET-2030,  Electro  Optics 
Technology  Inc.)  whose  output  voltage  is  read  using  an  oscillo¬ 
scope  (Tektronix  TDS3052B  Digital  Oscilloscope).  Two  mod¬ 
els  of  optical  dosimetry  probes  with  tips  made  of  nylon  and 
titanium  dioxide  were  designed  [37].  The  performance  of 
the  traditionally  used  nylon-based  probes  [19]  was  compared 
to  that  of  the  titanium  dioxide-based  probes  by  measuring 
the  response  to  different  angles  of  incidence  in  the  equatorial 
and  azimuthal  planes,  e.g.,  isotropicity  response  function. 

The  nylon  spheres  have  been  widely  used  for  photodynamic 
therapy  and  were  made  of  nylon  due  to  their  isotropic  response 
[19].  These  were  used  as  the  standard  of  reference  to  compare 
the  performance  of  the  novel  titanium-based  probes.  The  tita¬ 
nium-based  probes  have  the  scattering  tips  made  of  a  mixture 
titanium  dioxide  (Du  Pont  Ti-Pure  R-900)  and  a  clear  two-part 
epoxy  (Tra-Con  BA-F1 14)  mixed  in  a  ratio  of  9.1  mg  Ti02  to 
1  ml  of  the  epoxy  (mixed)  [37,38].  The  use  of  epoxy  as  a 
medium  for  Ti02  is  ideal  as  it  is  transparent  to  NIR  light 
and  ensures  that  light  scatter  in  the  probe  tip  is  entirely  due 
to  Ti02  particles.  The  titanium  probes  were  selected  since 
the  particle  size  of  titanium  dioxide  is  around  410  nm,  which 
results  in  optimal  scattering  at  820  nm  or  when  the  wavelength 
is  approximately  twice  the  particle  size.  Since  the  wavelength 
used  in  our  studies  is  800  nm,  the  titanium  probes  have  been 
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Fig.  3.  Design  and  working  of  the  optical  dosimetry  probe. 


designed  to  show  a  higher  sensitivity  than  the  nylon  probes. 
The  titanium-based  probe  tips  were  molded  within  a  specially 
designed  caste  and  a  fiber  optic  cable  (BFL48-400,  4  mm  diam¬ 
eter,  0.48  numerical  aperture)  was  inserted  within  the  cured 
spheres  at  a  depth  of  approximately  1/3  of  the  sphere’s  radius 
[37] .  The  nylon-based  probes  were  approximately  3. 173  mm  in 
diameter  and  were  machined  in  a  multistage  process  [37].  The 
other  end  of  the  fiber  optic  cable  was  connected  to  a  PIN  pho¬ 
todiode  (ET-2030,  Electro  Optics  Technology  Inc.)  through  a 
SMA-905  connector  and  the  voltage  output  was  read  using  an 
oscilloscope  (Tektronix  TDS3052B). 

Figures  4  and  5  show  the  setup  used  to  test  the  isotropicity 
of  the  dosimetry  probes  in  the  equatorial  and  azimuthal  planes. 
The  probe  tip  is  illuminated  by  laser  light  emanating  from  an 
integrating  sphere  (Melles  Griot  two-port  integrating  sphere), 
which  produces  a  uniform  broad  beam  of  diffuse  light,  which  is 
used  to  characterize  the  probe  response.  The  probe  was  rotated 
with  angular  increments  of  15°  in  the  equatorial  and  10°  in  the 
azimuthal  planes,  and  its  response  on  the  oscilloscope  was 
measured. 

The  center  of  the  probe  tip  was  maintained  at  the  center  of 
rotation  during  the  entire  procedure.  The  isotropicity  of  the 
probe,  measured  by  calculating  its  coefficient  of  variation,  de¬ 
fined  as  the  ratio  of  standard  deviation  to  the  mean  of  the  rms 
value  over  all  the  angles,  is  shown  in  Fig.  6(a).  The  range  of 


Optical  dosimetry  probe 


Fig.  4.  Setup  for  calibration  of  the  optical  dosimetry  probe:  meas¬ 
uring  isotropicity  (in  equatorial  and  azimuthal  planes)  of  the  probe. 
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Optical  fiberguide 
(connected  to  a 
laser) 


Fig.  5.  Setup  for  calibration  of  the  optical  dosimetry  probe: 
measuring  the  linearity  of  the  probe  to  incident  laser  fluence. 


angles  was  0  to  360  deg  for  equatorial  measurements  and  0  to 
130  deg  for  azimuthal  measurements.  The  angular  span  for  azi¬ 
muthal  measurements  was  limited  due  to  physical  limitations  of 
the  experimental  setup.  The  individual  rms  voltage  measured 
by  the  probe  was  normalized  to  the  power  measured  using  a 
calorimeter. 

In  the  test  of  linearity  (Fig.  5),  the  dosimetry  probe  response 
is  validated  using  a  calorimeter  (Ophir  CE,  Nova  2  Ophir: 
3A-P-SH-VI,  aperture  12  mm)  by  illuminating  with  a  pulsed 
laser  beam  of  wavelength  800  nm,  frequency  10  Hz,  output 
fluence  17.64  mW/cm2.  The  probe  and  the  calorimeter  are 
translated  along  the  beam  axis  (in  air)  and  the  photon  fluence 
is  measured  at  different  positions.  The  linearity  of  response  of 
the  1.5  mm  titanium  probe  (selected  due  to  its  isotropicity)  is 
shown  in  Fig.  6(b).  The  linearity  of  the  probe  measurement 
(Oprobe j  mV/mm2)  with  respect  to  the  calorimeter  measure¬ 
ment  (^calorimeter’  mV/mm2)  was  derived  as  follows  by  fitting 
a  line  on  the  measured  data: 

,  _  ^probe  +  1-7296 

calorimeter  590  46  "  '  ' 

The  probe  voltage  signal  was  filtered  using  a  bandpass 
Butterworth  filter  to  remove  high-frequency  noise  in  the  probe 
response.  The  photon  fluence  is  proportional  to  integral  of  the 
probe  voltage  output  over  the  time  of  the  laser  pulse  and  was 
calculated  by  summing  over  the  measured  voltages  over  the 
pulse  duration  (Fig.  3).  This  is  a  novel  technique  used  in  this 
study  and  was  found  to  be  a  much  better  representation  of  the 
probe  output  as  compared  to  measuring  the  absolute  voltage. 

The  probe  voltage  and  calorimeter  power  were  normalized 
to  the  respective  surface  areas  [probe  diameter  =1.5  mm, 
probe  surface  area  (sphere)  =  7.065  mm2,  calorimeter 
aperture  =12  mm,  calorimeter  surface  area  (circle)  = 
113.04  mm2]. 

C.  Absolute  Photon  Quantification  in  Tissue-Like 
Media 

The  linear  relationship  in  Eq.  (1)  was  modified  to  derive  the 
optical  fluence  in  liquid  brain  phantoms.  The  light  entering 
the  fiber  within  the  probe  depends  on  the  amount  of  multiple 
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Fig.  6.  (a)  Probe  isotropicity  of  titanium  versus  nylon  probes.  The 

coefficient  of  variation  (stdev/mean  rms  voltage)  was  calculated 
for  probes  made  of  nylon  and  titanium.  The  error  bars  show  minimum 
and  maximum  variation  in  individual  probe  response.  The 
Ti02  -based  probes  showed  more  isotropicity  in  equatorial  and  azimu¬ 
thal  directions,  (b)  Probe  linearity:  graph  of  total  voltage  (area  under 
pulse)  (mV/mm2)  measured  by  a  1 .5  mm  diameter  Ti02  probe  versus 
the  power  measured  by  a  calorimeter  per  unit  surface  area 
(mW/mm2),  The  probe  response  is  linear  over  a  dynamic  range  of 
0.01074  mW/mm2  to  0.031316  mW/mm2  as  measured  by  the 
calorimeter.  The  error  bars  show  a  3%  variation  (systematic  error) 
in  pulsed  laser  power  fluctuation. 


photon  scattering  within  the  spherical  probe  tip,  which  in  turn 
depends  on  the  ratio  of  the  refractive  index  of  the  probe  with 
respect  to  the  medium  [19].  The  refractive  index  of  the 
titanium  probe  tip  should  lie  close  to  that  of  the  epoxy 
(n  =  1.53)  [39]  since  the  concentration  of  titanium 

(n  =  2.73)  [38]  is  low  (<20  mg/mL).  Since  the  refractive  in¬ 
dex  of  air  is  less  than  that  of  water,  the  probe  response  in  water 
is  smaller  as  more  light  leaves  the  probe  in  water  compared  to 
air.  To  account  for  the  difference  in  the  reflectance  between 
probe-air  and  probe-water  interfaces  a  calibration  factor  was 
measured.  The  probe  was  placed  inside  a  cuvette,  0.2  cm  from 
the  surface.  First  the  probe  response  voltage  was  measured  in  air 
(<^air)  and  then  measured  again  by  filling  the  cuvette  with  water 
water)’  without  changing  the  position  of  the  probe  with  re¬ 
spect  to  the  cuvette  and  the  laser  source.  By  accounting  for 
the  attenuation  by  water  (wwater  —  0.020  cm-1;  0.2  cm  path 
length)  and  the  difference  in  the  beam  spread  (change  in  refrac¬ 
tive  index  for  cuvette-to-air  versus  cuvette-to-water),  the 
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calibration  factor  (CF)  was  determined  experimentally  to  be 
1.32,  using  the  following  equation: 

FWHMair  3 

CF  =  Z^Z_  ^FWHM^y 

& water 

The  FWHM  values  were  measured  for  both  these  setups 
using  the  probe  (FWHMair  =  40  mm,  FWFiMwater  =  42  mm) 
and  their  ratio  (0.9091)  denotes  the  proportional  change  in 
fluence  due  to  beam  spread.  The  experimentally  derived  cali¬ 
bration  factor  is  then  used  for  photon  quantification  within 
tissue-like  media: 


^calorimeter 


{CF  *  <£probe)  +  1.7296 
390.46 


(3) 


where  ^calorimeter  ls  the  absolute  energy  fluence  rate  in 
mW/mm2  and  Oprobe  is  the  fluence  measured  by  the  probe 
(mW/mm2).  The  absolute  energy  fluence  rate  can  be  con¬ 
verted  to  energy  fluence  (dTbsolute)  (J/mm2)  by  integrating  it 
over  the  illumination  time. 


D.  Beam  Profile  Modeling  Using  the  Monte  Carlo 

The  original  code  of  the  GPU-based  3D  MC  was  advanced  to 
emit  a  broad-beam  super-Gaussian  beam  profile  emanating 

from  an  optical  fiberguide:  (/)Sg  —  e~°'5(^n,  where  r  =  z 
tan  6,  6  is  the  divergence  angle,  z  is  the  depth  along  the  beam 
axis,  a  is  the  standard  deviation  of  the  super-Gaussian  beam,  n 
is  the  order  of  the  super-Gaussian  distribution.  For  these  stud¬ 
ies,  the  measured  broad  super-Gaussian  beam  was  simulated 
using  multiple  point  sources,  each  with  a  super-Gaussian  out¬ 
put  (a  =  3.8  and  n  =  2.9)  for  a  source  (fiberguide)  of  radius 
0.3  cm.  In  addition  to  the  super-Gaussian  beam,  a  flat  broad 
beam  with  radius  equal  to  the  output  port  of  the  integrating 
sphere  (port  radius  0.9  cm)  was  implemented,  similar  to  other 
studies.  The  beam  characteristics,  phantom  dimensions,  and 
optical  properties  of  the  phantom  were  inputted  to  the  MC 
code  (and  empirical  model  code)  to  emulate  the  illumination 
conditions  used  in  the  optical  dosimetry  setup  of  Fig.  7(a). 
A  total  of  108  photons  were  simulated  per  beam. 

E.  Validation  of  the  3D  Monte  Carlo  in  Brain 
Phantoms  Using  Optical  Dosimetry 

As  described  in  Section  2.A,  liquid  optical  phantoms,  resem¬ 
bling  white  matter,  gray  matter,  and  skull  bone,  were  designed 
using  predetermined  concentrations  of  India  ink  and  intralipid. 
The  setup  for  the  light  dosimetry  in  the  phantoms  is  similar 
to  that  shown  in  Fig.  7(a).  Cuvette  dimensions  were 
5.3  cm  x  4.6  cm  x  3.1  cm,  thickness  =  0.4  cm.  Two  beam 
profiles — flat-top  and  super-Gaussian — were  used  to  illumi¬ 
nate  these  phantoms  with  the  laser/optical  parametric  oscillator 
(OPO)  system  tuned  to  800  nm.  The  3D  energy  fluence  rate 
distribution  was  determined  by  translating  the  optical  dosim¬ 
etry  probe  positioned  within  the  brain  phantom  using  a  three- 
axis  micrometer  stage.  These  measurements  were  compared  to 
Monte  Carlo  generated  data  to  test  the  validity  and  accuracy  of 
the  Monte  Carlo  code  (Figs.  8-11).  The  error%  was  averaged 
over  profiles  measured  and  calculated  at  different  depths  within 
the  phantom  along  and  perpendicular  to  the  beam  axis. 

In  order  to  quantify  the  photon  energy  fluence  (in  J/mm2) 
brain  tissues,  we  converted  the  output  of  the  Monte  Carlo 
fluence  (Opbotons/mm2)  as  follows: 


Optical  Dosimetry  Probe 


Plexiglas  Cuvette  filled  with 
Liquid  Phantom 


N| 


/] 


NIR  Photon  source: 

1.  Broad  beam  (Integrating 
Sphere) 

2.  Super-Gaussian  beam 
(Optical  Fiberguide) 


(b)  Homogeneous  phantoms  (liquid)  Heterogeneous  layered 


Gray  matter  White  matter 

phantom  phantom 


Gray 
Matter  skull 
phantom  phantom 


Fig.  7.  (a)  Setup  for  optical  dosimetry  in  a  liquid  phantom.  The 

optical  dosimetry  probe  is  translated  within  a  Plexiglas  cuvette  filled 
with  liquid  phantom,  and  illuminated  by  an  NIR  source.  This  setup  is 
used  for  intralipid  calibration  as  well  as  for  actual  phantom  studies, 
(b)  Pictures  of  homogeneous  liquid  phantoms  (left)  and  heterogeneous 
layered  agar  phantoms  (right).  In  the  heterogeneous  phantom  a  white 
matter  liquid  phantom  (not  shown)  is  added  on  top  of  the  gray  matter 
solid  agar  layer  for  dosimetry  measurements. 


Probe  versus  Monte  Carlo  (White  Matter,  Gray  Matter) 
(Source:  Integrating  Sphere) 


o.ooi 


z 

Distance  (cm) 


X  MC  White  Matter 
Phantom 


■  Probe  White 
Matter  Phantom 


—  MC  Gray  Matter 
Phantom 


A  Probe  Gray  Matter 
Phantom 


Fig.  8.  Optical  dosimetry  in  white  and  gray  matter  phantoms 
(cuvette  dimensions:  3.3  cm  x  4.6  cm  x  3.1  cm,  thickness  =  0.4  cm). 
The  photon  source  used  is  an  integrating  sphere  connected  to  a  pul¬ 
sating  laser  source.  The  graphs  show  the  exponential  decrease  of  pho¬ 
ton  fluence  along  the  beam  axis  in  white  and  gray  matter  mimicking 
phantoms.  The  Monte  Carlo  generated  fluence  closely  matches  the 
fluence  measured  by  the  probe. 


48 


Research  Article 


Vol.  55,  No.  34  /  December  1  2016  /  Applied  Optics  9881 


Probe  versus  Monte  Carlo  (super-Gaussian  beam) 
(source:  Optical  fiberguide) 


Fig.  9.  Optical  dosimetry  in  white  matter,  gray  matter,  and  skull 
bone  phantoms  (cuvette  dimensions:  5.3  cm  x  4.6  cm  x  3.1  cm, 
thickness  =  0.4  cm).  The  photon  source  used  is  an  optical  fiberguide 
(connected  to  a  pulsating  laser  source)  with  a  super-Gaussian  beam 
distribution.  The  Monte  Carlo  generated  fluence  closely  matches 
the  fluence  measured  by  the  probe.  This  shows  that  the  Monte 
Carlo  beam  modeling  can  be  reliably  used  as  a  reliable  estimator  of 
photon  energy  distribution  in  tissue  phantoms. 


0/mm2)  =  (pMC  *  (4) 

where  h  is  Planck’s  constant,  c  and  A  are  the  speed  and  wave¬ 
length  of  light  in  the  medium.  The  energy  deposited  (Edose  in 
J/mm3)  in  the  tissue  is  obtained  by  multiplying  the  photon 
dose  with  the  absorption  coefficient  fia: 

he 

Eiost  0/mm3)  =  (puc  *na*—.  (5) 

F.  Empirical  Algorithm  for  Photon  Propagation 

The  empirical  model  approach  is  based  on  the  assumption  that 
the  photon  fluence  in  a  voxel  in  a  particular  layer  is  a  weighted 
sum  of  the  fluences  of  neighboring  voxels  in  the  previous  layer. 
The  schematic  of  this  algorithm  is  shown  in  Fig.  12.  The  steps 
of  the  algorithm  can  be  summarized  as  follows: 

1.  Assign  the  incident  photon  beam  distribution  (e.g., 
super-Gaussian  beam)  to  the  first  layer  of  voxels  entering 
the  tissue  phantom.  In  this  case,  the  direction  of  photon  propa¬ 
gation  is  along  the  +z  direction. 

2.  For  every  voxel  in  a  layer,  calculate  the  14  directional 
scatter  components  and  accumulate  total  fluence  by  summing 
the  directional  fluences  in  each  voxel: 

If  {i,j,  k)  =  k  -  1)  +  Idr(/'  -  1  ,j,  k  -  l)wt2 

+  Idf (i,j  -\,k  -  l)wt2  +  Idb (i,j  l)wt2 

+  Idl(z  +  1  ,j,  k  -  l)wt2,  (6) 


I(i,j,  k )  =  If  (i,j,  k )  +  Idr (i,j,  k) 

+  Idl (i,j,  k)  +  Idf  (i,j,  k)  +  Idb  (i,j,  k) 

+  Isr (i,j,  k)  +  Isl (i,j,  k)  +  Isf  (i,j,  k) 

+  Isb  (i,j,  k)  +  Idbr (i,j,  k )  +  Idbl (i,j,  k) 

+  Idbf  (i,j,  k )  +  Idbb (i,j,  k )  +  lb (i,j,  k ).  (7) 

The  nomenclature  for  the  14  directional  fluence  vectors  is 
shown  in  Table  3.  Equation  (6)  represents  an  example  of  how 
the  forward  scatter  component  is  calculated.  The  coordinates  /, 
j,  k  denote  the  x,  y,  and  z  location  of  the  voxel.  The  weights 
wtl  and  wt2  are  forward  and  diagonal  scattering  weights. 
The  sum  of  all  weights  is  always  equal  to  one  (i.e., 
wtl  +  4  *  wt2  =  1),  so  that  the  total  fluence  entering  a  voxel 
is  equal  to  the  sum  of  fluence  absorbed  and  exiting  the  voxel. 
The  equations  used  to  calculate  the  other  components  can  be 
similarly  determined.  The  sum  of  all  the  scatter  components  in 
a  voxel  equals  the  total  fluence  in  that  voxel  (arbitrary  units: 
photons/unit  volume),  as  shown  by  Eq.  (7). 

3.  Calculate  absorbed  fluence  (labs)  using  the  Beer- 
Lambert  law  for  a  path  length  dl  equal  to  the  voxel  length, 
based  on  the  absorption  coefficient  (jta)  of  the  voxel: 

(8) 

4.  Calculate  and  accumulate  the  fluence  over  all  voxels  in 
the  current  layer  and  increment  z  to  proceed  to  the  next  layer. 
Continue  propagation  until  the  fluence  exits  the  volume. 

The  weights  to  be  assigned  for  scatter  were  obtained  by  an 
iterative  optimization  phase,  which  compares  the  output 
fluence  of  the  empirical  algorithm  to  that  generated  by  the 
Monte  Carlo  routine  for  voxels  along  the  beam  axis,  as  calcu¬ 
lated  by  the  R-square  coefficient  (. R 2).  Table  4  shows  the 
weights  (wtl  and  wt2)  obtained  for  a  pencil  beam  simulation 
within  different  tissue  types.  Figure  11  shows  an  example  of 
the  match  in  fluence  profiles  between  the  Monte  Carlo  and 
the  empirical  approach. 

3.  RESULTS  AND  DISCUSSION 

A.  Optical  Dosimetry  Probe  Calibration 

The  isotropicity  (coefficient  of  variation)  of  nylon-  and 
titanium-based  optical  dosimetry  probes  are  shown  in 
Fig.  5(a).  The  titanium-based  probes  demonstrated  a  signifi¬ 
cantly  improved  isotropicity  with  an  equatorial  coefficient  of 
variation  0.038-0.072  versus  0.193  for  the  nylon  and  an  azi¬ 
muthal  coefficient  of  variation  ranging  from  0.21—0.29  versus 
0.44  for  the  nylon.  This  is  a  2-5  x  overall  better  isotropicity, 
depending  on  the  concentration  of  Ti02.  Figure  6(b)  shows  the 
linearity  of  response  of  the  1.5  mm  titanium  probe.  The  mea¬ 
surement  error  in  probe  response  was  1.7%,  and  1.1%  in 
the  calorimeter  readings,  where  the  majority  of  this  error 
was  due  to  variation  in  the  output  power  of  the  laser  (3% 
absolute  deviations).  Thus,  the  titanium  probes  proved  to  be 
better  suited  to  optical  dosimetry  studies  compared  to  nylon 
probes  and  hence  were  used  in  brain  phantom  studies. 
While  nylon  probes  have  been  used  before,  this  is  the  first  time 
that  titanium  probes  were  designed  and  shown  to  have  a  better 
isotropic  response.  Previous  studies  have  measured  the  isotro¬ 
picity  of  dosimetry  probes  made  of  Arnite  (1 1%),  nylon  (1 1%), 
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Fig.  10.  Lateral  beam  profile  measurement  using  optical  dosimetry  probe  and  validation  by  Monte  Carlo  in  homogeneous  phantoms.  The  laser 
source  used  was  an  optical  fiberguide  (0.3  cm  diameter)  connected  to  a  pulsed  NIR  laser  emitting  at  a  wavelength  of  800  nm.  Cuvette  dimensions: 
3.3  cm  x  4.6  cm  x  3.1  cm,  thickness  =  0.4  cm.  (a)  Beam  characterization  in  air  using  calorimeter.  Optical  dosimetry  to  measure  and  validate  lateral 
beam  profiles  in  phantoms  resembling:  (b)  white  matter,  (c)  gray  matter,  and  (d)  skull  bone,  i.e.,  perpendicular  to  the  beam  axis.  The  measurements 
have  been  normalized  to  the  maximum  value  measured  by  the  dosimetry  probe  in  each  of  the  respective  phantoms.  The  Monte  Carlo  generated 
fluence  closely  matches  the  probe-measured  fluence  and  effectively  models  the  beam  spread  in  tissue  phantoms. 


Probe  versus  Monte  Carlo:  Lateral  Profile  (Heterogeneous  Phantom) 
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Fig.  11.  Validation  of  Monte  Carlo  in  a  heterogeneous  phantom  (1.1  cm  skull  agar,  1.1  cm  gray  matter  agar,  and  3.9  cm  white  matter  liquid). 
Lateral  profile  measurements  (in  the  white  matter  phantom  solution)  are  shown  on  the  left  while  axial  comparison  is  on  the  right.  The  depth  shown 
by  z  is  the  distance  from  the  gray  matter  agar  surface. 
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Fig.  12.  Empirical  method  for  3D  photon  propagation.  The  fluence  at  each  voxel  in  a  layer  is  a  weighted  sum  of  the  neighboring  fluences  in 
adjoining  voxels.  Part  A  shows  the  simplified  schematic  of  fluence  accumulation  (in  2D),  while  the  actual  3D  model  with  directional  fluence  com¬ 
ponents  is  shown  in  part  B  of  the  figure.  The  fluence  in  each  voxel  consists  of  forward,  backward,  side,  and  diagonal  fluence  components  (mentioned 
in  Table  3).  The  directional  fluences  are  first  calculated  in  each  voxel  along  with  the  total  fluence,  before  being  propagated  to  neighboring  voxels. 


Helioseal  dental  sealant  (20%)  [17,19].  Our  study  measured 
the  %  variation  of  the  best  titanium  probe  to  be  13%  and 
the  best  nylon  probe  to  be  17%.  The  differences  in  probe  con¬ 
struction  and  calibration  setup  contribute  to  different  measured 
isotropicity  of  the  same  probe  material  (e.g.,  nylon)  across  vari¬ 
ous  studies. 

To  quantify  the  photon  fluence  in  tissues  and  tissue-like  me¬ 
dia,  the  probe  voltage  measured  within  the  medium  must  be 
multiplied  by  a  set  of  calibration  factors  to  account  for  the 
change  in  photon  acceptance  and  loss  in  the  probe-medium 
interface  with  respect  to  the  probe-air  interface.  Since  the 


refractive  index  of  tissues  is  nearly  the  same  as  that  of  water, 
the  calibration  factor  obtained  for  water  can  be  used  for 
biological  media.  The  probe  calibration  factor  was  determined 
to  be  1.32.  The  technique  presented  in  Section  2.C  [using 
Eqs.  (1) — (3)]  provides  a  simple  and  new  method  to  determine 
the  calibration  factor  experimentally. 

B.  Validation  of  Monte  Carlo  in  Homogeneous 
Phantom  Using  the  Optical  Dosimetry  Probe 

The  optical  fluence  measured  by  the  ODPs  in  white  and  gray 
matter  phantoms  was  compared  to  the  Monte  Carlo  generated 
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Table  3.  Vector  Directions  and  Nomenclature  Empirical 
Algorithm 


Vector 

Name 

Vector  Component 

5  (x,  y,  z) 

Description 

If 

(0  0  1) 

Upward  vector  in  -\-z  direction 

lb 

(0  0  -1) 

Downward  vector  in  -£  direction 

Isf 

(0  1  0) 

Side  vector  in  +y  direction 

Isb 

(0  -1  0) 

Side  vector  in  y  direction 

Isr 

(1  0  0) 

Side  vector  in  -\-x  direction 

Isl 

(-1  0  0) 

Side  vector  in  -.%•  direction 

Idf 

(0  0.5  0.5) 

Upward  diagonal  vector  in 
-\-z  and  +y  direction 

Idb 

(0  -0.5  0.5) 

Upward  diagonal  vector  in 
+z  and  -y  direction 

Idr 

(0.5  0  0.5) 

Upward  diagonal  vector  in 
-\-z  and  -\-x  direction 

Idl 

(-0.5  0  0.5) 

Upward  diagonal  vector  in 
-\-z  and  -x  direction 

Idbf 

(0  0.5  -0.5) 

Downward  diagonal  vector 
in  -£  and  +y  direction 

Idbb 

(0  -0.5  -0.5) 

Downward  diagonal  vector 
in  -£  and  -y  direction 

Idbr 

(0.5  0  -0.5) 

Downward  diagonal  vector 
in  -£  and  +x  direction 

Idbl 

(-0.5  0  -0.5) 

Downward  diagonal  vector 
in  -£  and  -at  direction 

Table  4.  Optimization  of  Empirical  Algorithm  for  Brain 
Tissue  Phantoms 


Tissue 

Phantom 

Resolution 

(mm) 

Wtl 

Wt2 

R2 

Coefficient 

White  matter 

0.1 

0.949461 

0.012635 

0.99999 

Gray  matter 

0.1 

0.87469 

0.031328 

0.999969 

Skull  bone 

0.1 

0.914828 

0.021293 

0.99979 

White  matter 

1 

0.4748 76 

0.131281 

0.99992 

Gray  matter 

1 

0.062097 

0.234476 

0.999828 

Skull  bone 

1 

0.245949 

0.188513 

0.999871 

fluence  maps  to  validate  the  Monte  Carlo.  These  experiments 
necessitate  an  accurate  determination  of  the  phantom’s  optical 
properties.  Figures  8  and  9  show  the  comparison  between  the 
probe  measurements  and  the  Monte  Carlo  in  brain  phantoms 
for  two  different  illumination  profiles,  flat  diffuse  beam  (inte¬ 
grating  sphere)  and  super-Gaussian  beam  (optical  fiberguide). 
A  comparison  was  also  made  in  the  direction  perpendicular  to 
the  beam  axis,  as  shown  in  Fig.  10.  From  these  data,  the  average 
relative  %  error  of  the  Monte  Carlo  simulation  relative  to  the 
probe  measurements  was  11.2%  for  white  matter,  13.2%  for 
gray  matter,  and  1 1.8%  for  skull  bone  along  the  beam  axis  for  a 
range  of  0.2-2. 0  cm  within  the  phantom,  while  perpendicular 
to  the  beam  axis,  the  %  relative  error  in  the  lateral  profile  was 
9.4%,  12.0%,  and  8.9%  for  white  matter,  gray  matter,  and 
skull  bone  phantoms,  respectively.  These  systematic  errors  re¬ 
present  the  level  of  confidence  with  which  the  Monte  Carlo 
matches  the  probe  measurements  and  validate  their  use  in  sim¬ 
ulating  optical  energy  distribution  in  media  (Tables  3  and  6). 

An  analysis  of  the  fluence  profile  within  the  skull  phantom 
showed  that  the  super-Gaussian  beam  fluence  decreases  with  an 


Table  5.  %  Mismatch  Error  in  Homogeneous  Phantom 
Along  Beam  Axis  (z) 


Phantom 

Error% 
(0.2-1  cm) 

Error% 

(1-2  cm) 

Average 

Error% 

White  matter 

12.13 

10.5 

11.27 

Gray  matter 

13.39 

30.12 

13.25 

Skull  bone 

11.39 

12.28 

11.81 

effective  attenuation  coefficient  of  4.81  cm1.  This  factor  was 
used  to  estimate  the  effectiveness  of  therapy  for  different  skull 
thicknesses.  Based  on  Monte  Carlo  simulated  data,  the  powers 
transmitted  by  the  fiberguide  source  for  various  skull  thick¬ 
nesses  are:  5.9  mW/cm2  (0.4  cm  thick),  2.16  mW/cm2 
(0.6  cm  thick),  0.83  mW/cm2  (0.8  cm  thick),  0.33  mW/ 
cm2  (1  cm  thick).  Thus,  the  %  power  transmitted  in  a 
1  cm  thick  adult  skull  is  1.8%,  while  in  a  pediatric  skull  with 
a  thickness  of  4.0  mm,  33.5%  optical  power  would  reach  the 
soft  tissue,  representing  18.5  times  more  fluence. 

C.  Validation  of  Monte  Carlo  in  Heterogeneous 
Phantom  Using  the  Optical  Dosimetry  Probe 

Energy  fluence  distribution  was  measured  in  the  white  matter 
volume  of  a  heterogeneous  head  phantom,  which  consisted  of  a 
1 . 1  cm  thick  skull  agar  phantom,  a  1 . 1  cm  thick  gray  matter 
agar  phantom,  and  a  5.9  cm  thick  white  matter  liquid  phantom 
(Fig.  7).  The  exponential  decreases  in  energy  fluence  along  the 
beam  axis  and  the  lateral  beam  profiles  at  different  depths  above 
the  gray  matter  agar  surface  are  shown  in  the  right  and  left  plots 
of  Fig.  11.  The  lateral  fluence  profiles  have  been  normalized  to 
the  intensity  at  a  depth  of  0.2  cm  above  the  gray  matter  agar 
surface,  while  the  fluence  profiles  along  the  beam  axis  were  nor¬ 
malized  to  the  fluence  at  z  =  0.1  cm  above  the  gray  matter 
agar  surface.  The  average  error  (relative  %  error)  between 
the  probe  and  Monte  Carlo  fluence  was  17.7%  along  the  beam 
axis  and  18.0%  perpendicular  to  the  beam  axis  (Tables  7  and 
8).  These  results  demonstrate  that  the  beam  spread  can  be  de¬ 
termined  accurately  by  the  Monte  Carlo  to  an  overall  depth  of 
nearly  4  cm  given  the  optical  coefficients  of  the  head  phantom. 

Past  studies  have  compared  the  Monte  Carlo  simulations  to 
dosimetry  probe  measurements  for  simple  beam  profiles,  such 
as  flat  beams  or  isotropic  sources  [14,15].  One  study  compared 
the  angular  radiance  measured  by  a  flat  cleaved  dosimetry  probe 
to  that  simulated  by  the  Monte  Carlo  and  determined  an  ac¬ 
curacy  of  around  25%  (see  Fig.  6  in  referenced  publication) 
[14].  Another  study  validated  the  Monte  Carlo  in  pig  bronchus 
tissue  using  dosimetry  probe  measurements  with  errors  up  to 
50%  at  4  cm  depth  and  systematic  variation  of  15%  [15].  In 
past  studies,  the  long  computation  times  (>24  h)  associated 
with  the  CPU-based  Monte  Carlo  prevented  a  rigorous  ap¬ 
proach  to  Monte  Carlo  validation  across  diverse  tissue  types 
and  under  different  illuminations.  Another  source  of  error  in 
these  studies  was  the  uncertainty  in  the  estimation  of  tissue/ 
phantom  optical  properties.  In  this  study,  we  have  minimized 
this  uncertainty  by  independently  validating  the  optical  proper¬ 
ties.  The  use  of  a  GPU-based  Monte  Carlo  reduces  the  develop¬ 
ment  time  (few  hours)  associated  with  optimizing  the  Monte 
Carlo  for  complex  illumination  and  boundary  conditions. 
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Table  6.  %  Mismatch  Error  in  Homogeneous  Phantom  Perpendicular  to  Beam  Axis 


Phantom 

Error%  z  =  0.2  cm 

Error%  z  =  0.6  cm 

Error%  z  =  1.0  cm 

Error%  z  =  1.4  cm 

Average  Error% 

White  matter 

10.34 

13.34 

9.14 

4.79 

9.4 

Phantom 

Error%  z  =  0.2  cm 

Error%  z  =  0.4  cm 

Error%  z  =  0.6  cm 

Error%  z  =  1  cm 

Average  error% 

Gray  matter 

13.45 

11.52 

14.29 

8.97 

12.06 

Phantom 

Error%  z  =  0.6  cm 

Error%  z  =  0.7  cm 

Error%  z  =  0.8  cm 

Average  error% 

Skull  bone 

8.41 

5.69 

12.64 

8.91 

Table  7.  %  Mismatch  Error  in  Heterogeneous  Phantom 
Along  Beam  Axis  (z) 


Phantom 

Error°/o 
(0.1-1  cm) 

Error% 
(1-2  cm) 

Average 

Error% 

White  matter  (over  skull  and  2.11 

gray  matter  layers) 

33.32 

17.71 

Table  8.  %  Mismatch  Error  in  Heterogeneous  Phantom 
Perpendicular  to  Beam  Axis 

Phantom 

Error%  Error%  Error% 

z  =  z  =  z  = 

0.2  cm  0.4  cm  0.6  cm 

Error% 

z  = 

1  cm 

Average 

Error% 

White  matter 
(over  skull 
and  gray 
matter  layers) 

18.52  15.09  15.37 

23.17 

18.04 

The  use  of  pulsed  laser  beams  has  been  shown  to  improve 
the  instantaneous  temperature  rise  and  release  of  drug 
molecules  tethered  to  nanoparticles,  compared  to  continu¬ 
ous-wave  lasers  [42] .  This  will  allow  more  instantaneous  power 
to  be  coupled  to  the  nanomolecules  while  reducing  the  chances 
of  laser-induced  adverse  tissue  heating.  Assuming  nanosecond 
pulse  duration,  the  maximum  permissible  exposure  limit  of  la¬ 
ser  for  skin  is  3.17  x  107  W  cm-2  [41,42].  Based  on  our  MC 
simulations  with  a  flat  2  cm  diameter  beam,  the  amount  of 


fluence  coupled  to  the  soft  tissue  for  a  1.1  cm  thick  skull  is 
4.3  x  105  W  cm-2  and  1.5  x  107  W  cm-2  for  neonatal  heads 
with  a  0.2  cm  thick  skull.  Thus,  a  drug  needs  to  be  designed 
to  activate  at  a  fluence  threshold,  less  than  the  maximum  flu¬ 
ence  coupled  to  the  tissues  at  a  specific  depth  where  the  meta¬ 
stasis  is  located.  For  example,  a  drug  molecule  with  an 
activation  threshold  of  1  W  cm-2  can  be  activated  up  to 
3.8  cm  in  the  white  matter  of  a  20  cm  diameter  adult  brain 
(1.1  cm  thick  skull  plus  2  cm  gray  matter)  and  up  to 
4.6  cm  in  a  8  cm  diameter  neonatal  brain  (0.2  cm  skull  plus 
0.2  cm  gray  matter).  These  results  demonstrate  that  NIR  acti¬ 
vated  drug  release  is  more  feasible  in  neonatal  brains  where 
4.6  cm  is  greater  than  half  of  the  brain  diameter.  The  temporal 
widening  of  the  laser  pulse  with  depth  in  the  brain  is  also  seen 
to  affect  the  instantaneous  fluence  (W  •  cm-2  •  s)  coupled  to 
tissues.  Based  on  our  Monte  Carlo  simulations,  the  pulse  width 
is  seen  to  be  4-5  times  at  depths  of  3-5  cm  for  a  1  ns  wide 
input  laser  pulse.  This  implies  a  decrease  in  instantaneous  flu¬ 
ence  by  that  amount.  The  Monte  Carlo  can  be  used  to  predict 
the  temporal  beam  widening  in  order  to  determine  the  instan¬ 
taneous  fluence  at  specific  depths. 

D.  Empirical  Versus  Monte  Carlo 

In  Figs.  12  and  13,  a  schematic  of  the  EM  used  in  these  studies 
is  presented  and  the  corresponding  weights  determined  by 
matching  photon  fluence  results  to  MC  data.  As  seen  in 
Table  4,  the  weights  obtained  during  optimization  depend 
on  the  voxel  size  and  the  scattering  coefficient  of  the  medium. 
The  scattering  coefficient  and  the  voxel  size  determine  the 
number  of  scattering  events  per  unit  volume,  which  affects 


♦  MC  White  Matter 
Phantom 

■  MC  Gray  Matter 
Phantom 

A  MC  Skull  Bone 
Phantom 

X  Empirical  White 
Matter  Phantom 

x  Empirical  Gray 
Matter  Phantom 

+  Empirical  Skull 
Bone  Phantom 


♦  Probe  White 
Matter  Phantom 

O  MC  White 
Matter  Phantom 
a  Empirical  White 
Matter  Phantom 
A  MC  Gray  Matter 
Phantom 
x  Empirical  Gray 
Matter  Phantom 

•  Probe  Gray 
Matter  Phantom 

+  Probe  Skull  Bone 
Phantom 
-  MC  Skull  Bone 
Phantom 
■  Empirical  Skull 
Bone  Phantom 


Fig.  13.  (a)  Optimization  phase:  a  flat  broad  beam  (2  cm  diameter)  is  simulated  using  Monte  Carlo  and  the  empirical  algorithm  is  optimized  to 

match  the  Monte  Carlo  with  the  identical  illumination  conditions,  (b)  Comparison  of  the  empirical  algorithm  with  Monte  Carlo  and  probe- 
measured  fluence  in  brain  phantoms  resembling  white  matter,  gray  matter,  and  skull  bone.  The  illumination  source  is  a  super-Gaussian  beam. 
The  empirical  approach  is  a  close  approximation  of  the  Monte  Carlo  and  the  probe  over  the  range  of  optical  properties  of  brain  tissues. 
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the  percentage  scatter  in  the  forward  versus  diagonal  directions, 
thus  impacting  the  weights  (wtl  and  wt2).  Based  on 
these  weights,  EM  simulated  data  for  homogeneous  and 
heterogeneous  phantoms  were  compared  to  Monte  Carlo 
and  ODP  measurements.  This  comparison  was  performed 
along  the  beam  axis,  see  Fig.  13(b),  in  homogeneous  brain  tis¬ 
sue  phantoms.  The  average  relative  %  error  between  the  em¬ 
pirical  and  probe  measurements  for  the  super-Gaussian  beam 
simulation  is  10.1%  for  white  matter,  43.2%  for  gray  matter, 
and  22.1%  for  skull  bone  phantoms;  while  the  comparison  be¬ 
tween  empirical  and  MC  is  22.6%  for  white  matter,  35.8%  for 
gray  matter,  and  21.9%  for  skull  bone  measurements.  The  er¬ 
rors  between  probe  measurements  and  the  empirical  algorithm 
are  cumulative  in  nature,  and  can  be  reduced  with  better 
optimization  methods.  No  systemic  trend  was  observed  in  these 
errors  with  %  differences  having  both  positive  and  negative 
values.  Figure  14  shows  a  comparison  between  the  lateral 


fluence  profiles  generated  by  a  Monte  Carlo  and  the  empirical 
models  for  a  flat  beam  of  2  cm  diameter  in  homogenous  white 
and  gray  matter  tissues.  A  significant  match  is  seen  between  the 
normalized  fluence  profiles  at  two  different  depths.  Figure  15 
shows  a  complete  head  simulation  and  the  comparison  between 
the  MC  and  empirical  fluence  versus  depth  along  the  beam 
axis.  The  empirical  fluence  closely  follows  the  Monte  Carlo 
but  has  significant  errors  at  object  boundaries. 

The  Monte  Carlo  code  is  computationally  expensive  in  spite 
of  being  based  on  a  GPU-based  platform.  In  order  to  simulate 
the  optical  properties  and  fine  structural  boundaries  of  brain 
tissues,  a  resolution  equal  to  or  less  than  1  mm  is  necessary. 
Even  though  an  optical  map  of  the  brain  can  be  obtained 
through  segmenting  and  translating  CT  and/or  MR  images, 
the  dataset  size  of  the  brain  becomes  very  large,  e.g.,  a  CT  scan 
requires  a512x512x512  matrix.  This  significantly  increases 
the  execution  time  of  the  Monte  Carlo  routine  to  iterate  over 


Fig.  14.  Lateral  beam  profile  comparison  between  normalized  lateral  fluence  profiles  of  empirical  (red  line)  and  Monte  Carlo  (blue  line)  for  gray 
matter  and  white  matter  tissue  at  depths  of  5  and  10  mm  in  a  homogeneous  phantom.  A  flat  broad  beam  (20  mm  diameter)  is  simulated.  A  close 
match  is  seen  between  the  Monte  Carlo  and  empirical  lateral  beam  profiles. 
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Fig.  15.  (a)  Simulation  of  a  spherical  head  phantom  with  brain  tissue  properties  shown  in  Table  1.  Skull  thickness  =  4  mm;  gray  matter 

thickness  =10  mm,  white  matter  diameter  =  72  mm.  Uniform  illumination  beam  diameter  =  20  mm  (800  nm  wavelength), 
(b)  Comparison  of  fluence  along  the  beam  axis  between  empirical  (red  line)  and  MC  (blue  line)  models.  The  empirical  model  achieves  this  sim¬ 
ulation  in  6  s  while  the  MC  takes  8  h  for  a  total  of  1010  photons. 


various  illumination  conditions.  The  empirical  approach 
provides  an  alternative  to  the  computationally  expensive 
Monte  Carlo,  and  can  be  used  to  select  the  best  illumination 
conditions  for  therapy.  The  speed-up  achieved  by  the  empirical 
algorithm  was  approximately  700  x  for  a  pencil  beam  simula¬ 
tion  and  1 6,  000  x  for  a  broad-beam  simulation  in  a  head 
phantom.  Thus  while  a  broad-beam  head  simulation  takes 
6  s  with  the  empirical,  it  can  take  in  excess  of  8  h  with  the 
GPU-based  Monte  Carlo  for  a  head  phantom  simulation  of 
volume  1000  cm3.  It  is  important  to  note  that  the  number 
of  photons  to  be  simulated  depends  on  the  depth  at  which 
the  photon  distribution  is  to  be  evaluated,  and  the  %  fluence 
penetration  (and  systematic  errors  between  MC,  probe,  and 
empirical)  is  independent  of  the  number  of  photons  simulated. 

4.  CONCLUSION 

We  validated  the  3D  Monte  Carlo  using  a  novel  optical  dosim¬ 
etry  probe  in  phantoms  resembling  white  matter,  gray  matter, 
and  skull  bone.  The  Ti02 -based  dosimetry  probe  was  shown  to 
have  superior  linearity  and  isotropicity  of  response,  and  was 
better  suited  to  validate  the  Monte  Carlo  using  localized  3D 
measurement  (<23%  error)  in  tissue  phantoms.  The  indepen¬ 
dent  characterization  of  the  probes  and  the  optical  phantoms 
using  the  calorimeter  “gold  standard”  and  spectrophotometer 
(respectively)  further  validate  our  results.  The  3D  validation 
of  the  Monte  Carlo  generated  fluence  for  two  different  types 
of  NIR  light  sources  in  different  optical  phantoms,  and  the  ease 
of  photon  quantification  to  absolute  dose  (in  mW/mm2  or 
J / mm2)  demonstrates  the  fidelity  of  the  Monte  Carlo  as  a  pre¬ 
dictive  tool  to  accurately  estimate  3D  photon  energy  distribu¬ 
tion  in  complex  media.  The  measurement  and  prediction 
accuracy  can  be  improved  with  better  instrumentation  for 
probe  localization  and  better  laser  beam  positioning  setups. 
The  voxel-based  Monte  Carlo  can  emulate  the  complex 


geometrical  shape  of  tissues  using  imaging  datasets  such  as 
CT/MRI  and  can  now  be  used  as  a  predictive  clinical  tool 
to  estimate  photon  dose  in  heterogeneous  brain  tissues. 

The  empirical  algorithm,  validated  by  the  3D  Monte  Carlo, 
increases  the  clinical  feasibility  of  optical  therapeutic  planning  to 
narrow  down  the  complex  possibilities  of  illumination  conditions, 
further  compounded  by  the  heterogeneous  structure  of  the  brain 
(e.g.,  varying  skull  thicknesses  and  densities).  A  limitation  of  us¬ 
ing  the  empirical  approach  is  that  the  weights  assigned  depend  on 
the  resolution  of  the  medium  and  changing  resolution  requires  an 
optimization  phase  to  derive  the  weights.  A  rigorous  validation  of 
the  empirical  approach  is  currently  under  development,  which  is 
necessary  for  its  applicability  in  heterogeneous  tissue  environ¬ 
ments  in  clinical  settings.  Despite  its  limitations,  it  can  be  used 
as  a  preliminary  step,  with  the  final  estimation  done  by  the  Monte 
Carlo  to  ensure  high  predictive  accuracy. 

Our  ultimate  goal  is  to  design  a  fast  Monte  Carlo-based  op¬ 
tical  therapeutic  protocol  to  treat  brain  metastasis.  This  requires 
voxel  by  voxel  determination  of  photon  energy  distribution  to 
accurately  predict  the  rate/quantity  of  drug  release.  This  study 
achieves  the  first  step  in  the  design  of  a  tool  to  validate  the  NIR 
photon  distribution  in  three  dimensions  in  optical  phantoms, 
where  tissue  structure  and  boundaries  can  be  determined  from 
imaging  modalities  such  as  MRI  or  CT. 
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ABSTRACT:  Nanoparticle-based  platforms  for  gene  therapy  and 
drug  delivery  are  gaining  popularity  for  cancer  treatment.  To 
improve  therapeutic  selectivity,  one  important  strategy  is  to 
remotely  trigger  the  release  of  a  therapeutic  cargo  from  a  specially 
designed  gene-  or  drug-laden  near-infrared  (NIR)  absorbing  gold 
nanoparticle  complex  with  NIR  light.  While  there  have  been 
multiple  demonstrations  of  NIR  nanoparticle-based  release 
platforms,  our  understanding  of  how  light-triggered  release 
works  in  such  complexes  is  still  limited.  Here,  we  investigate  the 
specific  mechanisms  of  DNA  release  from  plasmonic  nanoparticle 
complexes  using  continuous  wave  (CW)  and  femtosecond  pulsed 
lasers.  We  find  that  the  characteristics  of  nanoparticle-based  DNA 
release  vary  profoundly  from  the  same  nanoparticle  complex, 
depending  on  the  type  of  laser  excitation.  CW  laser  illumination 
drives  the  photothermal  release  of  dehybridized  single-stranded  DNA,  while  pulsed-laser  excitation  results  in  double- 
stranded  DNA  release  by  cleavage  of  the  Au— S  bond,  with  negligible  local  heating.  This  dramatic  difference  in  DNA  release 
from  the  same  DNA— nanoparticle  complex  has  very  important  implications  in  the  development  of  NIR-triggered  gene  or 
drug  delivery  nanocomplexes. 

KEYWORDS:  nanoshells ,  DNA ,  laser,  photothermal  heating,  oligonucleotide 


To  treat  and  prevent  disease,  light-triggered  release  using 
nanoparticle-based  complexes  has  been  shown  to  be  a 
promising  strategy  for  gene  therapy1-15  and  drug 
delivery.16-25  One  popular  complex  consists  of  thiolated  DNA 
bound  to  near-infrared  (NIR)  absorbing  gold  nanoparticles 
such  as  nanoshells,  nanorods,  or  nanocages.12-15  Resonant 
illumination  of  the  gold  nanoparticle— DNA  complex  results  in 
the  light-triggered  release  of  DNA  from  the  surface  of  the 
nanoparticle.7-  For  gene  therapy,  released  DNA,  or  in  some 
cases  RNA,  can  be  used  to  prevent  transcription  or  translation 
of  messenger  RNA  (mRNA)  or  to  inhibit  protein  function.7,26 
This  type  of  nanoparticle— DNA  complex  can  also  be  used  for 
drug  delivery,  where  therapeutic  molecules  can  be  non- 
covalently  bound  to  the  DNA  and  released  upon  resonant 
illumination.23-25  While  other  nanoparticle  platforms  have 
been  investigated  for  drug  delivery,  including  liposomes20,27 
and  hydrogels,28,29  as  well  as  polymeric/0-32  magnetic,33,34  and 
metallic  nanoparticles,25  these  nanosystems  typically  rely  upon 
changes  in  the  intracellular  environment  such  as  pH32,35-37  and 
temperature29,31'  -  7  to  induce  drug  release,  which  can  vary 
with  cell  type  and  location.  Externally  applied  triggers  such  as 


magnetic38,39  and  electric40,41  fields,  ultrasound,42  and, 
specifically,  light, 23,40-46  allow  for  precise  spatial  and  temporal 
release.  Controlled  release  of  chemotherapeutic  agents  once  the 
particles  have  been  delivered  specifically  to  a  tumor  site  is  a 
potential  strategy  to  reduce  the  systemic  toxicity  of  traditional 
chemotherapy  treatments. 

The  experimental  conditions  employed  for  light-triggered 
release  schemes  vary  widely  from  study  to  study.8-10,12,13,25'47 
Many  schemes  illuminate  metallic  nanostructures  coated  in 
double-stranded  DNA  (dsDNA)  with  continuous  wave  (CW) 
lasers  to  dehybridize  the  dsDNA,  resulting  in  release  of  single- 
stranded  DNA  (ssDNA),  while  the  complementary  strand 
remains  tethered  to  the  nanoparticle  surface  through  a  gold- 
thiol  bond.7,8,23,48,49  Other  schemes  utilize  pulsed  (femto¬ 
second,  nanosecond,  and  picosecond)  lasers  to  either 
dehybridize  the  dsDNA  or  break  the  gold— thiol  bond  between 
the  DNA  molecule  and  the  nanoparticle,  which  results  in 
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release  of  the  entire  dsDNA.9'10'13,50-54  Two  mechanisms  have 
been  proposed  to  explain  these  release  observations:  either 
simple  melting  of  the  dsDNA  from  the  increased  temperature 
of  the  nanoparticles  or  hot-electron  transfer  from  plasmon 
decay  causing  either  dehybridization  or  Au— S  bond  breaking. 
To  improve  control  and  tailorability  of  light-triggered  DNA 
release  in  these  systems,  it  is  critical  to  understand  the  specific 
DNA  release  mechanisms  and  how  they  relate  to  the  type  of 
laser  illumination  used  to  trigger  release. 

In  this  study,  we  examine  the  specific  mechanisms  of  NIR 
light- triggered  DNA  release  from  DNA-functionalized  plas- 
monic  nanoparticle  complexes.  Here,  we  use  nanoshells  (NS), 
composed  of  a  silica  core  surrounded  by  a  thin  gold  shell, 
functionalized  with  partially  thiolated  dsDNA.  NS  are  strong 
absorbers  of  NIR  light,  are  biocompatible,  and  are  easily 
functionalized.  -  We  employ  NIR  light  because  it  enables 
deep  tissue  penetration  due  to  the  low  absorption  of  biological 
tissue  in  this  region,60  while  high-energy  visible  or  UV  photons 
are  strongly  absorbed  by  bulk  tissue,  which  can  result  in  protein 
and  DNA  damage.  1  We  find  that  DNA  release  can  be  achieved 
in  two  distinct  regimes:  high-power  CW  illumination  and  low 
average  power,  high-intensity  pulsed-laser  excitation. 

Each  regime  is  controlled  by  a  distinct  mechanism  that  drives 
the  DNA  release  (Scheme  l).  We  have  observed  that  DNA 


Scheme  1.  NIR-Light-Induced  DNA  Releasea 


aCW  irradiation  results  in  dehybridization  and  release  of  fluorescently 
tagged  ssDNA,  while  pulsed  irradiation  results  in  Au— S  bond  breakage 
and  release  of  dsDNA. 


release  upon  CW  illumination  is  due  to  a  temperature  increase 
arising  from  the  collective  heating  effects  of  the  nanoparticles, 
causing  dehybridization  and  release  of  ssDNA.  This  limits  the 
use  of  CW  illumination  in  vivo  because  the  high  temperatures 
are  likely  to  result  in  nonspecific  cell  death  and  also  because  it  is 
highly  concentration- dependent.  Our  findings  clarify  previous 
work,  which  postulated  that  CW  illumination  could  drive 


release  through  hot-electron  injection.  Alternatively,  we  have 
also  investigated  femtosecond  pulsed-laser-induced  release  at 
low  (25  mW  range)  average  laser  powers.  DNA  release  is  found 
to  occur  through  a  hot-electron  transfer  process,  which  breaks 
the  Au— S  bond,  releasing  dsDNA.  While  this  mechanism  has 
been  demonstrated  with  high-energy  UV  pulsed  lasers,51  to  our 
knowledge,  no  previous  work  has  demonstrated  the  pulsed 
release  mechanism  with  NIR  light.  We  carefully  monitored  the 
bulk  temperature  under  these  conditions  and  found  essentially 
no  bulk  temperature  increase  under  illumination.  We  also 
found  that,  for  our  experimental  conditions,  the  process  is  also 
nearly  independent  of  nanoparticle  concentration.  These  results 
indicate  that  pulsed-laser-induced  release  is  substantially  more 
favorable  for  in  vivo  applications  than  CW  release,  where  it  is 
highly  desirable  to  minimize  nonspecific  cell  death  and  where 
local  nanoparticle  concentrations  and  distributions  are 
uncontrollable  and  essentially  unknown. 

RESULTS  AND  DISCUSSION 

Nanocomplex  Synthesis  and  Characterization.  dsDNA 
was  formed  by  hybridizing  a  thiolated  ssDNA  with  a 
nonthiolated  fluorescein-tagged  ssDNA.  The  melting  temper¬ 
ature  of  the  dsDNA  in  solution  was  characterized  by  circular 
dichroism  (CD).  Spectra  show  a  negative  peak  at  246  nm  and  a 
positive  peak  at  278  nm,  characteristic  of  B-form  DNA  (Figure 
S1A).62'  3  The  ellipticity  (6)  at  246  nm  was  monitored  as  a 
function  of  temperature  to  obtain  the  dsDNA  melting  curve 
(Figure  SIB).  The  first  derivative  of  a  Boltzmann  fit  of  the  CD 
curve  showed  a  melting  temperature  (Tm)  of  55.2  ±  0.7  °C, 
which  is  close  to  the  theoretically  calculated  Tm  of  51.3  °C  in 
33  mM  NaCl  TE  buffer.64 

Nanoshells  were  synthesized  with  a  120  nm  Si02  core  and  a 
21  nm  Au  shell  according  to  previously  reported  protocols.55'57 
The  NS(o)dsDNA  nanocomplex  was  formed  by  attachment  of 
the  dsDNA  to  the  NS  surface  via  a  gold— thiol  bond. 
Functionalization  of  the  NS  surface  was  verified  with  extinction 
spectroscopy,  transmission  electron  microscopy  (TEM),  and  f- 
potential  measurements.  The  bare  NS  extinction  spectrum  red- 
shifted  from  775  to  781  nm  after  dsDNA  functionalization, 
resulting  from  the  change  in  the  dielectric  environment  around 
the  NS  (Figure  1A).  TEM  images  show  the  formation  of  a  thin 
layer  of  DNA  around  the  NS  (Figure  1B,C).  The  f -potential 
decreased  from  —44.1  ±  0.7  to  —59.4  ±1.7  mV  due  to  the 
increased  negative  charges  associated  with  the  phosphate 
backbone  of  the  DNA  (Figure  ID). 

Thermal  and  CW  Laser-Induced  DNA  Melting.  To 
investigate  the  dehybridization  temperature  when  dsDNA  is 
tethered  to  the  NS  surface,  a  thermal  release  profile  at  5  °C 


Figure  1.  Nanocomplex  characterization.  (A)  Extinction  spectra  of  bare  NS  (green)  and  NS(5)dsDNA  (blue),  TEM  images  of  (B)  bare  NS  and 
(C)  NS(S)dsDNA,  and  (D)  f -potential  of  bare  NS  (green)  and  NS(S)dsDNA  (blue).  Scale  bars  are  50  nm. 
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Figure  2.  Investigation  of  CW-induced  DNA  release  mechanism.  (A)  Heating  profile  of  DNA  release  from  solutions  heated  via  a  CW  laser  at 
various  laser  powers.  Inset:  Thermal  image  of  DNA  release  experiment  at  1.5  W.  DNA  release  profile  plotted  as  a  function  of  (B)  time  and  (C) 
temperature.  A  Boltzmann  fit  estimates  a  melting  temperature  of  59.1  ±  4.5  °C.  (D)  DNA  release  from  solutions  in  a  cuvette  stage  set  to  4  °C 
under  1.6  (red)  and  3.4  W  (black)  laser  irradiation  showing  no  DNA  release.  Inset:  Corresponding  heating  profile  showing  that  the  solution 
temperature  does  not  exceed  room  temperature. 


increments  was  acquired.  The  nanoparticle  concentration  used 
was  1  X  109  NS/mL,  a  concentration  range  similar  to  the  local 
concentration  of  nanoparticles  in  tumors  (108— 109  particles/ 
mL).65'66  The  NS(S)dsDNA  solution  was  heated  in  a  temper¬ 
ature-controlled  stage.  After  20  min  of  equilibration  at  each 
chosen  temperature,  an  aliquot  was  removed  and  centrifuged.  A 
fresh  sample  was  used  for  each  temperature  point.  A  fluorescein 
tag  on  the  nonthiolated  strand  was  used  for  quantification  of 
the  released  ssDNA.  The  fluorescence  intensity  of  the 
supernatant  was  measured,  and  these  data  were  used  to  obtain 
the  amount  of  ssDNA  released  from  NS  using  a  linear  standard 
curve  (Figure  S2A).  The  ssDNA/NS  released  versus  temper¬ 
ature  was  found  to  follow  a  Boltzmann  response;  the  DNA 
melting  temperature  achieved  by  heating  the  samples  was  61.1 
±  2.3  °C  (Figure  S2B).  The  slightly  higher  Tm  from  thermal 
release  compared  well  to  the  Tm  of  55.2  ±  0.7  °C  observed  in 
CD,  where  the  difference  can  be  attributed  to  the  dsDNA  being 
anchored  to  the  NS  surface  versus  free  dsDNA  in  solution 
obtained  from  the  CD  measurements.  The  binding  of  the  DNA 
to  the  nanoparticle  results  in  increased  steric  hindrance  and  a 
difference  in  local  salt  concentration  relative  to  DNA  in 
solution;  these  differences  in  local  environment  are  sufficient  to 
alter  the  temperature  between  those  two  regimes,  consistent 
with  our  observations.67 

We  then  compared  the  thermal  release  of  DNA  from  the  NS 
to  NIR  release  induced  by  an  808  nm  wavelength  CW  laser. 
The  NS(S)dsDNA  solution  was  irradiated  at  0.4,  0.7,  0.8,  1.0, 
and  1.5  W  for  times  ranging  from  1.5  to  25  min.  During 
irradiation,  the  sample  temperature  was  monitored  with  a 
thermal  camera  (Figure  2A,  inset).  When  plotted  as  a  function 
of  time,  no  substantial  release  was  observed  except  at  laser 
power  levels  of  1.5  W,  which  resulted  in  the  bulk  solution 
temperature  reaching  the  DNA  dehybridization  temperature 
(Figure  2B).  When  plotted  as  a  function  of  temperature,  the 


CW  laser-induced  DNA  release  profile  is  consistent  with  the 
thermal  release  profile  with  a  Tm  of  59.1  ±4.5  °C  (Figure  2C). 

Unlike  previous  reports,  we  observed  no  increased  DNA 
release  below  the  dehybridization  temperature  under  CW 
irradiation.  It  is  likely  that  this  discrepancy  may  be  due  to  the 
lower  109  NS/mL  nanoparticle  concentration  used  in  our 
study,  compared  to  higher  1010  NS/mL  concentrations  used  in 
the  other  studies,  where  lower- temperature  DNA  release  was 
reported.7,8, 15,2^  An  earlier  study  reported  by  our  group  on 
collective  photothermal  heating  by  plasmonic  nanoparticles 
showed  that  particle  concentrations  of  1010  NS/mL  result  in 
light  absorption  profiles  with  a  decay  length  of  around  300  jum 
along  the  laser  propagation  direction.  It  was  also  shown  that,  at 
these  high  particle  concentrations,  where  the  scattering  mean 
free  path  is  ~830  //m,  multiple  scattering  of  photons 
significantly  increases  the  fraction  of  incident  light  absorbed. 
These  two  factors,  the  small  heating  volume  and  the  enhanced 
absorption  efficiency  due  to  multiple  scattering,  result  in  large 
temperature  increases  and  gradients  at  the  input  laser  powers 
studied.68  These  temperature  gradients  can  explain  the 
previously  observed  release  of  ssDNA  below  the  dehybridiza¬ 
tion  temperature,  where  the  temperature  was  probed  at  the 
bottom  of  the  temperature  gradient,  not  in  the  region  of 
incident  light.  The  results  we  report  here  clearly  show  that  one 
must  consider  the  full  temperature  profile  of  the  absorbing 
medium,  which  may  quite  possibly  be  significantly  nonuniform, 
in  regions  of  high  nanoparticle  concentration. 

The  similarity  of  the  CW  laser-induced  DNA  release  profile 
to  the  thermal  profile,  as  well  as  the  apparent  dependency  on 
nanoparticle  concentration,  suggests  that  thermal  effects  play  a 
dominant  role  in  DNA  dehybridization,  but  they  do  not 
unambiguously  separate  thermal  and  nonthermal  release 
mechanisms,  such  as  hot-electron  effects.10,69 
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We  performed  a  similar  set  of  CW  light-triggered  release 
studies,  this  time  keeping  the  entire  sample  below  the  DNA 
melting  temperature  by  using  a  thermoelectric  (Peltier  effect) 
stage.  The  NS(o)  dsDNA  solution  was  cooled  to  7  °C  and 
illuminated  from  the  side  with  an  808  nm  CW  laser  at  1.6  and 
3.4  W  for  durations  of  3—24  min.  A  fresh  sample  was  used  for 
each  time  point.  The  temperature  of  the  entire  sample  was 
maintained  below  20  °C  for  all  experiments  (Figure  2D,  inset). 
No  release  was  observed  for  either  laser  power,  indicating  that 
nonthermal  effects  such  as  non-equilibrium  charge  transfer  do 
not  play  a  measurable  role  in  release  (Figure  2D).  A  significant 
fraction  of  hot  electrons  from  plasmon  decay  should  have 
energies  equivalent  to  the  photon  energy  (1.53  eV),  which  is  3 
orders  of  magnitude  higher  than  the  energy  required  for 
dehybridization,  kh  (Tm  —  298.15  K)  ~  3  meV,  that  is,  the 
thermal  energy  difference  between  the  dehybridization  temper¬ 
ature  and  room  temperature.  Even  with  cooling  the  stage  to  7 
°C,  the  energy  requirement  for  dehybridization  is  only  4.6  meV, 
which  is  still  an  insignificant  fraction  of  the  energy  of  the  hot 
electrons.  Therefore,  the  absence  of  release  when  the  entire 
sample  is  maintained  below  room  temperature  supports  that 
thermal  effects  drive  the  DNA  dehybridization. 

Pulsed-Laser-Induced  Release.  The  DNA  release  ob¬ 
served  in  the  CW  laser  experiments  provides  strong  evidence 
that  the  release  mechanism  is  thermal  in  nature  and  requires 
significant  heating  to  dehybridize  DNA.  However,  it  would  be 
far  more  preferable  to  reduce  or  eliminate  this  heating 
requirement  to  avoid  nonspecific  damage  to  surrounding 
tissues  during  DNA  release  in  vivo.  The  combination  of  high 
pulse  energies  and  low  average  powers  achievable  with  pulsed 
lasers  is  a  drastically  different  illumination  regime,  capable  of 
generating  non-equilibrium  heating  effects  local  to  the 
nanoparticle  with  minimal  heating  of  the  surrounding  volume. 
Previous  studies  of  pulsed-laser-induced  DNA  release  have  used 
high-energy  UV  lasers  to  break  the  Au— S  bond51  or  large 
powers  in  the  NIR  to  melt  the  particles  in  order  to  release 
DNA.13,70,71  Here,  we  investigate  the  DNA  release  mechanism 
for  these  NIR  resonant  nanoshell-based  complexes  using  pulsed 
NIR  illumination. 

NS(S)dsDNA  solutions  were  irradiated  for  10  min  with  a 
~160  fs  800  nm  wavelength  pulsed  laser  at  powers  from  2  to  50 
mW  (Figure  3A).  DNA  release  increased  linearly  from  2  to  25 
mW.  However,  during  experimentation,  we  observed  above  25 
mW  there  was  color  change  in  the  solution  indicating  a  loss  of 
nanoshell  integrity  (blue  shaded  region).  We  performed 
experiments  in  triplicate  up  to  this  likely  damage  threshold, 
after  which  we  performed  one  experiment  at  each  laser  power 
to  verify  via  scanning  electron  microscopy  (SEM)  imaging  the 
breakdown  of  nanoshell  integrity.  The  linear  power  depend¬ 
ence  from  2  to  25  mW,  which  differs  distinctively  from  the 
thermal  release  profiles  in  the  CW  case,  indicates  that  the  DNA 
release  under  these  pulsed  excitation  conditions  is  not  likely 
due  to  a  thermal  process.  In  fact,  no  temperature  increase  of  the 
solution  was  observed  (Figure  3B). 

The  expected  transient  temperature  rise  of  the  DNA  due  to 
pulsed  irradiation  at  these  power  levels  is,  in  principle,  sufficient 
to  dehybridize  the  DNA.  To  understand  the  temperature 
experienced  by  the  DNA  during  illumination,  we  employed 
finite-element  method  (FEM,  Comsol;  see  Supporting 
Information)  to  calculate  the  spatial  and  temporal  dependence 
of  the  temperature  field  of  a  single  nanoparticle  in  aqueous 
solution.  We  employed  both  spatial  and  temporal  averages  to 
obtain  average  temperatures  (Figure  S3).  Averaging  the 


Figure  3.  Femtosecond  pulsed-laser-induced  ssDNA  release.  (A) 
DNA  release  profile  where  the  nonthiolated  DNA  strands  are 
fluorescently  tagged,  showing  linear  release  up  to  25  mW,  after 
which  release  decreases.  (B)  Heating  profile  of  femtosecond 
pulsed-laser  DNA  release  at  laser  powers  from  2  to  50  mW  showing 
no  bulk  solution  temperature  increase.  Inset:  Finite-element 
method  simulation  of  the  maximum  temperature  in  a  6.6  nm 
layer  around  the  NS.  (C)  DNA  release  profile  of  the  control  where 
the  thiolated  DNA  strands  are  fluorescently  tagged,  showing  the 
same  release  profile  as  (A),  indicating  Au— S  bond  breakage- 
mediated  DNA  release. 


temperature  in  a  shell  around  the  nanoparticle  with  dimensions 
corresponding  to  the  length  of  the  dsDNA  (6.6  nm)  gives 
average  temperatures  ranging  from  essentially  room  temper¬ 
ature  at  2  mW  to  nominally  140  °C  at  50  mW  (Figure  3B, 
inset).  For  laser  powers  where  considerable  release  is  observed 
(10—25  mW),  the  calculated  temperatures  fall  between  45  and 
80  °C,  which  is  consistent  with  the  Tm  results  from  both  the 
thermal  and  CW  laser  experiments. 

To  investigate  the  nature  of  the  DNA  release  under  pulsed- 
laser  excitation,  we  used  dsDNA  with  a  fluorescein  tag  on  the 
thiolated  DNA  strand.  With  the  thiolated  strand  tagged  with  a 
fluorophore,  if  thermal  or  CW  DNA  dehybridization  occurs,  no 
fluorophore  will  be  released  into  solution.  However,  for  pulsed- 
laser  illumination  conditions,  substantial  fluorescence  was 
observed  with  a  release  profile  very  similar  to  that  shown  in 
Figure  3A,  increasing  linearly  with  increasing  laser  power  up  to 
35  mW  (Figure  3C).  Because  the  binding  energy  of  the  Au— S 
bond  has  been  shown  to  be  ~0.1  eV  in  such  systems,  the 
increased  local  temperature  at  the  nanoparticle  surface  alone  is 
insufficient  to  induce  this  bond  breakage.51  Temperatures  in 
excess  of  890  °C  (khT  =  0.1  eV)  would  be  required  to  thermally 
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break  this  bond,  a  temperature  regime  where  nanoparticle 
melting  and  reshaping  would  also  be  expected  to  occur.  These 
results  strongly  suggest  that  a  nonthermal  mechanism,  such  as 
non-equilibrium  charge  transfer  of  hot  carriers  from  the  metal 
to  the  DNA  molecule,  may  induce  breakage  of  the  Au— S  bond, 
resulting  in  release  of  the  intact  hybridized  dsDNA  molecules. 

The  large  pulse  energies  generated  by  the  femtosecond  laser 
at  the  highest  power  levels  used  in  these  experiments  can 
considerably  increase  the  nanoparticle  temperature  such  that 
nanoparticle  reshaping  occurs.  Irreversible  nanoparticle  reshap¬ 
ing  is  the  reason  for  the  decrease  in  DNA  release  observed 
above  25  mW  of  pulsed  irradiation,  a  conclusion  supported  by 
extinction  spectroscopy  and  SEM  of  the  nanoparticles  before 
and  after  laser  irradiation.  The  extinction  spectra  of  samples 
irradiated  at  0—25  mW  display  a  dipole  peak  centered  at  ~785 
nm  and  a  well-defined  higher-energy  quadrupole  peak  at  ~600 
nm.  From  30  to  50  mW,  the  extinction  maximum  red  shifts  to 
688  nm,  accompanied  by  the  disappearance  of  the  quadrupolar 
plasmon  peak  and  the  appearance  of  a  second  peak  at  ~1015 
nm  (Figure  4A).  These  changes  indicate  loss  of  shell  integrity 


Figure  4.  Nanoparticle  characterization  after  pulsed-laser-induced 
DNA  release.  (A)  Extinction  spectroscopy  showing  a  red  shift  in 
the  plasmon  resonance  at  powers  above  25  mW.  (B)  FEM 
simulation  of  the  average  temperature  of  the  Au  shell  after  a  single 
150  fs  laser  pulse.  Inset:  Temperature  map  for  25  mW  incident 
power  100  ps  after  laser  pulse  (scale  in  °C).  SEM  images  after  laser 
irradiation  at  (C)  2,  (D)  25,  (E)  30,  (F)  35,  and  (G)  50  mW 
showing  particle  degradation  above  25  mW.  Scale  bars  are  200  nm. 


and  particle  aggregation.  SEM  images  show  no  change  in  the 
Au  shells  from  2  to  25  mW  but  reveal  an  increase  in  holes  and 
cracks  in  the  shells  at  30  and  35  mW,  with  widespread  holes  in 
the  shell  at  50  mW  (Figure  4C— G).  These  results  agree  well 
with  previously  reported  reshaping  of  nanoshells  under 
femtosecond  pulsed  irradiation/  The  fact  that  there  is 
considerable  DNA  release  in  the  range  from  2  to  25  mW 
accompanied  by  the  lack  of  Au  shell  breakdown  suggests  that 
the  release  is  not  due  to  reshaping  or  fragmentation  of  the  Au. 

FEM  calculations  also  support  the  observation  of  Au  NS 
reshaping  and  reduced  DNA  release  for  incident  laser  power 
levels  at  and  above  30  mW.  The  average  local  temperature 
increase  of  the  NS  as  a  function  of  incident  laser  power  is 
shown  in  Figure  4B.  At  the  highest  power  of  50  mW,  the  Au 
shell  reaches  calculated  temperatures  around  1000  °C,  close  to 
the  1064  °C  melting  point  of  bulk  Au.  It  is  well-known  that  at 
the  nanoscale  the  energy  densities  required  for  surface 


reshaping  are  a  fraction  of  the  bulk  melting  energy  density.7  ’ 
Our  calculations  show  that,  at  30  mW,  the  Au  shell  reaches  a 
maximum  temperature  of  ~550  °C,  which  appears  to 
correspond  to  the  onset  of  nanoscale  reshaping  for  this 
structure.  Over  the  range  of  laser  powers  used  in  our 
experiments,  the  Au  is  reaching  between  3  and  70%  of  the 
bulk  melting  energy  density  of  195  J/g  (Figure  S4),  with  the 
onset  of  reshaping  occurring  at  ~40%,  the  same  value  obtained 
for  reshaping  of  Au  nanorods.73  These  values  will  vary 
depending  on  nanoparticle  morphology  and  composition  as 
well  as  temporal  pulse  width. 

Time  and  Concentration  Dependence  of  DNA 
Release.  Using  25  mW  of  pulsed  irradiation,  the  time  to 
release  the  maximum  DNA  was  investigated  (Figure  5A). 


Figure  5.  Time  and  concentration  dependence  of  pulsed-laser- 
induced  DNA  release.  Variation  of  release  as  a  function  of  (A) 
irradiation  time  at  25  mW  with  1  X  109  particles/mL  fit  with  an 
exponential  that  reveals  a  time  constant  of  2.2  min  and  (B) 
concentration  of  NS(S)dsDNA  irradiated  at  25  mW  and  10  min 
irradiation  showing  no  appreciable  change  in  release  at  different 
concentrations. 


Approximately  40%  of  release  occurs  after  2  min  of  irradiation 
with  saturation  achieved  after  only  5  min.  Achieving  substantial 
DNA  release  over  short  irradiation  times  with  low  average 
powers  indicates  that  short  pulse-induced  DNA  release  is  a 
highly  attractive  strategy  for  therapeutic  applications. 

Nanoparticle  uptake  in  tumors  is  known  to  vary  widely,74 
indicating  the  desirability  for  a  release  method  that  is  insensitive 
to  nanoparticle  concentration.  For  pulsed-laser  illumination 
periods  of  less  than  10  min  at  25  mW,  no  appreciable  change  in 
the  amount  of  DNA  released  per  NS  was  observed  as 
nanoparticle  concentration  was  increased  (Figure  5B).  The 
relatively  low  average  power  of  the  laser  results  in  negligible 
bulk  heating.  Thus,  DNA  release  induced  by  ultrafast 
irradiation  differs  dramatically  from  CW-induced  release, 
which  depends  on  collective  nanoparticle  heating.  If  ultrafast 
induced  release  arises  from  hot  carrier  transfer,  it  should 
depend  only  on  the  independent  interaction  of  each  nano¬ 
particle  with  the  incident  light  field.  Therefore,  the  magnitude 
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of  release  per  nanoparticle  is  relatively  insensitive  to 
concentration.  At  high  concentrations,  some  dependence  is 
expected  due  to  the  attenuation  of  the  incident  laser  beam; 
however,  as  our  experiments  illustrate,  this  effect  is  not  relevant 
over  the  range  of  nanoparticle  concentrations  relevant  for  in 
vivo  therapeutic  applications. 

CONCLUSION 

In  this  study,  we  showed  that  NIR  CW  light- triggered  DNA 
release  is  the  result  of  thermally  induced  DNA  dehybridization, 
which  induces  release  of  ssDNA.  Collective  heating  effects 
characteristic  of  CW  irradiation  at  the  NS  concentrations 
relevant  for  therapy  make  the  DNA  release  highly  dependent 
on  nanoparticle  concentration.  At  the  low  particle  concen¬ 
trations  achievable  in  tumors,  we  demonstrated  that  the  global 
temperature  must  rise  above  the  DNA  dehybridization 
temperature  to  cause  DNA  release.  In  contrast,  we  observed 
that  a  NIR  ultrafast  pulsed  laser  induces  release  of  dsDNA  with 
no  measurable  bulk  temperature  increase  at  lower  nanoparticle 
concentrations.  For  pulsed-laser  illumination,  the  DNA  is 
released  by  inducing  Au— S  bond  breakage.  With  a  pulsed  laser, 
DNA  release  is  a  single-particle  phenomenon  and  is 
independent  of  particle  concentration.  These  results  clearly 
show  that  pulsed-laser-induced  DNA  release  is  a  highly 
promising  strategy  for  light- triggered  therapeutic  release 
without  inducing  bulk  heating,  making  this  approach  more 
likely  to  preserve  the  integrity  of  healthy  tissues  in  therapeutic 
applications  in  vivo. 

EXPERIMENTAL  METHODS 

Materials.  Potassium  carbonate  anhydrous  was  purchased  from 
Fisher  Scientific.  Hydrogen  tetrachloroaurate(lll)  trihydrate  (HAuC14- 
3H20),  DL-dithiothreitol  (DTT),  and  sodium  chloride  were  purchased 
from  Sigma- Aldrich.  Aminated  silica  cores  (120  nm  diameter)  in 
ethanol  (40  mg/mL)  were  purchased  from  nanoComposix,  Inc.,  San 
Diego,  CA.  HPLC-purified  DNA  strands  and  10  mM  Tris,  pH  7.5,  0.1 
mM  EDTA  buffer  (lx  TE  buffer),  were  purchased  from  Integrated 
DNA  Technologies.  Water  was  deionized  and  filtered  by  a  Milli-Q^ 
water  system  (18.2  MD-cm  at  25  °C,  Millipore). 

DNA  Purification.  Prior  to  nanoshell  functionalization,  thiolated 
ssDNA  [HS-C6H12-TATGATCTGTCACAGCTTGA]  was  incubated 
with  10  mM  DTT  reducing  agent  in  10  mM  TE  buffer  for  1  h  at  room 
temperature.  The  thiolated  DNA  was  purified  with  a  Sephadex  G-25 
column  (NAP  25,  GE  Healthcare),  eluting  with  TE  buffer.  The 
complementary  fluorescein-tagged  ssDNA  [  AT  AC  TAG  AC  AG  T- 
GTCGAACT-Fluor]  was  dispersed  in  600  fiL  of  TE  buffer.  DNA 
concentration  of  both  strands  was  determined  by  measuring  their 
absorbance  at  260  nm.  The  two  DNA  strands  were  mixed  in  a  1:1 
molar  ratio  with  33  mM  NaCl  and  boiled  for  4  min.  The  mixture  was 
cooled  slowly  overnight  to  hybridize  the  strands. 

Nanoshell  Synthesis  and  Functionalization.  Nanoshells  were 
synthesized  by  previously  published  procedures  with  120  nm  diameter 
Si02  cores  and  21  nm  thick  Au  shells.55,57  The  nanoshell  surface  was 
functionalized  by  mixing  dsDNA  with  bare  NS  in  a  ratio  of  20000 
dsDNA  strands/NS  and  stirred  overnight.  The  particles  were  washed 
three  times  by  centrifugation  and  redispersed  in  TE  buffer,  pH  7.5. 
Extinction  spectra  were  obtained  using  a  Cary  5000  UV/vis/NIR 
Varian  spectrophotometer.  ^-Potential  measurements  were  taken  with 
a  Malvern  Zetasizer  Nano  ZS.  TEM  images  were  obtained  on  a  JEOL 
1230  high-contrast  transmission  electron  microscope. 

Circular  Dichroism  Measurements.  CD  spectra  of  5  /aM  dsDNA 
in  33  ftM  NaCl  TE  buffer  (pH  7.5)  were  acquired  in  a  1  cm  path 
length  quartz  cell.  Measurements  were  made  using  a  JASCO  J-810 
spectropolarimeter  equipped  with  a  Peltier-type  temperature  control 
system  at  wavelengths  from  200  to  320  nm.  Thermal  dsDNA 
dehybridization  curves  were  obtained  by  monitoring  the  decrease  in 


ellipticity  with  temperature  in  the  desired  wavelength  range.  The 
temperature  was  varied  between  25  and  90  °C  in  5  °C  increments  with 
an  equilibration  time  of  180  s.  The  melting  CD  spectra  of  dsDNA, 
which  show  structural  transitions  from  the  double-strand  state  to  the 
single-strand  state,  are  shown  in  Figure  SIB.  The  melting  curve  was 
extracted  from  the  change  in  ellipticity  at  246  nm  ((9246),  and  the  first 
derivative  of  a  Boltzmann  fit  determined  the  Tm  to  be  55.2  ±  0.7  °C 
(Figure  SIB). 

Thermal  Studies.  While  being  stirred,  1  mL  of  1  X  109  particles/ 
mL  NS(a)dsDNA  was  heated  to  the  desired  temperature  in  a  1  X  1  X  4 
cm  methacrylate  cuvette  using  a  Peltier  temperature-controlled  cuvette 
stage.  After  20  min  of  heating,  a  300  ftL  aliquot  was  removed  and 
immediately  centrifuged  at  5000g  for  2  min  at  6  °C  to  separate  the 
released  ssDNA  from  the  particles.  The  process  was  repeated  in  5  °C 
increments  up  to  65  °C.  A  fresh  sample  was  used  for  each  temperature 
point.  Released  ssDNA  in  the  supernatant  was  quantified  by 
fluorescence. 

Temperature-Controlled  CW  Laser  Release  Studies.  NS@ 

dsDNA  (2.5  mL  of  1  X  109  p articles /mL)  was  cooled  to  7  °C  using  a 
Peltier  temperature-controlled  cuvette  stage.  While  being  cooled  and 
stirred,  the  NS(a)dsDNA  solution  was  irradiated  with  an  808  nm  CW 
diode  laser  (Diomed,  15Plus,  Angio  Dynamics)  at  1.6  and  3.4  W  with 
a  full  width  at  half-maximum  of  3.5  mm.  The  sample  was  irradiated 
from  the  side  to  minimize  vaporization.  Temperature  was  monitored 
with  an  infrared  thermocouple  (Omega  OS-80 1-MT-K)  with  a  0.5  cm 
spot  size.  After  the  desired  time,  a  300  ftL  aliquot  was  removed  and 
immediately  centrifuged  at  5000g  for  2  min  at  6  °C.  A  fresh  sample 
was  used  for  each  time  point.  Fluorescence  measurements  were  taken 
on  the  supernatant  to  quantify  released  ssDNA. 

Non-temperature-Controlled  CW  Laser  Release  Studies. 
NS(a)dsDNA  (l  mL  of  1  X  109  particles/mL)  was  stirred  on  a 
custom-built  plastic  cuvette  stage.  NS(a)dsDNA  solutions  were 
irradiated  from  the  side  with  the  CW  laser  at  0.4,  0.7,  0.8,  1.0,  and 
1.5  W/cm2.  For  each  laser  power,  the  saturation  temperature  was 
determined,  and  samples  were  then  irradiated  for  0,  2.5,  5,  7.5,  10,  and 
15  min  after  the  saturation  temperature  was  achieved.  Samples  were 
also  taken  2.5  and  5  min  after  the  laser  was  turned  on,  before  the 
saturation  temperature  was  reached.  Temperature  was  monitored  with 
a  thermal  camera  (FLIR  A5).  After  irradiation,  a  300  /aL  aliquot  was 
removed  and  immediately  centrifuged  at  5000g  for  2  min  at  6  °C. 
Spatial  and  temporal  average  temperatures  were  used  as  the 
temperature  value  in  the  analysis.  Fluorescence  was  measured  in  the 
supernatant  of  all  samples. 

Femtosecond  Pulsed-Laser  Release  Studies.  NS(S)dsDNA  (1 
mL  of  1  X  109  particles/mL)  was  stirred  on  a  cuvette  stage  and 
irradiated  from  the  side  with  a  Ti: sapphire  laser  amplifier  (Coherent, 
RegA)  operating  at  250  kHz  with  an  average  power  of  1  W  and  a  pulse 
length  of  ~  150— 180  fs.  The  laser  spot  had  a  full  width  at  half¬ 
maximum  of  44  ft m.  Temperature  was  monitored  with  a  thermal 
camera  (FLIR  A5).  Samples  were  irradiated  for  10  min  at  2—50  mW. 
After  irradiation,  a  300  fiL  aliquot  was  removed  and  immediately 
centrifuged  at  3000g  for  3  min.  Fluorescence  measurements  were 
taken  of  the  supernatants.  SEM  images  were  taken  of  the  particles 
before  and  after  release  (FEI  Quanta  650,  accelerating  voltage  of  25 
kV). 
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controlled  cuvette  stage  (Figure  S2);  detailed  description 
of  FEM  simulations;  FEM  simulation  results  of  spatially 
averaged  temperature  as  a  function  of  time  (Figure  S3); 
and  approximate  energy  density  in  the  Au  shell  as  a 
function  of  laser  power  (Figure  S4)  (PDF) 
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